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Technology  Transfer 

This  research  on  microtextures  has  drawn  the  interest  of  ExxonMobil  Research  and  Engineering  Company  in  Clinton  NJ,  and 
we  have  together  developed  a  microfluidic  chip  to  study  the  process  of  electrocoalescence  of  water  droplets  in  a  crude  oil.  In 
the  electrocoalescence  process,  an  electric  field  is  applied  across  an  emulsion  of  water  droplets  in  oil.  The  field  polarizes  the 
conducting  water  droplets  in  the  insulating  oil,  creating  dipole-dipole  attractive  forces  which  coalesce  the  droplets.  This 
coalescence  is  used  in  the  oil  industry  to  remove  water  droplets  from  the  oil.  The  microfluidic  cell  that  was  designed  uses  flow 
focusing  to  create  a  monodisperse  emulsion  of  the  droplets  in  an  oil  chamber  in  a  transparent  cell,  and  follows  the 
electocoalescence  process  visually,  through  optical  microscopy  in  the  chamber  as  a  field  is  applied  across  the  droplets  by 
electrodes  placed  on  opposite  ends  of  the  chamber.  The  microfluidic  design  allows  a  visual  picture  of  the  electocoalescence 
process  through  the  normally  opaque  crude  oil  because  the  height  of  the  microfluidic  chamber  is  small,  of  order  100  microns, 
and  the  monodisperse  size  distribution  allows  for  a  reproducible  study.  A  model  is  constructed  for  the  coalescence  process, 
which  solves  the  electrostatic  field  equations  in  the  emulsion,  to  compute  the  instantaneous  dipolar  forces  between  the  droplets. 
This  computation  is  done  directly  on  the  configuration  of  the  droplets  in  the  emulsion  as  obtained  from  the  optical  micrographs. 
From  a  map  of  these  forces,  evaluated  on  a  pairwise  basis,  predictions  are  made  as  to  which  droplets  coalesce  based  on 
whether  the  force  is  attractive  (relative  alignment  of  the  droplets  along  the  field)  or  repulsive  (side  alignment  relative  to  the  field), 
and  these  predictions  are  shown  to  agree  with  the  observed  coalescence  process.  The  cell  can  be  used  to  identify  field  strength 
conditions  necessary  to  electrocoalesce  the  water  droplets  in  a  particular  crude,  and  to  identify  surfactant  promotors  for  the 
particular  crude  which  can  accelerate  the  process.  This  research  is  described  in  the  last  chapter  of  the  detailed  progress  report 
which  is  attached.  ? 
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Chapter  1 


Introduction 


1.1  Scope  of  Proposed  Research 

The  overall  goal  of  this  research  program  is  to  study  experimentally  and  theoretically  the 
hydrodynamics  of  fluid  flow  along  microtextured  surfaces.  In  the  context  of  this  research, 
microtextured  surfaces  are  broadly  defined  as  interfaces  with  topological  features  of  char¬ 
acteristic  size  of  the  order  of  hundreds  of  nanometers  to  tens  of  microns.  The  features  in 
general  are  spatially  distributed  along  the  surface  in  either  an  irregular  or  organized  pattern, 
and  can  be  fabricated  from  hard  or  elastomeric  (soft)  materials.  When  fluid  flows  over  a  mi¬ 
crotextured  surface,  the  flow  within  the  texture  can  dramatically  affect  the  hydrodynamics 
above  the  surface  relative  to  the  case  in  which  the  surface  is  smooth,  and  as  such  can  lead 
to  innovative  applications  of  interest  to  the  Department  of  Defense  (DOD). 

The  research  proposal  was  originally  submitted  to  the  Air  Force  Office  of  Scientific  Re¬ 
search  (AFOSR)  under  a  announcement  for  historically  black  colleges  and  universities  and 
minority  serving  institutions  (HDCU /MI)  and  with  Professor  Charles  Maldarelli  of  the  City 
College  of  the  City  University  of  New  York  (CCNY-CUNY,  a  HBCU/MI)  as  PI  and  Profes¬ 
sor  Philip  Hall  of  Imperial  College  London  as  co-PI.  The  research  program  outlined  in  the 
proposal  consisted  of  two  parts:  The  first,  under  the  direction  of  the  PI  with  expertise  in 
materials  fabrication  and  microfluidics,  proposed  fabricating  textured  surfaces  with  a  regu¬ 
lar  pillar  (post)  pattern  using  lithographic  techniques,  and  soft  elastomeric  materials  (soft 
photolithography),  imaging  the  flow  through  the  post  array  with  particle  image  velocime- 
try,  measuring  flow  parameters  such  as  permeability  through  the  surface  with  the  surfaces 
assembled  in  microfluidic  flow  cells,  and  computational  modeling  the  flow  within  and  above 
the  surface.  The  second,  under  the  direction  of  the  co-PI  with  expertise  in  theoretical  hy- 
personics,  examined  the  application  of  these  surfaces  as  passive  control  elements  to  delay 
the  transition  to  turbulence  in  the  near  wall  laminar  boundary  layers  of  aerospace  vehicles 
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in  low  hypersonic  flow.  Hypersonic  boundary  layer  solutions  for  the  flow  over  the  surface 
were  to  be  coupled  to  the  flow  solutions  in  the  microtextures  layer  to  examine  the  effect  of 
the  flow  in  the  texture  on  altering  surface  vortices  and  the  transition  to  turbulence.  Upon 
the  terms  of  the  award  of  the  grant  (May,  2011),  a  three  year  effort,  the  proposal  was  moved 
to  the  Fluid  Dynamics  Program  at  the  Army  Research  Office  (ARO),  and  the  support  of 
the  Co-PI  was  disallowed  as  the  effort  was  intended  to  be  focused  at  the  HBCU/MI  institu¬ 
tion.  In  consultation  with  the  ARO  program  manager  (including  a  site  visit  in  May  2012), 
in  addition  to  the  fabrication  of  the  microtextures,  and  the  study  of  their  hydrodynamics, 
the  application  of  the  textured  surfaces  was  directed  to  the  following  slower  flow  regimes  of 
interest  to  the  ARO: 

1.  Hydrodynamics  of  Water  Droplet  Movement  in  Air  Over  Superhydrophobic  Micropost 
Textures  For  Self-Cleaning 

2.  Hydrodynamics  of  Water  Droplet  Movement  in  Oil  over  Particle  Microtextures  Im¬ 
pregnated  With  Oil  For  Reduced  Droplet  Adhesion 

3.  Hydrodynamics  of  Large  Interfacial  Slip  at  Porous  (Air-Permeable)  Polymer  Surfaces 
For  Drag  Reduction 

4.  Hydrodynamic  Capture  of  Microbeads  in  Microtexture  Obstacle  Courses  for  Biomolec- 
ular  Screening  in  Pathogen  Sensing  and  Point  of  Care  Diagnostics 

The  research  on  these  topics  was  undertaken  for  three  and  one/half  years,  and  concluded 
in  December  of  2014,  involving  four  graduate  students  (two  fully  supported  by  the  research 
and  two  partially  supported),  two  master’s  students  (not  supported  by  this  grant),  two 
undergraduates  and  two  high  school  students  supported  by  funds  supplemented  to  the  grant 
by  the  UA/HSAP  program.  Part  of  the  research  work  on  the  use  of  dielectrophoretic  forces  to 
measure  slip  over  a  polymer  surface  in  a  polymer  surface  has  been  utilized  at  the  ExxonMobil 
Research  and  Engineering  Company  (EMRE)  in  Annandale,  New  Jersey  to  study  the  process 
of  electrocoalescence  of  droplets  in  an  oil  phase.  This  technology  transfer  arose  as  one 
of  the  graduate  students  participating  in  the  ARO  effort  interned  at  ExxonMobil,  and  a 
collaborative  study  was  undertaken. 

This  report  is  organized  as  follows.  In  the  remainder  of  this  Introduction,  a  brief  sum¬ 
mary  is  given  on  the  results  on  each  of  the  major  topics  studied,  along  with  a  report  on  the 
work  on  the  technology  transfer.  The  next  five  chapters  describe  the  research  on  the  hy¬ 
drodynamics  water  droplet  motion  over  superhydrophobic  microtextures  in  air  (Chapter  2), 
under  oil  (Chapter  3),  large  interfacial  slip  over  porous  polymer  surfaces  (Chapter  4),  bead 
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array  obstacle  courses  for  biomolecular  screening  (Chapter  5),  and  the  measurement  of  elec¬ 
trocoalescence  of  water  droplets  in  oil  in  a  microfluidic  cell  collaboration  with  ExxonMobil 
(Chapter  6). 

1.2  Summary  of  Research 

A  brief  summary  of  the  results  on  each  of  the  above  topics,  and  a  report  on  the  technology 
transfer  is  provided  below. 


Hydrodynamics  of  Water  Droplet  Movement  in  Air  Over 
Superhydrophobic  Micropost  Textures  For  Self-Cleaning 

“Superhydrophobic”  surfaces  are  microtextured,  hydrophobic  interfaces  which  do  not 
allow  water  to  penetrate  into  the  texture.  Water  droplets  easily  glide  over  the  texture 
since  the  liquid  rides  over  a  cushion  of  air  in  the  texture.  Contaminant  particles  on  the 
“superhydrophobic”  surface  can  be  removed  from  the  surface  by  the  gliding  motion  of 


Time 


(a) 


(b) 


Figure  1.1:  Droplet  spreading  dynamics  over  a  topography  composed  of  solid  width  =  0.05  and  gap  =  0.05 
(relative  to  the  drop  radius)  for  advancing  contact  angle  of  90°.  (a)  Time  evolution  of  the  droplet  shape  near 
the  contact  line  region.  Droplet  achieves  a  final  equilibrium  shape  at  the  pinning  point  located  at  X  =  0.6. 
(b)  Left  hand  side  plot  shows  slip-stick-jump  motion  of  the  contact  line,  contact  line  slips  on  the  horizontal 
surface,  sticks  after  reaching  the  corner  point  and  jumps  from  one  pillar  to  the  next.  Corresponding  dynamic 
contact  angles  are  plotted  on  the  right. 

droplets  as  they  pick-up  the  particles  and  remove  them  from  the  surface.  To  understand 
this  motion,  hydrodynamic  solutions  were  obtained  for  thetwo  dimensional  spreading  of  a 
droplet  over  a  microtextured  surface  consisting  of  a  line  of  posts,  as  a  model  of  a  self-cleaning, 
superhydrophobic  surface.  The  droplet  is  found  to  execute  a  slip-stick  motion  in  which  it 
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moves  quickly  over  the  tops  of  the  pillars,  but  sticks  at  the  edges  of  the  pillar  before  jumping 
to  the  next  pillar  (Fig.  1.1).  Time-lapse  videos  for  forced  movement  of  droplets  placed  on 
and  moving  over  a  two  dimensional  grooved  surface  fabricated  by  soft  lithography  confirmed 
this  motion  (Fig.  1.2). 


Figure  1.2:  (a)  2-D  anisotropic  parallel  grooved/striped  surfaces  with  width  w  of  175  //in  and  375  //m 
respectively,  g  =  125/nri  and  h  =  100//m  in  both  cases,  (b)  Controlled  spreading  of  a  2-D  droplet  over  a 
striped  surface  with  w  =  175/.mi,  g  =  125/jin  and  h  =  100/.mi.  Flow  rate  through  the  slit  =  10  /.d/miii,  frame 
capture  rate  =  1  frame/sec.  Droplet  shows  pinning  at  the  corner  point  while  its  contact  angle  increases  till 
time  t=31  sec  at  which  time  the  contact  line  jumps  to  the  next  pillar 


Hydrodynamics  of  Water  Droplet  Movement  in  Oil  over  Particle 
Microtextures  Impregnated  With  Oil  For  Reduced  Droplet  Adhesion 

We  have  developed  microtextured,  hydrophobic  surfaces  which,  when  immersed  in  oil,  allow 
water  droplets  to  glide  over  the  surface  without  sticking.  The  textured  surfaces  are  con¬ 
structed  coating  a  nanoparticle/polymer  resin  suspension  on  a  flat  surface  and  curing  the 
resin,  which  is  scalable  to  texturizing  large  surface  areas.  The  textures  are  functionalized  to 
be  very  hydrophobic  so  that  the  oil  saturates  the  texture.  Water  droplets  under  oil  which 
become  attached  to  the  surface  are  nonadhesive  in  that  they  subtend  large  contact  angles  on 
the  surface,  and  slide  easily  along  the  surface  over  a  cushion  of  the  saturating  oil.  We  have 
measured  static  contact  angles  on  these  microtextured  surfaces  constructed  from  colloidally 
which  demonstrate  that  water  droplets  under  oil  on  these  surfaces  retain  large  contact  angles. 
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Water  droplets  drawn  over  the  surface  by  gravity  are  observed  to  move  with  high  mobility 
as  the  contact  angle  hysteresis  is  reduced.  These  surfaces  address  a  longstanding  problem 
of  how  to  prohibit  water  droplets  from  sticking  on  surfaces  under  oils.  Adhered  droplets  are 
the  nucleating  sites  for  corrosion  and  fouling  in  many  applications  involving  oils  on  surfaces, 
and  treatment  with  these  microtextures  can  represent  a  means  for  maintaing  water-droplet 
free  surfaces  under  oil. 


(c)  (d)  (e) 


Figure  1.3:  (a)  A  mixture  of  nanoparticles  and  binder  cures  to  form  a  microtextured  surface  (b)  Atomic 
force  microscopy  image  of  the  particle  surface  in  tapping  mode,  (c)-(e)  A  water  droplet  released  onto  the 
particle  textured  surface  under  oil  glides  over  the  surface 


Hydrodynamics  of  Large  Interfacial  Slip  at  Porous  (Air-Permeable) 
Polymer  Surfaces  For  Drag  Reduction 

Measurements  are  reported  of  “giant”  hydrodynamic  slip  of  the  order  of  one  micron  at  the 
interface  between  a  viscous  mineral  oil  and  the  elastomer  polydimethylsiloxane  (PDMS). 
The  surface  slip  is  quantified  by  a  new,  novel  technique  in  which  a  train  of  monodispersely- 
sized  aqueous  droplets,  approximately  100  microns  in  diameter,  are  formed  in  oil  by  flow 
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focusing,  and  are  streamed  through  a  wide  microfluidic  channel  with  a  height  only  slightly 
larger  than  the  droplets.  The  channel  material  is  PDMS,  as  fabricated  by  molding  over 
a  lithography  manufactured  master.  The  droplets,  heavier  than  the  mineral  oil,  descend 
to  the  bottom  channel  wall.  At  a  downstream  location,  a  uniform  electric  field  is  applied 
in  the  flow  direction  to  separated  pairs  of  droplets,  forcing  the  droplets  towards  each  by  a 
dielectrophoretic  force,  and  the  relative  droplet  velocity  is  measured  by  optical  microscopy. 
The  drag  of  the  droplets  against  the  channel  floor  (along  with  the  droplet-droplet  interaction) 
determines  the  relative  velocity,  and  this  drag  is  affected  by  the  slip  of  the  intervening  oil 
film  along  the  bottom  surface  of  the  channel.  The  distance  between  the  wall  and  the  droplet 
is  determined  from  the  transit  time  of  the  droplets  to  the  point  of  application  of  the  field. 
Measurements  of  the  velocity  are  compared  to  finite  element  hydrodynamic  solutions  for  the 
approach  velocity  to  determine  the  wall  slip  coefficient.  The  large  slip  at  the  PDMS  surface 
suggests  that  the  air  retained  in  PDMS  because  of  its  nanoporous  hydrophobicity  leads  to 
a  low  friction  layer  of  nanobubbles  at  the  interface,  and  this  represents  a  general  means  for 
generating  large  slip  without  (superhydrophobic)  surface  microtexturing.  These  results  are 
particularly  relevant  to  microfluidics  where  PDMS  is  the  standard  material,  and  the  use  of 
oil  streams  with  reagent  water  droplets  have  become  a  dominant  lab  on  a  chip  platforms. 


Figure  1.4:  (a)  Pancake-shaped  air  bubbles  form  at  the  oil/polydimethylsiloxane  interface  because  of  the 
large  air  solubility  of  PDMS  (b)  Flow-focusing  produces  a  train  of  water  droplets  in  oil  (c)  Application  of 
an  electric  field  parallel  to  the  train  creates  a  dielectrophoretic  force  which  causes  the  droplets  to  merge  and 
the  slip  on  the  surface  to  be  measured. 


Hydrodynamic  Capture  of  Microbeads  in  Microtexture  Obstacle  Courses 

FOR  BlOMOLECULAR  SCREENING  IN  PATHOGEN  SENSING  AND  POINT  OF  CARE 

Diagnostics 

We  have  constructed  (Fig.  1.6(a))  a  microtextured  surface  consisting  of  traps  (“V” 
shaped  open  enclosures),  integrated  as  the  bottom  surface  of  a  microfluidic  cell,  to  array 


Figure  1.5:  (a)  Microbeads  hosting  surface  probes  arrayed  in  a  microtexture  of  traps  for  screening  biomolec- 
ular  binding  interactions,  (b)  A  comparison  of  simulation  data  with  experiment  data  for  the  dynamic  response 
of  the  microbead  array.  The  blue  solid  line  with  error  bars  shows  the  normalized  intensity  of  the  fluorescence 
of  NeutrAvidin-FITC  bound  to  0.1  mol  %  NHS-  PEG4-Biotin  microbeads  as  a  function  of  time  for  a  solution 
concentration  of  500  ng/ml  and  a  flow  rate  of  100  pl/min.  The  red  line  shows  the  simulation  fit,  for  a  value 
of  ka  =  180  m3  mole-1  sec-1,  and  the  kinetic  limit  (red)  overlies  the  finite  element  simulation. 


microbeads,  and  demonstrated  the  use  of  the  cell  as  a  device  for  screening  the  binding  inter¬ 
actions  of  a  “target”  biomolecule  against  a  library  of  “probe”  molecules.  Sets  of  microbeads, 
with  each  set  displaying  a  probe  molecule  attached  to  its  surface,  are  streamed  through 
the  cell  sequentially  and  entrapped  by  the  microtextured  surface  to  form  the  probe  library. 
The  screening  step  is  undertaken  by  flowing  an  analyte  target  though  the  cell  and  over 
the  microbeads,  and  identifying  microbeads  which  bind  the  target,  usually  by  fluorescently 
labelling  the  target  and  scanning  for  luminescent  microbeads.  The  hydrodynamic  flow  of 
analyte  through  the  microtexture  is  solved,  and  the  solutions  are  used  to  simulate  the  mass 
transfer  of  the  target  to  the  microbead  surface.  Experiments  are  also  undertaken  using  the 
binding  of  a  target,  fluorescently  labeled  NeutrAvidin,  to  its  binding  partner  biotin,  dis¬ 
played  on  the  microbead  surface,  to  illustrate  the  use  of  the  device  as  a  screening  tool  (Fig. 
1.6(b)). 
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(a)  (b) 


Figure  1.6:  (a)  Schematic  of  the  microfluidic  cell  used  to  study  electrocoalescence,  and  (b)  a  dispersion  of 
monodisperse  droplets  of  water  in  a  crude  oil. 


Electrocoalescence  of  Water-in- Crude  Oil  Emulsions  in  Two  Dimensions: 
Technology  Transfer  with  ExxonMobil 

This  study  reports  on  the  technology  transfer  on  the  research  of  this  proposal  with  Exxon¬ 
Mobil.  Using  the  microfluidic  flow  studies  of  Chapter  4  on  dielectrophoretic  coalescence  of 
drops,  we  have  examined  the  electrocoalescence  of  water  droplets  in  crude  oils  using  mi¬ 
crofluidics,  In  the  electrocoalescence  process,  an  electric  field  is  applied  across  an  emulsion 
of  water  droplets  in  oil.  The  field  polarizes  the  conducting  water  droplets  in  the  insulating 
oil,  creating  dipole-dipole  attractive  forces  which  coalesce  the  droplets.  This  coalescence  is 
used  in  the  oil  industry  to  remove  water  droplets  from  the  oil.  The  microfluidic  cell  that 
was  designed  uses  flow  focusing  to  create  a  monodisperse  emulsion  of  the  droplets  in  an  oil 
chamber  in  a  transparent  cell,  and  follows  the  electocoalescence  process  visually,  through 
optical  microscopy  in  the  chamber  as  a  field  is  applied  across  the  droplets  by  electrodes 
placed  on  opposite  ends  of 

the  chamber.  The  microfluidic  design  allows  a  visual  picture  of  the  electocoalescence 
process  through  the  normally  opaque  crude  oil  because  the  height  of  the  microfluidic  chamber 
is  small,  of  order  100  /urn,  and  the  monodisperse  size  distribution  allows  for  a  reproducible 
study.  A  model  is  constructed  for  the  coalescence  process,  which  solves  the  electrostatic  field 
equations  in  the  emulsion,  to  compute  the  instantaneous  dipolar  forces  between  the  droplets. 
This  computation  is  done  directly  on  the  configuration  of  the  droplets  in  the  emulsion  as 
obtained  from  the  optical  micrographs.  From  a  map  of  these  forces,  evaluated  on  a  pairwise 
basis,  predictions  are  made  as  to  which  droplets  coalesce  based  on  whether  the  force  is 
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attractive  (relative  alignment  of  the  droplets  along  the  held)  or  repulsive  (side  alignment 
relative  to  the  held) ,  and  these  predictions  are  shown  to  agree  with  the  observed  coalescence 
process  (Fig.  1.6). 
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Chapter  2 


Hydrodynamics  of  Water  Droplet 
Movement  in  Air  Over 
Superhydrophobic  Post  Microtextures 
For  Self-Cleaning 

2.1  Background 

A  microtextured  interface  presents  topological  features  of  elevations  and  depressions  with 
length  scales  of  the  order  of  lOOnm  -  100/im  in  height  and  width.  When  the  surface  of  the 
features  is  hydrophobic  (nonpolar),  water  is  unable  to  penetrate  into  the  topology.  From  a 
materials  viewpoint,  surface  textures  which  are  hydrophobic  can  be  fabricated  by  either  tex¬ 
turing  (lithographically,  etching,  etc.)  nonpolar  materials  (e.g.  hydrocarbon  or  fluorocarbon 
(teflon)  polymers),  layering  hydrophobic  particles  or  depositing  polymers  on  flat  surfaces,  or 
coating  already  formed  textures  with  nonpolar  molecules  (e.g.  methyl-  or  fluoro-silanes)  or 
polymers.  The  inability  of  water  to  penetrate  into  the  texture  leads  to  striking  static  and 
dynamic  interfacial  effects  when  an  aqueous  drop  is  placed  on  the  surface:  As  shown  in  Fig. 
2.1(a),  when  a  water  droplet  is  placed  on  a  smooth  flat  surface  either  made  of  either  the 
same  material,  or  having  the  same  surface  coat  as  the  microtexture,  the  droplet  moves  over 
the  solid  surface  and  comes  to  equilibrium.  The  contact  angle  6  at  equilibrium  along  the 
contact  line  is  large  (90°  or  larger)  because  of  the  surface  hydrophobicity.  For  very  hydropho¬ 
bic  materials  or  coatings  (e.g.  teflon),  6  can  be  as  large  as  120°  (but  no  larger)  and  drops 
appear  to  bead-up  on  the  surface.  When  a  water  droplet  is  placed  upon  and  spreads  along  a 
hydrophobic  microtexture  in  air  (Fig.  2.1(a)),  the  droplet  moves  along  elevations  and  above 
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Figure  2.1:  (a)  Equilibrium  configurations  and  contact  angles  of  an  aqueous  drop  on  a  flat,  hydrophobic 
surface  and  a  microtextured  interface  with  a  hydrophobic  surface  in  which  water  does  not  wet  into  the 
texture  (Cassie-Baxter  wetting),  (b)  Contact  angle  hysteresis  on  a  textured  surface  in  which  water  wets  into 
a  texture  not  as  hydrophobic  or  hydrophilic  (Wenzel  state),  and  gliding  motion  without  hysteresis  over  a 
hydrophobic  microtexture  in  the  Cassie-Baxter  state. 
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pockets  of  air  in  the  microtexture,  rather  than  saturating  into  the  texture,  and  stops  as  the 
contact  line  becomes  pinned  on  the  edges  of  the  texture.  Because  of  the  pinning,  for  droplets 
of  equal  volume  as  on  the  solid  surface,  the  static  equilibrium  contact  angle  on  the  texture 
is  much  larger  than  that  over  the  flat  surface  of  the  same  hydrophobic  material  or  surface 
coat,  and  the  droplet  appears  to  be  nearly  non-adherent,  resembling  pearls  or  marbles  on 
a  surface,  as  contact  angles  0*  approaching  180  °  can  be  achieved  (Cassie-Baxter  wetting 
and  super- hydrophobic  behavior)1-19.  Droplet  dynamics  is  also  significantly  different  along 
hydrophobic  microtextured  surfaces  relative  to  flat  non-polar  surfaces:  On  the  flat  surface, 
the  drop  drags  along  the  surface  with  large  surface  friction  to  the  liquid  viscosity,  while  along 
the  microtextured  surface  the  air  cushion  significantly  reduces  the  surface  skin  friction  and 
droplets  glide  along  the  surface.  Droplet  movement  is  accompanied  by  the  motion  of  the 
advancing  and  receding  contact  lines  at  the  leading  and  trailing  edge  of  the  droplet  surface. 
On  an  ideally  flat  surface  which  is  perfectly  smooth,  the  contact  line  advances  with  the 
equilibrium  contact  angle  on  both  edges,  the  air/liquid  interfacial  tension  force  (7  on  the 
droplet  is  equal  on  both  edges  and  there  is  no  net  contact  line  resistance  to  the  motion. 
But  if  the  surface  is  textured,  and  the  texture  is  not  very  hydrophobic  (or  is  hydrophilic) 
so  that  the  aqueous  phase  wets  into  the  texture  (Wenzel  wetting),  contact  angles  9wettexture 
are  much  smaller  than  over  the  superhydrophobic  surface,  but  the  contact  lines  at  the  front 
and  back  of  the  drop  impede  the  droplet  motion  (Fig.  2.1(b)).  For  the  droplet  to  move, 
at  the  leading  edge  large  contact  angles  (>  9)  are  necessary  (relative  to  the  mean  plane)  to 
move  down  an  elevation,  while  small  angles  (6)  are  necessary  at  the  back  end,  leading  to  a 
net  air /water  interfacial  tension  force  on  the  droplet  retarding  the  forward  motion  (contact 
angle  hysteresis)20-22.  This  resistance  immobilizes  droplets  on  tilted  surfaces  from  moving 
downward  due  to  gravity,  and  large  tilt  angles  can  be  necessary  for  drop  motion.  Water 
droplets  on  superhydrophobic  microtextured  surfaces  exhibit  a  reduced  resistance  to  contact 
line  motion  as  the  liquid  perimeter  of  the  drop  either  at  the  forward  or  trailing  edges  of  the 
drop  jumps  from  one  pinning  point  on  the  texture  to  the  next  (Fig.  2.1(b).  The  advancing 
and  receding  contact  angles  are  equal  to  9*,  and  there  is  no  resultant  interfacial  tension 
traction  force23-28.  Tilt  angles  of  only  a  few  degrees  are  necessary  to  make  drops  move  over 
inclined  superhydrophobic  surfaces. 

Nature  provides  several  examples  of  water-repellant,  superhydrophobic  surfaces,  among 
them  the  lotus  leaf29,30  (Fig.  2.2(a))  and  the  legs  of  a  water  strider31  (Fig.  2.2(c)),  where 
scanning  electron  microscopy  (SEM)  examinations  detail  a  microstructure  with  waxy  surface 
coats  which  create  superhydrophobicity  (2.2(b)  and  2.2(d)).  In  the  case  of  the  water  strider, 
the  high  contact  angles  due  to  the  pinning  allow  the  strider  to  support  its  weight  on  water. 
In  the  case  of  the  lotus  leaf,  droplets  glide  over  inclined  leaf  surfaces,  collecting  particles  on 
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Figure  2.2:  Examples  of  superhydrophobic  surfaces  in  Nature:  (a)  the  lotus  leaf  and  and 
SEM  of  its  hierarchical  microtexture  of  posts  ( papillae )  with  nanoscale  waxy  tubules  coating 
the  projections  and  (b)  the  non-wetting  legs  of  a  water  strider  and  an  SEM  of  its  elongated 
microstructures  (microsetae)  22 
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the  surface  and  allowing  the  lotus  leaf  to  be  self-cleaning  (2.3(a))  in  removing  dirt  that  can 
block  sunlight  and  photosynthesis.  Note  that  on  flat  surfaces  in  which  liquid  droplets  drag 
over  the  surface,  particles  are  redistributed  by  the  droplet  and  do  not  self-clean  (2.3(a)).  The 
self-cleaning  property  of  superhydrophobic  surfaces  can  be  used  in  many  technologies  such 
as  self-cleaning  coatings  for  satellite  dishes,  solar  energy  panels,  photovoltaics  and  textiles32. 
As  mentioned  above,  the  reduced  skin  friction  when  a  liquid  flows  over  the  surface  suggests 
applications  as  drag  reduction  coatings  for  open  and  closed  channels  (pipeline)  flows 31,33-35 
(2.3(b)),  and  water  vehicles  such  as  ships  and  submersibles.  In  addition,  because  water  at 
the  interface  with  the  microtexture  only  minimally  contacts  the  solid  part  of  the  texture 
at  the  elevations,  these  microtextures  are  effective  as  interfaces  which  resist  ice  formation 
(icephobic  surfaces)  and  the  adsorption  of  biomolecules  and  retention  of  microorganisms  such 
as  algae  (biofouling) 36 
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Figure  2.3:  Applications  of  superhydrophobic  surfaces:  (a)  Self-cleaning  ability  used  by  the 
lotus  leaf.  Particulates  on  the  flat  surface  (c)  are  redistributed  by  a  sliding  drop,  while  (d) 
on  a  textured  surface  droplet  rolls  as  opposed  to  slide  and  collects  the  particulates  on  its 
way.  (b)  Drag  reduction:  Air  cushion  on  the  textured  surface  of  a  wall  reduces  the  fluid  drag 
exerted  on  the  wall  since  the  fluid  interface  has  zero  shear  stress.37 


2.2  Research  Objectives 

Most  of  the  research  on  superhydrophobicity  has  centered  on  the  calculation  of  the  equi¬ 
librium  states  of  droplets  resting  on  a  surface.  The  initial  studies  of  Cassie  and  Baxter38 
modeled  the  surface  topology  with  elevations  all  of  the  same  height  with  flat  tops  separated 
by  air  gaps  and  a  solid  area  fraction  given  by  os  (see  Fig.  2.4(a)).  A  variational  approach 
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was  used  to  compute  the  differential  change  in  interfacial  energy  over  this  idealized  surface 
as  the  drop  contact  line  spreads  over  an  area  5a  along  the  solid  tops  and  across  the  air  gaps, 
with  5a  assumed  to  be  much  larger  than  the  sizes  of  the  solid  areas  and  the  air  gaps.  Thus 
the  contact  line  samples  several  microstructures,  and  the  following  macroscopic  or  “average” 
equilibrium  contact  angle  9*  is  obtained: 

COS#*  =  —1  +  0s(cos  9 planar  +  1)  =  - - -  (2-1) 

7 

where  7 sv  is  surface  energy  of  the  dry  solid,  jsl  is  the  surface  energy  of  the  wetted  solid,  7 
is  the  surface  tension  of  the  liquid  and  9pianar  is  the  (obtuse)  contact  angle  the  water  makes 
on  a  planar  surface  of  the  (hydrophobic)  texture  material  or  a  planar  surface  with  the  same 
surface  coat.  Since  the  surface  is  hydrophobic  7 sv  is  less  than  7 sl  and  the  contact  angles  on 
the  microtextured  surface  are  obtuse  and  approach  180°  as  the  area  fraction  tends  to  zero. 
If  the  aqueous  phase  wets  the  topology  (Wenzel  wetting),  i.e.  the  surface  is  hydrophilic  or 
slightly  hydrophobic,  then  using  a  similar  derivation39  (see  Fig.  2.4(b))  the  contact  angle  is 
given  by 

a*  _  a  _  7sv  “  7 SL  (n  oN 

cos  0wcttcd  r  •  cos  0pianar  (2-2) 

where  r  is  the  surface  roughness,  defined  as  r  =  (actual  surface  area) /(geometric  surface 


Figure  2.4:  Variational  approach  for  obtaining  equilibrium  contact  angle  for  a  water  droplet 
spreading  over  (a)  a  microtextured  surface  in  which  the  water  does  not  wet  the  texture 
(Cassie-Baxter  superhydrophobic  regime)  and  (b)  the  water  wets  the  texture  (Wenzel  state). 


area)  and  6  is  the  (acute)  contact  angle  of  water  on  the  planar  hydrophilic  surface.  Since  the 
water  wets  the  solid,  7 sv  is  greater  than  7 sl  and  the  angles  on  the  wetted  texture  is  acute 
and  decrease  as  the  roughness  increases  (see  the  discussion  in40  42 .  More  recent  research 
on  the  calculation  of  the  equilibrium  states  have  obtained  global  free  energy  minima  for  a 
droplet  on  a  given  surface  topography  without  averaging  the  texture  properties.  In  these 
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procedures,  the  free  energy  of  a  given  volume  of  a  droplet  in  the  form  of  a  spherical-cap 
(gravity  is  ingnored)  is  calculated  for  different  contact  angles  corresponding  to  different 
wettiing  surface  coverage.  A  global  minimum  is  obtained  in  the  free  energy,  corresponding 
tothe  equilibrium  contact  angle,  and  this  angle  was  shown  to  be  in  close  agreement  with 
the  predicted  values  from  Wenzel  and  Cassie-Baxter  equations3,43-45.  The  energetics  of  the 
transitions  between  the  Wenzel  wetting  and  Cassie-Baxter  nonwetting  states  have  also  been 
studied  both  experimentally  as  well  as  theoretically7,43,46-48.  Again  free-energy  minimization 
principles  are  used  to  describe  which  wetting  state  is  favorable  and  phase  diagrams  are 
constructed  based  on  the  surface  geometric  properties. 

Although  the  above  theories  based  on  thermodynamic  grounds  do  elucidate  the  effect 
of  the  geometric  properties  of  the  texture  on  the  achieved  equilibrium  contact  angles,  the 
dynamics  of  the  spreading  process  and  how  the  texture  affects  the  rate  of  spreading  -  all 
important  for  applications  in  self-cleaning  -  have  not  been  studied  in  detail.  The  aim  of 
this  study  is  to  examine  theoretically  and  experimentally  the  hydrodynamics  of  the  spread¬ 
ing  process  as  a  water  droplet  spreads  along  a  microtextured  surface  in  superhydrophobic 
(Cassie-Baxter)  wetting,  and  focus  on  the  interaction  of  the  fluid  with  the  surface  topology. 
Only  a  few  articles  have  appeared  in  the  literature  studying  the  dynamics,  and  these  have 
focused  on  theory  and  have  used  the  Lattice-Boltzmann  method  to  model  the  spreading48-51. 
Our  modeling  effort  will  focus  on  droplet  motions  in  which  the  fluid  inertia  is  negligible,  and 
we  will  use  a  boundary  integral  method  to  solve  the  hydrodynamic  equations  for  the  gravity- 
driven  motion  of  the  fluid  in  the  drop.  Our  simulations  are  two-dimensional,  and  model  the 
surface  topology  as  a  regular  pattern  of  microposts  as  in  Fig.  2.4(a).  In  Section  2.3,  the 
formulation  of  the  problem  is  presented.  Section  2.4  describes  the  numerical  solution  to  the 
problem  where  we  demonstrate  the  use  of  Biharmonic  Boundary  Integral  Method  along  with 
the  code  validation,  followed  by  numerical  results  and  discussion  in  section  2.5.  Following 
these  calculations,  we  present  a  series  of  related  experiments  on  visualizing  the  two  dimen¬ 
sional  motion  of  droplets  across  a  topology  of  parallel  grooves  that  are  fabricated  using  soft 
lithography.  The  flow  is  visualized  using  photography  and  a  high  speed  camera.  In  section 
2.6.1  we  review  the  soft  lithography  methods  used  for  fabricating  regular /periodic  micron 
scale  patterns  on  a  surface.  To  demonstrate  in  general  how  surfaces  with  regularly  arrayed 
features  can  be  superhydrophobic,  in  the  following  section  2.6.2  we  fabricate  surfaces  with 
two  dimensional  arrays  of  pillars  and  holes,  and  measure  the  equilibrium  contact  angles  of 
aqueous  droplets  placed  on  the  surface  by  imaging  the  droplets  on  the  surfaces.  This  imaging 
also  reveals  the  Cassie-Baxter  wetting  conditions  on  the  surface.  The  measured  contact  an¬ 
gles  are  then  compared  with  the  Cassie-Baxter  equation  (2.1).  In  the  next  section  2.6.3,  the 
experiments  on  the  grooved  surfaces  are  described,  and  compared  to  the  results  of  the  two 
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dimensional  simulations  of  the  previous  section.  Experimentally  droplet  dynamics  is  stud¬ 
ied  by  two  procedures:  firstly  under  gravity  driven  spreading  and  secondly  with  controlled 
spreading,  i.e.  by  increasing  the  volume  of  the  droplet  gradually  with  time. 

2.3  Model  Formulation 

We  study  the  gravity  driven  spreading  motion  of  a  two-dimensional  aqueous  droplet  on  a 
super  hydrophobic  surface  in  the  inertialess  Stokes  flow  regime.  The  time  evolution  of  the 
droplet  profile  can  be  obtained  by  solving  the  Stokes  equations  and  computing  the  drop 
shape  quasi-statically.  Consider  a  fluid  droplet  on  a  surface  with  a  viscous  fluid  region  (Q) 
separated  from  an  inviscid  air  overlayer  by  a  boundary  (dfl).  The  Stokes  equations  obeyed 
by  the  fluid  are, 

VP  =  /rV2u  +  pg 
V  •  u  =  0 

where  u(x,y)  is  the  fluid  velocity,  P(x, y)  is  the  pressure,  p  and  ji  are  the  density  and  the 
viscosity  of  the  fluid,  and  g  is  the  gravitational  acceleration.  Boundary  conditions  on  the 
air /fluid  interface  are  normal  and  tangential  stress  balances.  The  inviscid  phase  is  assumed 
to  be  passive,  hence  zero  tangential  stress  is  exerted  on  the  droplet  surface,  and  the  normal 
stress  exerted  by  the  droplet  liquid  is  balanced  by  the  surface  tension  (Laplace  pressure)  and 
the  constant  pressure  in  the  inviscid  phase. 


t  •  a  •  n  =  0 

(2.4) 

nan  —  —7  k 

(2.5) 

where  7  is  the  surface  tension  of  the  liquid  and  k,  is  the  curvature  defined  as  positive  if  the 
center  lies  within  the  liquid,  t  and  n  are  counter-clockwise  tangential  and  outward  normal 
unit  vectors  and  a  is  the  stress  tensor 

<x  =  -PI  +  p  (Vu  +  ( Vu)T)  . 

On  the  solid  wall  the  classical  no-slip  boundary  condition  is  applied,  n(x,y)  =  0. 

2.3.1  Singularity  at  the  contact  point  and  slip  velocity 

As  the  droplet  spreads,  the  contact  line  has  to  move  or  slide  along  the  surface.  One  common 
problem  associated  with  the  hydrodynamic  solution  of  problems  with  a  contact  line  move- 


(2.3) 
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ment  is  the  singularity  that  arises  due  to  the  discontinuity  in  the  boundary  conditions  at 
the  contact  point.  On  the  free  air/water  interface  we  have  a  finite  velocity  but  on  the  solid 
wall  the  fluid  velocity  is  zero  due  to  the  no-slip  condition.  This  gives  rise  to  an  unrealistic 
infinite  force  at  the  contact  line52.  As  a  consequence,  although  the  solid  surface  remains 
at  rest,  fluid  in  contact  with  it  has  to  move  violating  the  no-slip  boundary  condition.  For 
highly  wetting  fluids,  to  eliminate  this  singularity,  a  precursor  film  models  have  been  used 
widely,  where  it  is  assumed  that  a  thin  layer  of  liquid  (precursor  film)  is  present  ahead  of  the 
contact  line  and  the  free  fluid-fluid  interface  then  merges  smoothly  with  the  precursor  film 
at  the  contact  line53,54.  For  partially  wetting  fluids,  this  stress  singularity  can  be  avoided 
by  relaxing  traditional  no-slip  boundary  condition,  i.e.  by  allowing  some  finite  slip  velocity 
at  the  contact  line  or  by  not-so-realistic  approach  of  setting  the  contact  angle  to  7r  (rolling 
motion  instead  of  sliding)55.  Slip  velocity  at  the  contact  line  can  be  imposed  in  several 
ways.  Navier-type  slip  model  is  based  on  the  concept  of  finite  slip-length,  b,  (distance  be¬ 
yond  liquid-solid  interface  in  a  direction  perpendicular  to  the  wall  where  velocity  approaches 
zero),  where  slip  velocity  at  a  wall  is  proportional  to  the  wall  shear  stress  with  propor¬ 
tionality  constant  as  b.  Greenspan  used  Navier-type  slip  model  near  the  contact  line  while 
smoothly  reverting  to  no- slip  condition  outside  the  contact  line  region56.  More  simplified 
approach  is  using  a  numerical  slip  as  demonstrated  by  Mazouchi  et  al57, where  slip  at  the 
contact  point  is  implemented  implicitly  by  discretizing  domain  with  grid  points  around  the 
contact  point  but  not  including  the  contact  point.  The  governing  equations  are  then  solved 
with  traditional  no-slip  condition  at  all  grid  points  along  the  solid  substrate  and  the  position 
of  the  new  contact  line  is  determined  by  extending  the  calculated  fluid-fluid  interface  till  the 
solid  surface  with  independently  prescribed  contact  angle.  In  this  work  Mazouchi  et  al  have 
assumed  the  prescribed  contact  angle  is  constant  during  the  simulations  and  is  equal  to  the 
static  advancing  contact  angle,  independent  of  the  contact  line  velocity. 

We  use  a  slightly  modified  slip  model  than  Mazouchi  et  al,  where  we  compute  slip  velocity 
by  taking  into  account  dynamic  contact  angles  and  contact  angle  hysteresis.  Contact  angle 
hysteresis  is  a  phenomena  during  which  contact  line  remains  fixed  at  a  position  (zero  slip 
velocity)  for  a  range  of  contact  angles,  9rec  <  9  <  9adv,  where  0rec  and  9adv  are  the  static 
receding  and  static  advancing  contact  angles  respectively.  When  there  is  a  relative  motion 
between  a  contact  line  and  solid  surface,  the  contact  angle  liquid  makes  with  the  solid,  6dyn , 
is  a  function  of  the  relative  velocity.  This  contact  angle  hysteresis  and  a  relationship  between 
dynamic  contact  angle  and  slip  velocity  is  shown  in  Fig.  2.5.  Only  when  the  dynamic  contact 
angle  9dyn  becomes  larger  than  0a()v ,  droplet  interface  slips  forward  with  a  velocity  V  >  0  and 
similar  way  when  9dyn  <  9rec ,  the  interface  recedes,  i.e.  V  <  0.  There  is  a  body  of  literature 
where  experimentally  measured  dynamic  contact  angles  for  different  slip  velocities  are  fitted 
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to  an  empirical  relation,  one  such  relationship  is  the  Hoffman- Jiang  correlation58: 


cos 6adv  c°s^  =  tanh(4.96C'a0'702)  (2.6) 

cos  0adv  T 1 

where  Ca  is  a  capillary  number,  Ca  =  fjV/'y,  V  is  the  slip  velocity.  We  use  numerical  slip 
model  where  new  contact  line  position  is  obtained  iteratively  till  above  relation  is  satisfied 
for  the  contact  angle  and  slip  velocity. 


Figure  2.5: 


Relationship  between  dynamic  contact  angle  and  slip  velocity. 


2.4  Numerical  Solution 

2.4.1  Biharmonic  Boundary  Integral  Method 

We  use  the  Boundary  Integral  Method  for  solving  this  problem.5960.  The  advantage  of  using 
the  Boundary  Integral  Method  is  that  we  reduce  the  dimensionality  of  the  problem  by  one, 
by  making  use  of  Green’s  second  identity;  i.e.  transforming  the  Stokes  equations  which  are 
valid  on  the  entire  fluid  domain  to  integral  equations  which  contain  only  information  from 
the  boundaries.  It  is  convenient  to  use,  as  flow  variables,  the  streamfunction  T  and  vorticity 
uj  defined  as  follows. 

n x  ~n  i  y 

oy  ox 

u j  =  —V  x  u 
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(2.7) 

(2.8) 


Then  the  Stokes  equations  in  terms  of  these  new  flow  variables  can  be  written  as  a  bihar¬ 
monic, 

V4T  =  0 


or  in  split  form, 

V2T  =  u,  V2lj  =  0  (2.9) 

Now  applying  Green’s  identity  to  these  biharmonics,  the  streamfunction  and  vorticity  at 
any  arbitrary  point  p  can  be  written  in  terms  of  information  only  from  the  boundary  points 
q  E  <9Q. 

V (p)  '  *(p)  =  J  (*(Q)  '  ^  •  dti(q) 

dn 

+i  /  (w<?) '  ^  '  °2)  ' dm  (2'10) 

dn 

v(p)  ■  <*>(p)  =  J  (u(q)  ■  •  G^j  ■  dh(q)  (2.11) 

an 

where  Q(q)  is  a  differential  increment  along  the  boundary  at  q,  r/(p)  and  is  given  by 

r](p)  =  0,  if  p  £  0  +  dfl, 

=  tp,  if  p  E  dQ,  (2-12) 

=  27 r,  if  p  E  fl. 

and  where,  <p  is  the  angle  between  the  two  tangents  on  either  side  of  the  point  p.  G\  and  Go 
are  the  fundamental  solutions  of 

V2Gi  =5(\p-q\),  and  V4G2  =  5  (\p  -  q\) 

which  gives, 

Gi  =  log  \p  —  q\  ,  and  G2  =  \p  -  q\2  (log  \p  -  q\  -  1)  (2.13) 

We  have  four  variables  ip,  u  and  their  normal  derivatives  involved  in  the  governing  equa¬ 
tions  (2.10)  and  (2.11).  Thus  we  have  to  write  boundary  conditions  in  terms  of  these  four 
variables.  Rigorous  manipulation  starting  from  Eq.  (5.5)  yields  the  tangential  stress  balance 
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Figure  2.6:  Schematic  of  approximating  droplet  with  a  polygon. 


and  eliminating  the  pressure  term  from  equations  (2.3)  and  (5.6)  we  get  the  normal  stress 
balance  condition  as, 

duj  d2  ( d d'  \  d  (  <9T  \  7  8k  p 

dn  ds2  \dn )  ds  \  ds  J  pds  p^y) 

where  s  is  the  arc-length.  The  no-slip  condition  at  the  solid  surface  can  be  written  as  T 
=  0  and  |^  =  0.  These  boundary  conditions  and  the  two  biharmonic  boundary  integral 
equations  are  need  to  be  solved  simultaneously  to  obtain  T  and  uj  for  the  boundary  points. 

2.4.2  Implementation  of  Boundary  Integral  Method 

As  an  analytical  solution  of  this  problem  is  not  feasible,  a  numerical  approach  has  to  be 
adopted.  We  approximate  the  drop  boundary  by  a  polygon  with  N  linear  segments  such 
that  flow  variables  have  constant  values  over  each  segment.  As  mentioned  earlier,  we  employ 
a  modified  numerical  slip  model  to  solve  this  problem.  In  this  model  we  do  not  evaluate  flow 
variables  exactly  at  the  contact-line  but  at  the  midpoints  of  each  segments  as  shown  in  the 
Fig.  (2.6),  where  q/s  are  the  vertices  of  the  polygon  and  midpoints  are  the  collocation 
points  or  nodes,  where  we  evaluate  the  flow  variables.  Hence,  we  apply  the  stress  balance 
conditions  on  every  node  at  the  air/water  interface  and  the  no-slip  condition  is  applied  to 
the  nodes  which  are  on  the  solid-surface.  By  avoiding  solving  equations  at  the  contact-point 
we  naturally  avoid  the  singularity. 
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The  governing  equations  in  discrete  form  can  be  written  as, 


N  (  r 

vi*i  =  Y,  ^  J 

j=1  \  dQj 


dGlr^~  d*i  f  ^  0~  ,  1  f  dG2n~  1  duj  f  „  ^ 
^ — 7) —  /  Gid7l-\--ujj  /  -x — 371--— —  /  G237l 
on  on  4  on  4  on 


\ 


dQi 


dQi 


N 

V iUJi  =  ^2 
3= 1 


/ 


LJi 


dG'-dtt  -  f  GidQ 


\ 


V 


dn 


dn 


dQn 


8Qn 


) 


which  can  be  written  in  matrix  form  as 


where, 


Aij  — 


Bij 


3 T  dcu 

AT  +  B  — — b  Clu  +  D—  =  0 
3n  3n 

Aw  +  B^  =  0 
3n 


3 

—  log  | pi  -  q |  3Q,(q)  -  rjiSij, 


qGdQj 


q£dO,j 


Cii  =  X 


log  \pi-q\3fl(q), 


3 


0^  [I  Pi  -  q\2  (log  I  Pi  -q\-  1)]  3n(q), 


Dij  — 


I  Pi  -  q\2  (log  |  Pi  -q\-l)  dn(q). 


q^dQj 


(2.16) 


) 


(2.17) 


(2.18) 

(2.19) 


(2.20) 

(2.21) 

(2.22) 

(2.23) 


Taking  advantage  of  the  polygonal  geometry,  we  evaluate  these  integrals  analytically  follow¬ 
ing  Kelmanson’s  approach  Consider  a  triangle  formed  by  node  p  and  endpoints  ./,  and  ./,+J 
of  any  boundary  dilp  as  shown  in  the  Fig.  (2.7)  where, 


o>=\Pi-Qj  h  b=\pi-qj+1\,  h  qj  -  r/,- .  ,  | . 

P  =  ^PiOjQj+u  and  C  =  ^qjPiQj+i- 

Using  these  notations,  we  can  evaluate  integrals  in  Eqs.  (2.20)-(2.23)  exactly  as, 

J  ^  log  | Pi  ~  q\  3&(q)  =  C, 

qGdflj 


(2.24) 


(2.25) 
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Figure  2.7:  Triangle  geometry  for  evaluating  integrals  in  eqns.  (2.20)-(2.23). 


J  log  | pi  —  q\  dQ(q)  =  a(log  a  —  log  b )  cos  (3  +  h  log  b  —  h  +  a(  sin  (3  =  1,  say,  (2.26) 

qEdQj 

j  iL  ^pi  - q^2  ^iog  \pi  -q\-  ^(9) = a(2/  - k)  sin/^>  (2-27) 

qGdflj 

J  [\Pi  -  q\2  (log  \Pi  -q\-  1)]  dQ(q)  =  (asin/3)2  j  I -^h-  ^ a(sm/3 

(2.28) 

loga-i)}. 

Also  we  discretize  the  tangential  and  normal  stress  balances,  Eqs  (2.14)  and  (2.15),  along 
the  arc-length.  The  calculations  of  the  curvature,  and  the  central  difference  scheme  used  for 
discretization  are  described  in  Appendices.  Making  use  of  the  symmetry  of  the  problem,  we 
solve  for  only  half  of  the  droplet  (x  >  0).  On  the  axis  of  symmetry  we  provide  boundary 
conditions  as  xp  =  0  and  u>  =  0.  Then  for  a  given  initial  shape  we  solve  the  following  matrix 
system  of  4 N  unknowns  and  AN  equations  to  obtain  the  flow  variables. 


1 


+  -  <  (h  —  a  cos  (3)  I  log  b - )+(acos/3) 


A  B  C  D  \ 

(  T  ^ 

f  o\ 

0  0  A  B 

dAt/dn 

0 

R  S  T  0 

UJ 

0 

£ 

0 

> 

y  dco/dn  J 

w 

A,  B,  C  and  D  are  dense  matrices  and  R,  U,  V  are  tri-diagonal  and  S,  T  and  W  are  diagonal 
matrices  at  the  free  surface.  At  the  solid  wall,  BtJ  =  Vij  =  Sij  and  S,.j  =  Tij  =  Uij  =  W,:I  = 
0.  The  matrix  system  is  then  solved  with  LU  decomposition  to  obtain  the  streamfunction 
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and  its  normal  derivative  at  each  boundary  node,  which  is  used  to  compute  the  velocities. 
Instantaneous  velocities  obtained  are  then  used  for  generation  of  a  new  droplet  profile  by 
marching  in  time  with  the  help  of  the  kinematic  condition, 

x(t  +  At)  =  x(t)  +  ux  At ,  and  y(t  +  At)  =  y(t)  +  uy  At,  (2.30) 

where,  ux  and  uy  are  the  components  of  velocity  u  in  x  and  y  direction.  To  get  a  position 
of  the  new  contact  line,  the  calculated  new  interface  is  made  to  intersect  the  solid  surface  at 
a  position  where  dynamic  contact  angle-velocity  correlation  (Eq.  2.6)  is  satisfied. 

During  the  time  evolution  of  the  droplet,  boundary  nodes  may  get  convected  along  the 
surface  and  thus  distribution  of  the  nodes  is  no  longer  uniform.  To  avoid  crowding/ piling  of 
nodes  in  some  region  as  opposed  to  other,  a  cubic  spline  is  fitted  after  every  time  iteration  to 
the  advanced  droplet  profile  and  boundary  nodes  are  again  redistributed.  Also,  often  a  high 
curvature  region  is  obtained  near  the  contact  line  during  the  spreading  process,  especially 
during  the  droplet  spreading  over  a  topography  due  to  very  large  contact  angles.  In  such 
cases,  the  air /fluid  interface  is  smoothed  out  by  fitting  a  cubic  polynomial  for  the  first  few 
nodes  from  the  contact  line. 

To  summarize,  a  FORTRAN  code  is  developed,  which  uses  only  an  initial  shape  of  a 
droplet.  From  this  initial  shape,  using  Boundary  Integral  Method,  it  computes  matrix  ele¬ 
ments  of  Eq.  2.29.  Matrix  system  is  solved  with  LU  decomposition  to  obtain  flow  variables. 
Using  flow  variables  boundary  nodes  are  advanced  in  time  to  get  the  new  droplet  profile. 
Contact  point  is  moved  with  a  slip  velocity,  dependent  upon  the  dynamic  contact  angle. 
Boundary  nodes  are  re-meshed  and  the  process  is  repeated  to  obtain  the  time  evolution. 

2.4.3  Code  Validation 

The  developed  fortran  code  is  validated  by  comparing  calculated  instantaneous  flow  variables 
with  those  obtained  from  an  exact  analytical  solution  of  a  simpler  problem.  Eq.  (2.9),  i.e. 
V4-0  =  0,  has  a  separable  solution  in  polar  co-ordinates 

ip(r,0)  =  rx+1fx(9),  (2.31) 

where,  A  is  any  constant  (real  or  imaginary)  and  the  general  form  of  f\(0)  is  given  by, 

f\{9)  =  Ai  cos  [(A  +  1)0]  +  A2  sin  [(A  +  1)6*]  +  A3  cos  [(A  —  1)0]  +  A4  sin  [(A  —  1)6] ,  (2.32) 
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Figure  2.8:  (a)  Schematic  of  flow  in  a  corner  between  a  free  and  a  no-slip  surface,  (b) 
Numerical  grid  employed  to  solve  this  problem  using  the  Boundary  Integral  Method 


Ai,  A‘2 ,  A3,  and  A4  are  constants  to  be  evaluated  from  boundary  conditions.  We  consider 
a  problem  of  flow  in  a  corner,  formed  by  a  no-slip  surface  and  a  free  surface  (Fig.  2.8),  for 
which  an  analytical  solution  can  be  obtained  Normal  and  tangential  stresses  are  applied  on 
the  free  boundary  and  a  no- slip  condition  at  the  solid  surface,  which  are  written  as, 


No  —  slip  : 

Tangential  Stress  Balance  : 
Normal  Stress  Balancce  : 


/( 0)  =  /'( 0)  =  0, 
/"(«)-/(«)( A2-1)  =  0, 
/'"(a)  +  (3A2  +  !)/'(«)  =  0. 


(2.33) 


A(«)  can  be  obtained  from  the  normal  stress  balance  condition.  We  use  a  =  1,  for  which 
A  =  0.814.  Now  we  have  three  boundary  conditions  and  four  unknown  constants.  Thus,  we 
arbitrarily  fix  A\  and  evaluate  the  other  constants  as, 


Ao  — 


A4  (A  +  1)  cos  [(A  +  l)a]  —  (A  —  1)  cos  [(A  —  l)a] 


A  +  1  sin  [(A  —  l)a] 

A3  =  —Ai,  and  A4  = 


sin  [(A  +  l)a] 

A 


(2.34) 


For  solving  this  problem  numerically  using  the  boundary  integral  method,  we  have  to  have 
a  closed  calculation  domain  which  is  achieved  by  joining  the  two  surfaces  smoothly  with  a 
circular  arc.  As  boundary  conditions  on  the  curved  surface,  we  use  the  actual  ip  and 
from  the  known  theoretical  solution.  The  matrix  system  in  Eq.  2.29  is  then  solved  with 
appropriate  boundary  conditions  to  obtain  ip,  u ;  and  their  normal  derivatives.  Fig.  2.9 
compares  the  flow  variables  obtained  from  the  code  with  analytical  values,  and  a  very  good 
agreement  between  the  two  is  observed.  The  singularity  at  the  apex  (in  between  26th  and 
27th  node)  has  been  handled  well  by  the  numerical  method. 
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Figure  2.9:  Comparison  of  numerical  results  (symbols)  with  analytical  solution  (solid  line) 
for  flow  in  a  corner  problem. 

2.5  Numerical  Results  and  Discussion 

The  first  step  in  understanding  the  dynamics  of  the  spreading  of  a  droplet  over  a  surface  with 
topology  is  to  study  the  equilibrium  shapes  of  the  drops  after  they  have  come  to  rest.  These 
equilibrium  shapes  represent  the  final  resting  states  of  the  trajectory  of  drops  advancing 
over  a  surface.  As  such,  in  what  will  become  clear  further  on,  these  shapes  are  important 
in  constructing  a  correlation  between  the  initial  shape  and  location  of  a  drop  placed  on  a 
surface  with  topology,  and  its  final  position  on  that  surface.  The  static  equilibrium  shapes 
are  discussed  in  Section  2.5.1.  Once  we  know  all  the  possible  equilibrium  shapes  for  a  droplet 
of  a  given  volume,  from  integration  of  Young-Laplace  equation,  we  consider  the  dynamics 
of  the  spreading  motion  over  surfaces  and  we  expect  that  the  equilibrium  shape  obtained  at 
the  end  of  the  dynamic  process  should  correspond  to  one  of  the  equilibrium  shape  predicted 
by  the  above  procedure.  Numerical  results  for  the  spreading  dynamics  are  obtained  from 
the  Boundary  Integral  Method.  In  all  simulations  presented  in  further  sections,  we  use  the 
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representative  case  of  a  2  mm  in  radius  (circular)  water  droplet.  The  initial  shape  of  the 
droplet  is  assumed  to  be  section  of  a  circle  with  radius  a,  and  the  droplet  is  symmetric  about 
Y-axis.  Time  evolution  is  obtained  by  solving  the  Stokes  equations  by  using  the  Boundary 
Integral  Method.  Equations  are  non-dimensionalized  with  initial  drop  radius,  a,  and  the 
characteristic  time  and  velocity  are  given  by; 

_/  _  /'  T//  _  pga2 

T  5  y 

pga  \i 

Thus,  by  writing  the  governing  equations  and  boundary  conditions  in  dimensionless  form, 
the  only  parameter  which  remains  is  the  Bond  number,  Bo  =  pga2 /q,  which  is  a  ratio  of  the 
body  forces  to  the  surface  tension  forces.  For  a  water  droplet  of  radius  a  =  2 mm,  density 
p  =  1000 kg/m3  and  with  gravitational  acceleration  g  =  9.81  m/s2,  the  Bond  number  is 
Bo  =  0.545  and  this  value  is  used  in  the  simulations.  First,  we  will  describe  the  numerical 
results  on  a  planar  hydrophobic  surface  (without  topography,  Section  2.5.2),  to  provide 
base  cases  to  compare  with  spreading  over  a  surface  which  has  topology.  We  then  describe 
the  spreading  dynamics  of  a  droplet  over  a  superhydrophobic  surface  consisting  of  regular 
patterns  of  rectangular  elevations  separated  by  rectangular  valleys  (Section  2.5.3). 

2.5.1  Static  equilibrium  shapes 

Consider  a  2-dimensional  droplet  of  liquid  with  surface  tension  7  and  density  p  resting  on  a 
surface  in  a  gravitational  field  of  g.  The  forces  acting  on  the  droplet  are  the  gravitational  or 
body  force,  and  an  interfacial  or  capillary  force  due  to  the  surface  tension  of  the  liquid.  On 
a  planar  surface,  which  is  perpendicular  to  gravity,  and  of  infinite  extent,  the  equilibrium 
shape  of  the  droplet  on  the  surface  is  uniquely  defined  if  two  parameters  are  specified:  volume 
of  the  droplet  and  a  contact  angle  it  makes  with  a  surface.  These  are  the  familiar  ’’sessile” 
forms.  Thus  for  a  given  volume  and  contact  angle,  the  droplet  interface  has  a  unique  shape 
as  described  by  the  Young-Laplace  equation, 


(2.35) 


AP  =  7/C  =  ^  (2.36) 

IX 

where  AP  is  the  pressure  difference  across  the  interface,  which  is  balanced  by  the  surface 
tension  times  the  curvature  n  of  the  interface,  where  R  is  the  radius  of  curvature.  This 
Young-Laplace  equation  can  be  integrated  to  obtain  the  interfacial  shape  at  the  equilibrium. 
Details  of  the  procedure  are  provided  in  the  Appendix  C. 

For  a  planar  surface,  for  a  given  volume  V,  by  changing  the  contact  angle  we  can  obtain 
different  equilibrium  shapes,  which  intersect  the  solid  surface  at  different  X  positions.  From 
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this  exercise,  we  can  obtain  the  extent  of  spreading  on  a  flat  surface,  presumed  to  be  infinite 
in  extent,  as  shown  in  Fig.  2.10  for  two  different  volumes  (VI  and  V2,  V2>V1),  where 
equilibrium  contact  angles  are  plotted  against  equilibrium  X  positions  for  contact  angles  (as 
measured  through  the  liquid  phase)  up  to  180°.  Figure  2.10  shows  that  on  a  flat  surface,  as 
expected,  a  droplet  which  makes  a  larger  contact  angle  on  a  flat  surface  spreads  to  a  lesser 
extent  as  compared  to  a  droplet  which  makes  a  smaller  contact  angle.  The  reason  for  this 
is  the  larger  static  contact  angle  generates  very  large  curvatures  and  capillary  pressures  at 
the  air/water  interface  near  the  contact  line,  which  can  support  the  weight  of  the  droplet 
and  restrict  its  spreading  extent.  Also,  as  shown  in  the  figure,  for  the  same  contact  angle 
a  droplet  with  a  larger  volume  spreads  to  greater  extent  as  the  curvature  forces  are  not 
sufficient  to  support  the  larger  weight  of  the  droplet.  Fig.  2.10  also  presents  the  equilibrium 
shapes  obtained  by  integrating  the  Young-Laplace  equation  for  different  equilibrium  contact 
angles,  namely  100°,  150°  for  the  volume  VI.  It  is  clear  from  the  interfacial  loci  depicted  in 
the  figure  that  the  larger  the  contact  angle  the  smaller  the  radius  of  curvature  at  the  apex 
of  the  sessile  drop. 

In  the  context  of  the  study  of  the  droplet  spreading  dynamics  on  surfaces  with  topography, 
the  final  resting  shapes  of  drops  placed  on  the  surface  contact  the  surface  on  sharp  90°  edges 
or  corners  where  they  become  pinned.  These  final  states  represent  sessile  drops  on  finite 
planar  surfaces  in  which  the  droplet  hangs  freely  over  the  edge  of  the  surface.  They  can 
be  computed  by  integrating  the  Young-Laplace  equation  through  the  point  at  which  the 
tangent  to  the  surface  becomes  parallel  to  gravity  and  the  shape  becomes  multi-valued  in 
the  coordinate  x.  For  a  fixed  volume  of  droplet  V,  as  the  equlibrium  X  (edge  or  corner) 
position  of  these  hanging  shapes  decreases  from  the  X  position  of  the  droplet  (with  volume 
V)  with  a  contact  angle  of  180°,  the  hanging  contact  angle  increases  above  180°  as  shown  in 
Figure  2.10.  The  shape  for  a  hanging  contact  angle  equal  to  192.5°  is  shown,  and  the  radius 
of  curvature  at  the  apex  is  larger  than  that  of  the  droplet  which  intersects  the  surface  at 
180°. 

2.5.2  Spreading  dynamics  over  a  planar  surface 

Here  we  present  numerical  results  for  a  droplet  spreading  on  a  planar  surface,  i.e.  without 
topography.  As  mentioned  before,  a  singularity  is  associated  with  the  three-phase  contact 
line  sliding  over  the  surface.  The  first  case  we  consider  here  is  when  the  prescribed  contact 
angle  that  the  liquid  makes  with  the  solid  is  equal  to  7r  or  180°.  In  this  case  the  droplet 
rolls  over  the  surface  (tank-treading  motion)  as  opposed  to  sliding  and  thus  the  contact 
line  singularity  is  absent.  Fig.  2.11(a)  shows  the  time  evolution  of  a  droplet.  Initially,  the 
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Figure  2.10:  Graph  on  the  left  shows  extent  of  spreading  on  a  flat  surface  obtained  from 
integration  of  the  Young-Laplace  equation  for  two  different  volumes  V2>V1,  Vl=12.46  mm2 
and  V2=28.11  mm2.  Plots  on  the  right  shows  different  equilibrium  shapes  obtained  for  same 
volume  VI  at  different  equilibrium  contact  angles  on  a  flat  surface  (shaded  in  grey).  These 
three  equilibrium  points  are  also  marked  on  the  extent  of  spreading  curve  for  volume  VI. 
Equilibrium  angles  larger  than  180°  are  obtained  if  droplet  is  pinned  indefinitely  at  a  corner 
point.  Droplets  are  symmetrical  about  Y  axis. 


shape  of  a  droplet  is  prescribed  to  be  a  perfect  circle  (of  radius  a  =  2mm)  just  touching 
the  flat  solid  surface  only  at  a  point.  For  a  water  droplet  of  radius  a  =  2 mm,  density 
p  =  1000 kg/m3  and  with  gravitational  acceleration  g  =  9.81  m/s2,  the  Bond  number  is 
Bo  =  0.545  and  this  value  is  used  in  the  simulations.  The  initial  shape  of  the  droplet  is 
approximated  by  a  polygon  of  N  =  280  sides.  The  matrix  system  generated  by  boundary 
integral  code  (1120x1120)  is  solved  in  a  loop  to  obtain  the  time  evolution  of  the  droplet. 
The  contact  line  is  advanced  by  making  the  air/ water  interface  to  intersect  the  solid  surface 
at  the  constant  contact  angle  of  180°,  which  achieved  by  fitting  a  parabola  through  the  last 
three  points.  As  a  function  of  time,  under  the  action  of  gravity  the  droplet  spreads  on  the 
flat  surface  with  tank-treading  motion  until  the  equilibrium  is  reached,  where  gravitational 
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forces  are  balanced  by  the  curvature  and  surface  tension  as  shown  in  the  Fig. 2. 11(a).  At  this 


(a)  Droplet  shape  evolution  with  time,  solid  (b)  Droplet  volume  change  during  the  spreading  dy- 
lines  represents  numerical  simulations  at  di-  namic  simulations, 
mensionless  time,  r  =  0,  0.05,  0.15,  0.3, 

0.5,  0.75,  1.05,  1.4,  1.8,  and  symbols  repre¬ 
sents  equilibrium  shape  obtained  from  inte¬ 
gration  of  Young-Laplace  equation.  Excel¬ 
lent  agreement  between  the  two  equilibrium 
shape  is  observed. 

Figure  2.11:  Droplet  spreading  over  a  flat  surface  with  constant  contact  angle  of  180°. 
Droplet  rolls  as  opposed  to  sliding. 

equilibrium  shape,  Young-Laplace  equation  is  satisfied  at  every  point  on  the  boundary,  i.e. 
pressure  drop  across  the  interface  is  equal  to  the  surface  tension  times  the  curvature.  This 
predicted  Young-Laplace  equilibrium  shape  is  also  plotted  in  the  Fig. 2. 11(a),  and  excellent 
agreement  between  the  two  equilibrium  shapes  is  obtained.  This  agreement  validates  our 
time  dependent  code.  As  a  check  on  the  incompressibility,  the  volume  change  during  the 
spreading  process  is  also  plotted  in  Fig. 2. 11(b)  and  a  negligible  volume  change  is  observed 
(about  0.01%). 

However,  as  mentioned  before,  on  a  flat  surface  a  water  droplet  cannot  achieve  a  contact 
angle  larger  than  120°.  Thus,  this  tank-treading  motion  is  not  possible  for  such  a  low  contact 
angles  and  the  droplet  rather  slides  or  slips  on  the  flat  surface.  Fig.2.12  shows  spreading 
dynamics  of  a  droplet  with  specified  static  advancing  contact  angle  of  120°.  The  drop  initial 
shape  prescribed  to  be  a  section  of  circle  with  initial  contact  angle  of  160°.  To  provide  slip  at 
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Figure  2.12:  Droplet  spreading  over  a  flat  surface  (shaded  in  gray)  with  static  advancing 
contact  angle  of  120°.  Numerical  results  are  shown  with  solid  line  at  dimensionless  time  r  = 
0,  0.21,  0.67,  1.49,  2.8,  4.8,  7.73,  11.87,  1757,  and  30.  Symbols  represents  equilibrium  shape 
obtained  from  integration  of  Young-Laplace  equation.  Excellent  agreement  between  the  two 
equilibrium  shape  is  observed. 


the  contact  line  Hoffman-Jiang  correlation  (Eq.  2.6)  is  used,  which  is  solved  iteratively  during 
each  time-step.  Again  the  final  equilibrium  shape  obtained  by  the  simulations  matches  very 
well  with  the  predicted  Young-Laplace  shape;  in  particular  the  final  position  of  the  contact 
line  is  predicted  from  Figure  2.10. 

2.5.3  Spreading  dynamics  over  topography 

Modeling  the  Hydrodynamic  Movement  of  a  Contact  Line  Over  a  Surface  To¬ 
pography  of  Pillars 

For  studying  the  spreading  dynamics  of  a  droplet  over  a  topography,  we  need  to  consider  the 
interaction  of  the  fluid  flow  with  the  surface  microstructures  from  a  hydrodynamic  point  of 
view.  The  droplet  interface  interacts  with  the  surface  through  the  contact  angle  hysteresis, 
and  velocity  dependent  dynamic  contact  angles.  Fig.  2.13  demonstrates  pictorially  these 
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interactions  with  different  surfaces.  Consider  first  a  fiat  surface,  Fig.  2.13(a).  It  is  well 
known  that  in  order  for  the  liquid  to  spread  on  the  flat  surface,  the  liquid-air  interface  has  to 
exceed  a  minimum  critical  contact  angle,  which  is  the  static  advancing  contact  angle  (9adv) 
measured  through  the  liquid  till  the  solid  surface.  This  9adv  is  an  intrinsic  property  of  the 
liquid-solid  system  in  consideration  and  is  independent  of  the  orientation  of  the  flat  surface. 
That  is,  if  the  liquid  is  spreading  downwards  on  a  vertical  surface,  Fig.  2.13(b),  the  minimum 
critical  contact  angle  required  for  the  on-set  of  the  downward  spreading  motion  is  still  0af[v 
but  measured  with  respect  to  the  vertical  wall.  To  extend  this  concept  of  the  static  advancing 
contact  angle  to  the  surface  with  a  topography  in  the  form  of  pillars,  consider  a  case  in  Fig. 
2.13(c)  where  liquid  spreads  on  the  top  horizontal  surface  of  a  pillar  with  minimum  critical 
advancing  angle  of  9crit}H{=  9adv)-  But  when  droplet  reaches  the  sharp  corner  point  of  this 
pillar,  the  contact  line  can  slip  vertically  downwards  only  if  it  exceeds  9crity  i=9adv  +  90°) 
measured  with  respect  to  the  horizontal  surface.  Thus,  for  all  dynamic  contact  angles  in  the 
range  9cruy  <  Odyn  <  9cruy  the  contact  line  does  not  move,  when  at  the  corner,  even  though 
interface  makes  contact  angle  larger  than  9adv  with  the  horizontal  plane.  In  other  words, 
these  sharp  corner  points  act  as  pinning  points  -  they  pin  the  contact  line  at  the  corner  - 
and  offers  multiple  static  contact  angles.  The  rest  of  the  droplet,  however,  is  free  to  move, 
and  if  the  droplet  can  come  to  an  equilibrium  shape  with  a  contact  angle  within  the  above 
mentioned  limits  then  further  droplet  spreading  motion  is  arrested.  On  the  other  hand,  if 
such  an  equilibrium  shape  does  not  exist  for  a  given  volume  of  the  drop  at  a  given  pinning 
position,  then  dynamic  contact  angle  exceeds  9cray  and  the  contact  line  begins  its  motion 
again.  In  this  process,  the  free  interface  of  the  droplet  bends  so  much  that  it  is  possible  for 
the  interface  to  touch  the  neighboring  post/pillar  and  then  that  touching  point  becomes  the 
new  contact  line.  The  droplet  then  spreads  further  on  the  top  horizontal  surface  of  the  next 
pillar. 


Figure  2.13:  Minimum  critical  contact  angles  for  the  on-set  of  droplet  contact  line  motion: 
(a)  forward  on  a  horizontal  surface,  (b)  downward  on  a  vertical  surface,  and  (c)  horizontal 
till  the  corner  and  downward  after  pinning  at  the  corner  for  contact  angles  in  the  range 

@crit,H  —  9dyn  —  @crit,V 
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The  Indefinite  Pinning  of  an  Advancing  Droplet  At  a  Corner 


Figure  2.14:  Droplet  initially  spreads  over  a  flat  surface  (shaded  in  gray)  with  static  advanc¬ 
ing  contact  angle  of  120°  and  suddenly  gets  pinned  at  a  corner  point  located  at  X  =  0.4. 
Numerical  results  are  shown  at  dimensionless  time  r  =  0,  0.1,  0.3,  0.6,  1,  1.5,  2.1,  2.8,  10.25 
and  28.25.  Apparent  equilibrium  contact  angle  observed  was  #*?  =  189.98°,  which  can  be 
located  on  extent  of  spreading  curve  (Fig.  2.10)  for  volume  VI  at  pinning  position  of  X  =  0.4 


As  an  illustration  of  the  pinning  of  the  droplet  contact  line  at  a  corner,  we  choose  the 
same  initial  shape  for  the  droplet,  a  section  of  a  circle  with  a  contact  angle  of  160  °  and  a 
volume  V=V1  (from  Figure  2.10)  =12.46  mm2,  and  an  advancing  contact  angle  (120°)  as  the 
previous  simulation  on  a  flat  surface  (Figure  2.12).  However,  an  edge  is  located  at  a  distance 
of  X  =  0.4  from  the  centerline  of  the  drop.  We  choose  this  value  for  the  edge  location 
because  it  is  clear  from  Figure  2.10  that  for  this  volume  a  hanging  drop  pinned  at  X  =  0.4 
can  exist  at  equilibrium  with  a  hanging  angle  of  approximately  190°.  Importantly  this  angle 
is  smaller  than  the  angle  ( 9crit,v  =  9adv  +  90  =  210°)  required  for  the  drop  to  slide  down  the 
vertical  surface.  Hence  we  anticipate  the  drop  to  move  first,  and  subsequently  get  pinned  at 
the  corner  until  equilibrium.  Incorporating  the  interaction  between  the  topography  and  the 
fluid  flow  over  a  textured  surface  as  described  above,  Fig.  2.14  shows  the  numerical  results 
for  a  droplet  spreading  on  a  flat  surface  up  to  the  distance  of  X  =  0.4  where  it  encounters 
the  corner  or  pinning  point.  As  the  simulation  begins,  the  droplet  initially  spreads  on  the 
flat  surface  as  the  dynamic  contact  angle  is  larger  than  9adv-  At  time  r  =  0.1455  the  contact 
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line  reaches  the  corner  point  (at  X  =  0.4)  and  the  contact  line  gets  pinned  indefinitely  at  this 
point.  But  as  the  rest  of  the  droplet  shape  is  still  not  an  equilibrium  shape,  the  interfacial 
nodes  continue  to  move.  The  interface  bends  over  while  remaining  pinned  at  the  corner  point 
until  finally  an  equilibrium  configuration  is  reached  where  all  the  forces  are  balanced.  The 
final  equilibrium  shape  is  realized  at  the  corner  point  with  the  equilibrium  contact  angle  of 
about  189.98°  in  agreement  with  Fig.  2.14.  This  very  large  contact  angle  is  possible  because 
it  still  lies  in  the  range  of  0crit,H  <  @eq  A  Verity  ■  The  contact  line  position  and  the  contact 
angle  as  a  function  of  time  indicate  that  the  contact  line  slips  untill  it  reaches  the  corner 
point,  where  it  sticks  until  equilibrium.  Similarly,  the  contact  angle  increases  monotonically 
when  the  contact  line  is  pinned  at  the  corner  and  reaches  an  equilibrium  value  of  189.98°. 

The  Spreading  Dynamics  of  Droplets  Over  Consecutive  Pillars 

With  reference  to  the  previous  simulation  of  movement  over  an  isolated  pillar  and  pinning 
at  the  downstream  edge,  advancement  over  a  microtextured  surface  of  pillars  involves  the 
fundamental  steps  of  movement  over  the  horizontal  surface  of  the  pillar  to  the  corner  edge 
to  form  a  hanging  drop,  followed  by  attachment  of  part  of  the  meniscus  to  the  next  pillar. 
For  this  attachment  to  occur,  the  next  pillar  must  be  close  enough  (that  is  the  gap  small 
enough)  so  that  the  hanging  drop  -  as  it  evolves  over  the  edge  -  contacts  the  next  pillar 
with  its  hanging  interface  before  coming  to  an  equilibrium  or,  if  it  cannot  reach  equilibrium, 
develop  a  contact  angle  with  respect  to  the  vertical  larger  than  the  advancing  angle  so  it 
moves  down  the  vertical  wall  of  the  pillar.  (In  the  latter  case,  attachment  can  still  occur 
as  the  movement  down  the  vertical  wall  allows  the  meniscus  to  come  in  contact  with  the 
next  pillar.)  Fig.  2.15(a)  shows  droplet  spreading  motion  over  a  surface  composed  of  a 
topography  in  the  form  of  a  width  of  the  pillar  =  0.05  and  gap  =  0.1.  Initial  shape  is  similar 
as  before  but  the  prescribed  advancing  contact  angle  is  100°.  Droplet  spreads  initially  on 
the  flat  surface  till  the  corner  at  X  =  0.4  and  contact  line  remains  pinned  while  the  interface 
bends  over.  And  before  it  can  come  to  an  equilibrium,  contact  line  touches  the  neighboring 
pillar  at  X  =  0.5  and  begins  its  spreading  motion  till  it  is  pinned  at  the  next  pinning  corner 
at  X  =  0.55  and  finally  comes  to  an  equilibrium  as  all  the  external  forces  and  the  contact 
angle  boundary  conditions  are  satisfied.  We  see  distinct  slip-stick-jump  motion  of  the  contact 
line  over  the  topography  while  trapping  an  air  underneath  the  droplet.  Slipping  over  the 
horizontal  surface,  sticking  to  the  corner  point  and  jumping  from  one  pillar  to  the  next.  Fig. 
2.15(b)  shows  this  motion  when  we  plot  the  contact  line  X  position  as  a  function  of  a  time 
and  contact  angle  also  varies  accordingly.  Thus,  just  by  prescribing  an  intrinsic  angle  of  only 
100°  and  having  a  topography  on  the  surface,  we  see  that  the  final  equilibrium  angle  reached 
by  the  droplet  is  very  large  (about  174.3°).  This  large  equilibrium  angle  is  possible  because 
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the  multiple  contact  angles  offered  by  the  sharp  pinning  point.  Final  equilibrium  position 
of  the  contact  line  and  the  contact  angle  can  be  located  on  the  extent  of  spreading  curve 
generated  before. 

Note  that  in  these  simulations,  after  droplet  contact  line  jumps  to  the  next  pillar,  liquid- 
air  meniscus  in  between  the  two  pillar  is  considered  to  be  flat.  That  is,  although  zero 
tangential  stress  condition  is  applied  on  the  free  air-water  interface  riding  over  trapped  air 
pocket,  normal-stress  balance  is  not  solved  for.  Instead,  ip  =  0  is  specified  on  the  flat 
interface.  This  simplification  is  valid,  as  long  as  the  radius  of  curvature  of  the  free  interface 
is  much  larger  than  the  distance  between  the  two  pillars,  which  is  usually  the  case.  We  have 
also  investigated  the  spreading  dynamics  while  considering  the  curvature  of  the  interface, 
and  found  that  the  curvature  is  indeed  small  and  the  final  equilibrium  shape  and  overall 
dynamics  of  the  droplet  is  not  affected.  Spreading  dynamics  of  the  droplet  is  dominated 
by  the  contact  line  motion  and  free  interface  above  the  solid  but  it  is  unaffected  by  solid 
topography  or  shape  of  the  liquid  menisci  in  the  interior  of  the  droplet.  Some  of  the  numerical 
results  capturing  dynamics  of  the  liquid  meniscus  can  be  found  in  Appendix  D  at  the  end  of 
the  chapter. 

For  the  topography  considered  in  our  simulation,  if  we  compute  the  average  solid  area 
fraction  (<fis)  available  for  a  droplet  to  spread  on  a  composite  surface  made  of  pillar  width 
=  0.05  and  gap  =  0.1,  we  get  os  =  0.333  and  corresponding  predicted  contact  angle  by 
Cassie-Baxter  relation  (Eq.  ??)  is  136.43°.  There  is  clearly  a  mismatch  between  the  two 
equilibrium  contact  angle  on  a  similar  surface  topography.  But  the  equilibrium  shape  we 
obtained  by  our  simulation  is  a  perfectly  valid  shape  for  a  given  volume  and  it  falls  exactly 
on  the  straight  line  predicted  by  integrating  Young-Laplace  equation  (Fig.  2.10).  That 
is,  in  the  equilibrium  shape  obtained  in  our  simulation  the  pressure  difference  across  the 
interface  is  exactly  equal  to  the  curvature  times  the  surface  tension  at  every  single  point 
on  the  interface.  Reason  behind  this  discrepancy  is  that  the  relations  based  on  variational 
energy  principles  consider  the  surface  is  composed  of  micro-structures  so  small  that  surface 
is  assumed  to  have  some  average  properties,  i.e.  effective  solid-vapor  surface  energy  (Ysv) 
and  effective  solid- liquid  surface  tension  (7 *SL),  as  mentioned  before.  And  droplet  is  assumed 
to  be  spreading  over  a  macroscopically  flat,  but  with  an  average  lower  energy  surface.  This 
variational  energy  approach  clearly  does  not  account  for  interaction  of  the  flow  with  the 
topography.  This  interaction,  as  described  earlier,  can  result  in  any  equilibrium  contact 
angle  if  0crit)H  <  0*q  <  0crity  inequality  is  satisfied.  And  final  equilibrium  contact  angle 
achieved  depends  on  the  history  of  the  droplet. 

In  simulation  presented  in  Fig.  2.15(a)  droplet  never  feels  the  presence  of  the  second 
pillar  situated  at  X  =  0.65.  However,  if  we  push  the  droplet  down  forcefully  -  on  the 
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(a)  Time  evolution  of  the  droplet  shape  near  the  contact  line  region.  Droplet  achieves 
a  final  equilibrium  shape  at  the  pinning  point  located  at  X  =  0.55. 


Time 


Slip 


(b)  Left  hand  side  plot  shows  slip-stick-jump  motion  of  the  contact  line,  contact  line  slips  on  the  horizontal 
surface,  sticks  after  reaching  the  corner  point  and  jumps  from  one  pillar  to  the  next.  Corresponding  dynamic 
contact  angles  are  plotted  on  the  right. 


Figure  2.15:  Droplet  spreading  dynamics  over  a  topography  composed  of  solid  width  =  0.05 
and  gap  =  0.1  for  advancing  contact  angle  of  100°. 


same  textured  surface  as  above  -  it  sees  or  explores  the  pillars  located  at  farther  distances. 
For  example,  in  Fig  2.16(a),  Fig  2.16(b),  Fig  2.16(c)  and  Fig.  2.16(d)  where  droplet  initial 
shape  itself  has  a  contact  line  already  located  on  more  distant  pillars.  In  such  cases  droplet 
still  spreads  ahead  as  the  contact  angle  it  makes  with  the  horizontal  surface  is  larger  than 
Oadv  —  100°  and  comes  to  an  equilibrium  at  the  next  available  pinning  point.  Thus,  depending 
upon  where  the  final  pinning  point  is  located,  the  apparent  equilibrium  contact  angle  varies. 
Figures  show  that  for  pinning  points  located  at  X  =  0.7,  X  =  0.85,  X  =  1.0  and  X  =  1.15, 
apparent  equilibrium  contact  angles  obtained  on  the  exactly  same  textured  surface,  for  the 
same  volume  droplet,  are  161.98°,  148.29°,  135.23°  and  122.57°  respectively.  And  all  these 
equilibrium  shapes  fall  on  the  extent  of  spreading  curve,  as  shown  in  Fig.  2.17.  Fig.  2.17 
also  plots  the  Cassie-Baxter  equation  predicted  angle  on  the  curve  and  we  see  that  contact 
angles  larger  or  smaller  than  Cassie-Baxter  prediction  are  possible,  provided  they  are  within 
the  range  0crit^H  <  0*eq  <  9crity  and  we  have  pinning  points  at  respective  locations.  In 
other  words,  for  a  textured  surface,  there  is  no  one  single  apparent  equilibrium  contact  angle 
(like  obtained  from  global  energy  minimization  or  Cassie-Baxter  relation)  that  describes  the 
surface  fully. 

In  above  simulations,  droplet  contact  line  never  starts  its  downward  motion  once  it  is 
pinned.  In  order  to  see  this  penetrating  motion,  contact  angle  has  to  exceed  the  critical 
angle  of  9crity.  To  demonstrate  this  behavior  we  consider  a  case  where  intrinsic  advancing 
contact  angle  is  90°.  Therefore,  when  contact  angle  exceeds  180°,  contact  line  will  resume 
its  downward  motion  and  equilibrium  angles  larger  than  180°  cannot  be  realized  as  in  the 
case  of  Fig.  2.14.  Fig. 2. 18  and  Fig  2.19  shows  slip-stick-penetration-jump  behavior  for  two 
different  pillar  densities,  where  contact  line  after  sticking  penetrates  down  the  pillar  before 
the  interface  touches  the  next  surface.  Extent  of  downward  penetration  is  not  significant  but 
can  be  seen  in  the  simulations.  And  again,  final  equilibrium  position  and  the  contact  angle 
can  be  located  on  the  extent  of  spreading  curve. 

2.6  Experimental  Studies 

This  section  presents  experimental  investigations  of  superhydrophobic  surfaces.  In  section 
2.6.1  we  demonstrate  techniques  used  for  fabricating  regular/periodic  patterns  on  the  surface, 
which  are  of  the  micron-scale,  and  characterizing  these  surfaces  for  their  micro-structures  as 
well  as  static  equilibrium  contact  angles.  In  section  2.6.3,  dynamics  of  a  droplet  spreading 
over  a  micro-textured  surface  is  presented  while  maintaining  conditions  somewhat  similar 
to  the  numerical  section  so  that  a  qualitative  comparison  could  be  made  between  the  two. 
Experimentally  droplet  dynamics  is  studied  by  two  procedures:  firstly  under  gravity  driven 
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(a)  Droplet  pins  at  X  =  0.7  9*  =  161.98°  (b)  Droplet  pins  at  X  =  0.85  9*  =  148.29° 


(c)  Droplet  pins  at  X  =  1.0  9*  =  135.23°  (d)  Droplet  pins  at  X  =  1.15  9*  =  122.57° 


Figure  2.16:  Droplet  spreading  on  a  topography  of  pillar  width  =  0.05  and  gap  =  0.1  with 
different  initial  shape,  shows  different  apparent  equilibrium  contact  angels. 
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Figure  2.17:  Extent  of  spreading  curve  for  volume  VI,  on  which  different  equilibrium  points 
are  plotted  for  a  given  surface  texture  obtained  from  the  dynamics  of  the  spreading  problem 
with  different  initial  shapes.  Shaded  area  shows  the  range  of  contact  angles  possible  (0critH  < 
6*eq  <  Ocrity)  for  the  prescribed  intrinsic  surface  energy  of  the  liquid-solid  system  0a(iv .  V  is 
the  Cassie-Baxter  relation  predicted  apparent  contact  angle. 


spreading  and  secondly  with  controlled  spreading,  i.e.  by  increasing  volume  of  the  droplet 
gradually  with  time. 

2.6.1  Microfabrication  and  Static  Contact  Angle  Measurements 

Fabrication  of  Patterned  Surfaces 

Soft-lithographic  technique  was  used  for  fabricating  superhydrophobic  surfaces  from  PDMS, 
i.e.  poly  (dimethylsiloxane).  Fig.  (2.20)  illustrates  the  major  fabrication  steps.  Microstruc¬ 
tures  were  imposed  onto  the  PDMS  by  casting  it  over  a  master.  For  preparation  of  the 
master  from  silicon  wafer  we  first  clean  the  wafer  with  standard  technique  of  series  of  chem¬ 
ical  washes,  acetone,  iso-propyl  alcohol  and  de-ionized  water  and  blow  dried  with  nitrogen 
gas.  We  spin-coated  a  thin  layer  of  negative-toned  photoresist  (SU-8  2050,  MicroChem  Cor- 


41 


(a)  Time  evolution  of  the  droplet  shape  near  the  contact  line  region.  Droplet  achieves 
a  final  equilibrium  shape  at  the  pinning  point  located  at  X  =  0.55. 


(b)  Left  hand  side  plot  shows  slip-stick-jump  motion  of  the  contact  line,  contact  line  slips  on  the  horizontal 
surface,  sticks  after  reaching  the  corner  point  and  jumps  from  one  pillar  to  the  next.  Corresponding  dynamic 
contact  angles  are  plotted  on  the  right. 


Figure  2.18:  Droplet  spreading  dynamics  over  a  topography  composed  of  solid  width  =  0.05 
and  gap  =  0.1  for  advancing  contact  angle  of  90°.  Liquid  starts  its  downward  motion  on  the 
first  corner  point  before  the  interface  touches  the  adjacent  pillar. 
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(a)  Time  evolution  of  the  droplet  shape  near  the  contact  line  region.  Droplet  achieves 
a  final  equilibrium  shape  at  the  pinning  point  located  at  X  =  0.6. 


(b)  Left  hand  side  plot  shows  slip-stick-jump  motion  of  the  contact  line,  contact  line  slips  on  the  horizontal 
surface,  sticks  after  reaching  the  corner  point  and  jumps  from  one  pillar  to  the  next.  Corresponding  dynamic 
contact  angles  are  plotted  on  the  right. 


Figure  2.19:  Droplet  spreading  dynamics  over  a  topography  composed  of  solid  width  =  0.05 
and  gap  =  0.05  for  advancing  contact  angle  of  90°.  Downward  motion  of  the  interface  on 
the  second  pillar  is  more  clearly  seen  before  the  interface  touches  the  adjacent  pillar. 
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poration,  Newton,  MA)  onto  a  cleaned  silicon  wafer.  The  desired  film  thickness  is  achieved 
by  tuning  the  spinning  conditions  such  as  rpm,  acceleration  and  spinning  time.  The  silicon 
wafer  along  with  the  SU-8  layer  is  then  baked  (pre-exposure  bake)  over  a  hot  plate  before 
exposing  to  UV  light  (at  365  nm)  through  a  transparency  mask.  High  resolution  trans¬ 
parency  mask  were  printed  from  Pageworks,  Cambridge,  MA,  which  consist  of  the  negatives 
of  required  designs  to  be  patterned  on  the  PDMS  surface.  Exposing  the  SU-8  layer,  on 
silicon  wafer  substrate,  to  UV  light  at  prescribed  exposure  energy  causes  the  photoinduced 
cross-linking  reaction  and  only  part  of  the  SU-8  which  is  exposed  through  the  transparency 
gets  cross-linked.  After  UV  exposure  the  photoresist  layer  is  baked  again  on  a  hot  plate 
(post-exposure  bake).  Patterns  were  developed  on  the  master  by  immersing  the  wafer  in 


Figure  2.20:  Fabrication  of  patterned  PDMS  surface:  (a)  UV-exposure  of  SU-8  resist  through 
transparency  mask,  (b)  Master  prepared  after  washing  off  unpolymerized  SU-8,  (c)  PDMS- 
prepolymer  mixture  poured  over  the  master  and  cured  in  oven,  (d)  PDMS  surface  with 
micro-patterns  is  obtained  after  peeling-off  from  the  master. 


the  developing  solution  (SU-8  developer,  MicroChem  Corporation,  Newton,  MA),  to  remove 
uncross-linked  SU-8  from  the  wafer.  Optimum  parameters  for  the  SU-8  fabrication  process 
(including  spin-coating  rpm,  baking  times,  UV  light  exposure  energy)  are  obtained  from  the 
datasheet  provided  by  MicroChem  corporation61.  To  make  the  PDMS  peel-off  easy,  the  ex¬ 
posed  silicon  wafer  on  the  master  was  made  hydrophobic  by  reacting  it  with  perfluorosilane. 
The  surface  modification  reaction  was  carried  out  by  dissolving  10/d  of  (Heptadecafluoro- 
l,l,2,2-tetrahydrodecyl)trichlorosilane  (Gelest  Inc,  Morrisville,  PA))  into  25ml  of  chloroform 
and  then  the  master  was  immersed  into  the  solution.  After  2  hours,  the  master  was  removed, 
blow  dried  with  N2  and  heat  treated  for  20  min  at  115°C. 

For  the  the  fabrication  of  the  PDMS  surface,  the  PDMS  prepolymer  mixture  is  prepared 
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by  mixing  the  prepolymer  with  a  curing  agent  (Sylgard  184  Silicone  Elastomer  Kit,  Dow 
Corning,  Midland,  MI)  in  a  10:1  weight  ratio.  The  prepolymer  mixture  was  then  thoroughly 
mixed  and  defoamed  in  a  centrifugal  Thinky  Mixer  (AR-100,  Thinky  Corporation,  Tokyo, 
Japan).  The  PDMS  prepolymer  mixture  was  then  poured  onto  the  master  and  trapped  air 
in  the  patterns  was  removed  by  applying  vacuum  for  30  min.  The  PDMS  mixture  was  then 
cured  in  an  oven  at  65°C  for  2  hours.  The  PDMS  replica  was  then  gently  peeled  off  from 
the  master. 

Different  structured  patterns  of  PDMS  surfaces  were  prepared  like  a  regular  array  of 
square  pillars,  square  holes,  and  circular  pillars.  SEM  images  of  the  patterned  PDMS  surfaces 
were  taken  (Fig.  2.21)  with  a  Zeiss  Evo  40  operated  at  lOkV.  Images  were  recorded  under 
variable  pressure  (VP)  mode  at  40  Pa.  Samples  were  also  characterized  under  an  optical 
microscope  to  get  quantitative  data  about  the  structures. 


Figure  2.21:  Scanning  Electron  Microscopy  images  of  the  fabricated  PDMS  surfaces:  (a) 
Square  pillars  (SP),  (b)  Square  wells  (SW),  (c)  and  (d)  Circular  pillars  (CP).  Scale  bar  is 
shown  at  the  bottom-right  corner  of  each  image. 


2.6.2  Contact  Angle  Measurements 

As  a  measure  of  the  hydrophobicity,  water  contact  angle  measurements  were  performed  on 
different  PDMS  surfaces.  Measurements  were  performed  with  deionized  water  and  surfaces 
were  cleaned  to  remove  dust  particles  simply  by  blowing  nitrogen  gas  over  it.  A  pendant 
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Figure  2.22:  Schematics  of  Video  Microscopy  apparatus  used  for  imaging  the  droplets:  A, 
Objective  lens;  B,  collimating  lens;  C,  pinhole;  D,  focussing  lens 


Table  2.1:  Comparing  experimentally  observed  contact  angle  with  that  predicted  from 
Wenzel  and  Cassie-Baxter  theory;  where  a  is  side  of  a  square,  d  is  diameter  of  a  circle,  p  is 
square  pitch,  and  height  h  =  100  pm. 


No. 

Surface  Topology  (in  pm) 

r 

0S  (in%) 

ec 

0W 

@exp 

1 

Flat 

i 

100 

- 

- 

105 

2 

SP,  a  =  40,  p  =  100 

2.6 

16 

151.8 

132.3 

144 

3 

SW,  a  —  170,  p  =  200 

2.85 

27.75 

142.6 

137.5 

145 

4 

CP,  d  =  110,1?  =  250 

1.55 

15 

152.7 

113.6 

141 

5 

CP,  d  =  130,  p  =  250 

1.65 

21 

147.6 

115.3 

128.4 

bubble/sessile  drop  video  microscopy  setup  was  used  to  image  the  droplet,  from  which  the 
contact  angle  was  then  obtained.  Fig.  (2.22)  shows  the  video  microscopy  setup,  in  which  light 
is  collimated  by  using  a  series  of  optics,  and  then  imaged  by  a  CCD  camera  which  is  connected 
to  a  computer.  Images  recorded  are  shown  in  the  Fig. (2.23)  with  respective  contact  angles 
for  a  constant  droplet  volume  of  6 pi.  Light  coming  from  beneath  the  drop  is  an  indication 
of  the  Cassie-Baxter  kind  of  wetting  and  the  droplet  resting  on  a  composite  surface  of  solid 
and  air.  Approximately  a  40°  jump  in  the  water  contact  angle  was  observed  between  the  flat 
and  the  roughened  surfaces.  Table  2.1  compares  the  experimentally  observed  water  contact 
angles  ( 6exp )  with  the  predicted  values  from  Wenzel’s  approach  (9W)  and  the  Cassie-Baxter’s 
approach  (9C).  We  can  see  that  the  experimental  values  are  closer  to  that  predicted  by  the 
Cassie-Baxter’s  relation.  This  also  conforms  that  the  droplet  is  riding  on  a  composite  surface 
rather  than  having  a  fully  wetted  contact  with  the  solid.  In  these  experiments  droplet  was 
placed  on  the  surfaces  with  a  hand-held  syringe  and  even  with  a  little  shaky  motion  of  hand 
it  is  possible  that  the  droplet  contact  line  is  forced  to  jump/explore  the  neighboring  pillars 
and  thus  contact  angle  is  much  lower  than  180 - °. 
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Figure  2.23:  Contact  angles  as  measured  from  the  images  obtained  from  video  microscopy 
on  (a)  Square  pillars  (SP);  (b)  Square  wells  (SW);  (c)  and  (d)  Circular  pillars  (CP). 

2.6.3  Experimental  Results  of  Droplet  Spreading  over  a  Superhy¬ 
drophobic  Surface 

Gravity  driven  spreading  of  a  droplet 

To  study  the  gravity  driven  spreading  dynamics  of  a  droplet  and  to  have  somewhat  similar 
conditions  like  our  simulations,  we  fabricate  2-D  surfaces  i.e.  surfaces  have  microstructures 
in  the  form  of  long  parallel  grooves  as  opposed  to  pillars.  SEM  images  of  such  fabricated 
anisotropic  surfaces  are  shown  in  Fig.  2.24.  Also,  this  spreading  dynamics  is  a  very  fast; 
for  water  if  we  compute  characteristic  velocity,  V  =  pga2 / p,  for  water  V  ~  0(10)  m/s 
and  corresponding  Reynolds  number,  Re  =  aV'p/p  ~  O(105).  Thus,  stokes  flow  assumption 
is  definitely  not  valid.  Therefore,  we  use  glycerol  instead  of  water  because  glycerol  exhibit 
very  similar  contact  angles  on  the  flat  PDMS  surface  as  that  of  water  (9giy  =  100°  and 
9wat  =  103°).  The  advantage  using  glycerol  over  water  is  that  its  viscosity  is  1.2  Pa.s  (1200 
times  viscosity  of  the  water).  Hence,  characteristic  velocity  for  glycerol,  V  ~  0(0.01)  and 
corresponding  Reynolds  number  is  Re  ~  0(0.01). 

To  capture  the  dynamics  of  the  droplet,  glycerol  drop  was  gently  touched  the  solid  surface 
while  still  attached  to  the  needle  and  needle  was  gently  pulled  up  so  that  the  drop  snaps  off 
from  the  needle  and  spreads  under  gravity  over  the  parallel  grooved  surface.  Its  dynamics 
is  visualized  along  the  transverse  direction  (perpendicular  to  the  grooves)  with  a  fast  video 
camera  capable  of  capturing  few  thousand  frames  per  second.  Fig.  2.25  shows  an  experiment 
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(a) 


Figure  2.24:  2-D  anisotropic  parallel  grooved/striped  surfaces,  (a)  Schematic  of  striped 
surface  with  width  w,  height  h  and  gap  between  the  two  stripes  g.  Length  of  these  stripes 
is  1  cm,  much  longer  as  compared  to  w,g  and  h.  (b)  and  (c)  Scanning  Electron  Microscopy 
images  with  width  w  of  175  /im  and  375  /im  respectively,  g  =  125/un  and  h  =  100/un  in 
both  cases. 


of  gravity  driven  dynamics  of  glycerol  droplet  over  a  grooved  surface  of  solid  width  175/iin 
and  gap  between  the  two  pillars  125/im.  Images  were  captured  by  the  fast  video  camera  at 
2000  frames/second.  Initially  at  time  t=0  s  droplet  is  attached  to  the  needle  and  touching 
the  top  of  three  solid  surfaces  and  contact  line  is  pinned  at  the  two  corners  (one  on  the  either 
sides).  During  the  snap-off  process,  as  the  needle  was  pulled  up,  interface  bends  over  while 
still  remaining  pinned  at  the  same  corner  points  and  while  doing  so  interface  just  touches  the 
neighboring  pillar  (on  the  left)  at  time  t=0.0175  s.  This  initial  contact  point  then  becomes 
the  new  contact  line  position  and  new  contact  line  slips  over  the  top  surface  of  the  pillar 
from  time  t=0.0175  s  to  t=0.0265  s,  at  which  time  reaches  the  next  pinning  point  located 
at  the  corner  of  the  pillar.  After  time  t=0. 0265s  contact  line  remains  pinned  (stuck)  to  the 
corner  point  and  the  contact  angle  increases  till  the  droplet  finally  achieves  an  equilibrium 
at  the  new  pinning  position.  During  this  process  neck  connecting  the  droplet  with  the  needle 
thins  and  finally  breaks.  Equilibrium  contact  angle  measured  at  the  end  of  this  process  is 
168.3°.  The  light  coming  from  underneath  the  droplet  again  indicates  the  Cassie-Baxter 
kind  of  wetting  where  droplet  rests  on  a  composite  surface  of  air  and  solid.  For  this  surface 
if  we  calculate  the  average  solid  area  fraction  in  the  composite  surface,  we  get  6S  =  0.5833. 
By  using  Cassie-Baxter  relation  Eq.2.35,  predicted  equilibrium  contact  angle  for  this  is 
9*  =  121.19°.  This  clearly  indicates  that  Cassie-Baxter ’s  relation  fails  to  explain  such  a  large 
equilibrium  contact  angle.  Because  droplet  while  spreading  does  not  see  an  average  solid 
area  fraction  of  the  composite  surface  but  local  surface  properties.  On  the  other  hand,  in  our 
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Figure  2.25:  Gravity  driven  spreading  motion  of  a  glycerol  droplet  over  a  striped  surface 
of  solid  width  175//m  and  gap  125//m.  Images  were  recorded  by  a  fast  video  camera  at 
2000  frames/sec.  Droplet  contact  line  shows  distinct  slip-stick  motion  after  jumping  to  the 
adjacent  pillar.  Droplet  interface  touches  the  pillar  on  the  left  and  slips  over  the  surface  till 
it  remains  stuck  at  the  corner  point. 


simulations  we  have  observed  multiple  equilibrium  contact  angles  depending  upon  the  initial 
base  radius  and  for  smaller  initial  base  radius  we  have  observed  contact  angles  of  about  170°. 

Similar  experiment  of  gravity  driven  spreading  is  shown  in  Fig.  2.26  on  a  grooved  surface 
of  solid  width  3 75 pm  and  gap  125pm.  This  experiment  more  clearly  shows  contact  line’s 
slip-stick  motion  after  contact  line  jumps  from  its  initial  pinning  point  to  neighboring  pillar 
on  the  left  at  time  t  =  0.0417  s.  Images  were  recorded  at  6000  frames/sec.  After  going 
through  similar  steps  droplet  again  comes  to  equilibrium  at  the  corner  point  with  final 
contact  angle  of  163  ±  2°.  Average  solid  area  fraction  of  the  composite  surface  is  </>s  =  0.75 
and  Cassie-Baxter  relation  predicts  9*  =  113.14°. 
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Figure  2.26:  Gravity  driven  spreading  motion  of  a  glycerol  droplet  over  a  striped  surface 
of  solid  width  375pm  and  gap  125pm.  Images  were  recorded  by  a  fast  video  camera  at 
6000  frames/sec.  Droplet  contact  line  shows  distinct  slip-stick  motion  after  jumping  to  the 
adjacent  pillar.  Droplet  interface  touches  the  pillar  on  the  left  and  slips  over  the  surface  till 
it  remains  stuck  at  the  corner  point. 


Controlled  spreading  of  a  droplet 

Slip-stick  motion  of  the  contact  line  is  demonstrated  in  section  2.6.3  where  dynamics  was 
too  fast  and  it  is  captured  with  a  fast  video  camera.  In  order  to  have  a  closer  look  at  the 
contact  line  pinning  and  available  multiple  contact  angles  as  well  as  to  minimize  the  effect  of 
the  second  curvature  in  3-D  droplet  a  controlled  spreading  experiments  were  performed.  In 
these  experiments  droplet  contact  line  does  not  advance  by  gravity  driven  spreading  motion 
but  due  to  increase  in  volume  of  the  droplet.  Experimental  arrangement  is  as  follows. 
Microtextured  surface  used  is  again  2-D,  parallel  grooved  surface.  Water  is  used  as  a  liquid 
in  these  experiments  as  by  controlling  the  rate  of  increase  in  the  volume  of  the  droplet 
we  can  still  have  stokes  flow  conditions.  Volume  of  a  droplet  is  increased  with  the  help 
of  a  syringe  pump  which  pumps  the  water  at  a  very  small  rate  through  a  needle  which  is 
attached  at  its  end  to  a  rectangular  slit  arranged  above  the  surface  with  slit  length  aligned 
along  the  grooves.  Water  coming  out  of  this  slit  forms  a  cylindrical  droplet  with  axis  of 
cylinder  parallel  to  the  grooves.  This  cylindrical  shape  of  the  droplet  minimized  effect  of 
3-D  curvature  and  resembles  more  closely  with  the  simulation  conditions.  Schematic  of  this 
experimental  arrangement  is  shown  in  Fig.  2.27. 

Fig.  2.28  shows  one  such  typical  experiment  where  volume  of  the  droplet  is  increased 
gradually  as  the  cylindrical  drop  forms  over  a  surface  consisting  of  stripes  of  275  pm  solid 
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Figure  2.27:  Schematic  of  the  experimental  arrangement  for  the  controlled  spreading  of  a 
cylindrical  droplet  over  a  parallel  grooved  surface.  Droplet  volume  in  increased  by  flow 
coming  in  through  the  slit. 


width  and  125  /im  gap.  Flow  rate  to  the  droplet  was  maintained  at  15  /tl / min  and  droplet 
spreading  is  observed  in  a  transverse  direction  through  a  CCD  camera.  Camera  is  focused 
in  the  middle  of  the  cylindrical  droplet  in  order  to  avoid  any  end  effects  and  as  a  result  the 
surface  topography  appears  blurred  or  washed  out.  (When  camera  is  focused  on  the  front 
end  of  the  surface  topography  looks  very  sharp,  as  in  the  previous  results  of  gravity  driven 
spreading.)  Images  were  captured  at  the  rate  of  1  frame/sec.  Figure  shows  time  dependent 
droplet  shape  in  the  vicinity  of  the  contact  line,  which  is  initially  pinned  at  a  corner  point. 
When  the  flow  into  the  droplet  is  started  contact  line  remains  pinned  at  the  same  location, 
however  rest  of  the  droplet  is  free  to  move  and  as  a  function  of  time  droplet  goes  through 
multiple  contact  angles  from  157°  till  very  close  to  180°  after  47  sec  where  the  free  interface 
just  touches  the  next  pillar.  As  soon  as  the  interface  touches  the  neighboring  pillar  on  the 
left,  contact  line  quickly  slips  and  sticks  to  the  downstream  end  of  the  corner  point  within  a 
second  where  contact  angle  is  again  about  157°  and  the  process  continues.  Fig.  2.29  shows  a 
very  similar  experiment  on  a  surface  with  topography  of  solid  stripes  width  175  /iin  and  gap 
between  the  stripes  125  /im.  In  this  experiment  flow  rate  through  the  slit  was  maintained  at 
10/d/miii  and  images  were  captured  at  1  frame/sec.  These  two  experiments  clearly  confirms 
the  presence  of  sharp  pinning  points  and  droplet  can  achieve  any  contact  angle  from  the 
range  permissible. 
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Figure  2.28:  Controlled  spreading  of  a  2-D  droplet  over  a  striped  surface  with  w  =  275/iin, 
g  =  125/rm  and  h  =  100/rm.  Flow  rate  through  the  slit  =  15  //I / min,  frame  capture  rate  = 
1  frame/sec.  Droplet  shows  pinning  at  the  corner  point  while  its  contact  angle  increases  till 
time  t=47  sec  at  which  time  the  contact  line  jumps  to  the  next  pillar. 
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Figure  2.29:  Controlled  spreading  of  a  2-D  droplet  over  a  striped  surface  with  w  =  175/iiri, 
g  =  125/rm  and  h  =  100/rm.  Flow  rate  through  the  slit  =  10  //I / min,  frame  capture  rate  = 
1  frame/sec.  Droplet  shows  pinning  at  the  corner  point  while  its  contact  angle  increases  till 
time  t=31  sec  at  which  time  the  contact  line  jumps  to  the  next  pillar. 


2.7  Conclusion 

From  our  numerical  simulations,  what  we  see  is  the  superhydrophobicity  is  attributed  to 
the  nature  of  the  pinning  points  and  the  final  equilibrium  contact  angle  is  dependent  upon 
where  the  pinning  points  are  located.  Closer  the  distance  of  the  final  pinning  point  to  the 
axis  of  symmetry,  for  a  given  volume,  larger  is  the  observed  equilibrium  contact  angle.  Once 
the  droplet  attains  equilibrium  at  a  corner,  it  doesn’t  know  the  presence  of  the  pinning 
points  which  are  ahead  of  the  contact  fine.  In  our  simulations,  we  also  observe  very  large 
equilibrium  contact  angles,  much  larger  than  those  predicted  by  Cassie-Baxter  equation  for 
similar  surface  features.  Major  difference  between  the  two  is  that  Cassie-Baxter  equation  is 
dependent  upon  effective  surface  properties  while  our  hydrodynamic  model  is  based  on  the 
interaction  of  fluid  droplet  contact  fine  with  local  topography  of  the  surface.  On  the  textured 
surface  there  are  number  of  equilibrium  shapes  of  the  droplet  are  possible  at  different  pinning 
locations  with  corresponding  contact  angles.  Final  equilibrium  shape  attained  by  the  droplet 
is  based  on  the  history  of  the  droplet,  i.e.  how  droplet  achieved  that  equilibrium.  When 
droplet  spreads  on  the  surface  with  a  very  small  initial  base  radius,  it  first  sees  pinning  points 
that  are  close  to  the  symmetry  axis  and  hence  final  equilibrium  contact  angle  obtained  is 
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larger.  If  for  some  reason  (often  in  experiments  due  to  inertia  and  relative  motion  of  the 
depositing  needle  with  the  surface)  the  droplet  initially  starts  spreading  with  a  larger  base 
radius,  then  as  it  spreads  it  will  come  across  the  pinning  points  that  are  ahead  of  the  initial 
base  radius.  In  that  case,  droplet  will  be  pinned  at  a  farther  distance  from  the  symmetry 
axis  and  for  a  given  volume,  it  will  achieve  lower  equilibrium  contact  angle  corresponding 
to  the  location  of  the  corner  point.  Thus  apparent  equilibrium  contact  angle  achieved  is 
dependent  on  location  of  pinning  points,  volume  of  the  droplet,  and  history  of  the  droplet 
spreading. 

Super  hydrophobic  surfaces,  of  PDMS  material,  consisting  of  regular  array  of  pillars  were 
fabricated  using  soft-lithography  technique.  On  2-D  array  of  posts  and  pores,  a  water  droplet 
placed  clearly  shows  Cassie-Baxter  kind  of  wetting  and  measured  water  contact  angle  value 
also  fairly  matches  with  the  predictions  from  the  Cassie-Baxter  relation.  Water  droplet 
placed  on  these  surfaces,  from  a  hand-held  syringe,  were  associated  with  large  oscillations 
due  to  inherent  inertia  of  the  depositing  and  droplet  snap-off  process.  This  oscillations 
causes  the  droplet  to  overshoot  local  minima  and  arrive  at  an  equilibrium  which  is  closer  to 
the  global  minimum,  i.e.  near  Cassie-Baxter  equilibrium.  However,  when  droplet  is  more 
carefully  placed  and  is  of  higher  viscosity,  inertial  effects  are  quickly  damped  and  droplet 
motion  is  arrested  at  the  first  local  minimum  available.  Dynamics  of  these  viscous  droplets 
studied  with  a  fast  video  camera  reveal  the  same  slip-stick  motion  (in  agreement  with  the 
simulations)  of  the  contact  line  and  droplet  achieves  equilibrium  while  remaining  pinned  at  a 
downstream  end  of  the  post  of  1-D  parallel  grooved  surface.  Observed  apparent  equilibrium 
contact  angle  is  much  larger  than  the  Cassie-Baxter  prediction  as  the  droplet  remains  pinned 
in  the  local  minimum.  Similarly,  controlled  spreading  experiment  reveals  multiple  contact 
angles  provided  by  the  corner  point  before  the  contact  line  jumps  to  the  next  available  pillar. 
And  the  droplet  interface  indeed  makes  a  contact  angle  near  180°  before  it  advances  to  the 
next  pillar.  These  observations  are  consistent  with  our  numerical  simulations  and  reiterate 
the  fact  that  wetting  is  a  dynamic  phenomena  and  equilibrium  droplet  contact  angle  on 
a  textured  surface  is  dependent  upon  the  history  of  the  droplet  as  well  as  available  local 
minima  in  the  form  of  pinning  points. 
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Chapter  3 


Hydrodynamics  of  Water  Droplet 
Movement  in  Oil  over  Particle 
Microtextures  Impregnated  With  Oil 
to  Reduce  Droplet  Adhesion 

3.1  Background 

This  research  is  a  fundamental  study  on  texturing  a  surface,  so  that,  when  the  surface  is 
submerged  in  a  hydrocarbon  (oil)  phase,  the  adhesion  of  water  droplets  in  the  oil  to  the 
surface  is  low,  and  the  mobility  of  the  droplets  along  the  surface  is  high.  The  rationale  for 
the  study  is  based  on  the  lotus  leaf  paradigm  or  super  hydrophobic  effect  (see  the  previous 
chapter  and  the  references 1-19 ),  in  which  water  droplets  in  air  placed  on  a  hydrophobic 
textured  surface  are  observed  to  be  nonadherent,  and  to  easily  roll  off  the  surface.  The 
water  droplet  does  not  penetrate  into  the  texture  due  to  the  hydrophobicity,  and  instead 
beads- up  on  the  surface,  and  under  a  lateral  force  (e.g.  gravity),  glides  along  the  tops  of  the 
texture.  As  reviewed  in  these  references,  techniques  for  the  materials  fabrication  of  surfaces 
with  the  required  hydrophobic  texture  have  been  developed,  and  have  been  applied  to  many 
applications,  in  drag  reduction  and  de-icing,  and  as  anti-sticking,  anti-fouling  surfaces  and 
self-cleaning  surfaces.  Virtually  all  of  the  research  has  centered  on  water  droplets  or  water 
flow  over  textured  surfaces  in  air.  We  study  here  how  this  interfacial  physics  can  be  applied 
to  design  surfaces  which,  under  oil,  allow  water  droplets  to  be  non-sticking  and  interfacially 
mobile.  A  primary  application  of  surfaces  which  are  repellant  to  water  droplets  in  oil,  and 
which  motivates  our  work,  is  to  prevent  corrosion  of  surfaces  in  oils  due  to  water  droplet 
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adhesion  as  for  example  in  oil  pipelinning. 

We  provide  first  a  brief  background  on  the  rationale  for  using  superhydrophobic  surfaces 
under  oil  to  reduce  the  adhesion  of  water  droplets  in  the  oil  to  the  surface,  and  increase  their 
mobility  along  the  surface.  We  then  detail  our  preliminary  progress  on  fabricating  nano 
and  micro-textured  surfaces  through  particle  assembly,  and  our  contact  angle  measurements 
(advancing,  receding  and  static)  which  verify  that  water  droplets  are  non-adhesive  on  these 
surfaces.  Finally  we  detail  our  measurements  of  the  gliding  (rolling-off)  motion  of  the  drops 
under  oil  over  particle-textured  surfaces,  in  which  the  surfaces  are  tilted  to  force  the  droplet 
motion. 


continuous  oil 
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water  drops  stick  on  tilted  surfaces 
due  to  contact  angle  hysteresis 


Figure  3.1:  A  water  droplet  in  a  hydrocarbon  oil,  separated  by  a  thickness  h  from  a  wall  with  a  planar 
surface,  settles  to  the  wall  surface  where  it  sticks  with  a  finite  contact  angle  0pianar  on  a  surface  that  is 
polar  or  slightly  hydrophobic.  The  surface  is  nominally  flat  with  a  surface  roughness.  Due  to  contact  angle 
hysteresis,  the  droplet  remains  stuck  when  the  interface  is  inclined,  or  moves  with  a  retarded  motion. 


For  illustration,  in  Fig.  3.1  is  shown  a  water  droplet  (of  order  1-100  /im  in  (equivalent) 
spherical  diameter),  entrained  in  the  flow  of  a  continuous  oil  phase  in  a  channel  viewed  in 
cross  section.  Buoyancy  forces,  owing  to  density  differences  between  the  water  and  oil  phases, 
bring  the  droplet  in  the  vicinity  of  a  channel  wall.  As  gravity  presses  the  droplet  against  the 
wall,  it  drains  the  oil  in  the  gap  between  the  droplet  and  the  surface,  and  flattens  the  fluid 
interface  to  create  a  thin  intervening  film  of  oil  which  can  rupture  and  adhere  the  droplet 
to  the  surface  (for  reviews  see,  for  example,62-67  for  film  stability  and22,68-70  for  surface 
wetting).  Adhesion  occurs  when  the  surface  is  polar  enough  so  that  ryw/s  <  j0/s  +  7 W/Q 
(where  7 Q/w  is  the  oil/ water  fluid  tension  and  7,,, 7  and  7 a/s  are  the  water/solid  and  oil/solid 
surface  energies)  and  the  balance  of  tensions  at  the  contact  line  insures  a  finite  contact  angle 
(lw/s  —  lo/s  —  lo/w  cos6eq,  where  9eq  is  measured  through  the  water  phase).  This  is  the 


56 


usual  case  in  oil  contacting  equipment,  where  the  interfaces  are  usually  polar,  e.g.  stainless 
steel.  When  the  drop  attaches  with  a  finite  contact  angle  (6*),  the  contact  line  perimeter  can 
immobilize  the  drop  because  of  contact  angle  hysteresis  due  to  surface  roughness 20-22 ,  Fig. 
3.1.  In  order  for  the  drop  to  move  along  the  surface,  the  interface  at  the  advancing  perimeter 
has  to  deform  to  a  larger  contact  angle  then  the  equilibrium,  and  the  trailing  perimeter  to  a 
smaller  angle  (cf.  the  discussion  in  the  previous  chapter).  Continuous  phase  stresses  above 
a  critical  value  can  move  the  droplet  (see  for  example71-73);  however,  even  in  this  case,  the 
drop  mobility  is  retarded  by  this  hysteresis  and  usually  moves  in  a  slow,  episodic,  jerky 
motion. 

As  discussed  earlier,  a  super  hydrophobic  surface  in  air  is  a  textured  interface  in  which 
the  texture  can  be  regular  (e.g.  pillars  or  wells)  or  random  (as  for  example  textures  obtained 
by  etching,  polymer  adsorption  or  particle  coating).  The  texture  is  required  to  be  hydropho¬ 
bic  (either  the  texture  material  is  hydrophobic,  or  the  texture  surface  is  modified  with  a 
hydrophobic  coating) .  Because  of  this  hydrophobicity,  when  a  water  droplet  is  placed  on  top 
of  the  surface  in  air  ,  it  does  not  penetrate  into  the  texture,  but  resides  over  a  cushion  of  air, 
forming  a  contact  line  periphery  on  the  surface  by  pinning  to  the  texture  edges  (“Cassie- 
Baxter”  wetting)1-19.  This  pinning  results  in  a  contact  angle  (#*)  which  is  much  larger  than 
the  contact  angle  ( 9 )  of  a  water  droplet  on  a  smooth  surface  of  the  texture  material  (or  a 
smooth  surface  with  the  same  hydrophobic  coating).  Thus  the  droplet  contact  area  on  the 
super  hydrophobic  surface  is  reduced  relative  to  the  smooth  surface.  A  simple  expression  for 
the  larger  contact  angle,  in  general  agreement  with  data,  is  based  on  minimizing  the  fluid  and 
solid  energies,  and  is  given  by  the  Cassie-Baxter  equation  cos  0*  =  (1  +  cos  9)0  s  —  1  where 
<ps  is  the  area  fraction  of  the  solid.  Hence  for  an  area  fraction  of  <p=.l  and  a  contact  angle 
of  9= 90  °,  for  example,  an  angle  of  #*=150  °  is  subtended  on  the  superhydrophobic  surface. 
A  second  feature  of  the  Cassie-Baxter  wetting  state  is  that  water  droplets  in  air  can  glide 
over  the  surface.  Measurements  of  contact  angle  hysteresis  on  superhydrophobic  surfaces 
show  that  microtextures  of  posts  in  particular  (in  contrast  to  wells)  have  small  contact  angle 
hysteresis  which  tends  to  zero  as  the  area  fraction  tends  to  zero23-28.  In  the  case  of  posts, 
the  advancing  and  receding  menisci  move  similarly  over  the  tops  of  the  posts,  and  hysteretic 
effects  only  arise  from  differences  in  motions  over  the  solid  surfaces,  which  become  reduced 
as  (j)  — >  0.  The  receding  meniscus  over  wells  moves  along  the  matrix  around  the  wells,  while 
the  advancing  meniscus  jumps  over  the  wells  leading  to  a  larger  hysteretic  effect.  Within  the 
drop,  the  viscous  frictional  resistance  as  the  liquid  moves  along  the  surface  is  also  reduced 
because  of  the  air  cushion.  This  motion  is  characterized  by  a  slip  coefficient  A  defined  as 
u  •  t  =  An  -  {Vu  +  Vu*}  -  t  where  u  is  the  velocity,  n  and  t  are  the  unit  normal  and  tan¬ 
gent  vectors  to  the  surface,  Vu  is  the  rate  of  strain  tensor  (evaluated  at  the  surface)  and 
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A  — >  0  describes  the  usual  no-slip  condition.  While  slip  coefficients  are  usually  of  the  scale 
of  nanometers,  for  super  hydrophobic  surfaces  the  slip  coefficient  can  exceed  microns 34,74-76 . 


continuous  oil 
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y 

textured  very  hydrophobic 
surface 


water  drops  stick  on  tilted  surfaces 
due  to  contact  angle  hysteresis 


Figure  3.2:  Superhydrophobicity  and  Cassie-Baxter  Wetting:  (a)  microtextured  hydrophobic  surface  of 
posts  under  oil  is  impregnated  by  the  oil,  and  water  droplet  on  the  surface  subtends  a  large  contact  angle 
because  of  the  pinning  of  the  interface  to  the  geometry  as  is  the  case  in  air.  More  importantly  contact  angle 
hysteresis  is  eliminated  because  of  the  reduced  contact  of  the  water  with  the  solid  surface  of  the  texture  and 
the  drop  glides  over  the  texture  and  does  not  adhere. 

We  extend  the  concepts  of  Cassie-Baxter  wetting  and  superhydrophobicity  to  develop 
surfaces  which  allow  water  droplets  in  oil  and  in  contact  with  the  surface  to  move  more 
smoothly  over  the  interface,  rather  than  sticking  and  (or)  moving  episodically.  In  this  way, 
the  droplets  have  no  residence  time  on  the  surface  of  the  wall  as  they  are  flushed  downstream 
with  the  flow,  and  this  should  prevent  nucleation  phenomena  initiated  by  stagnant  drops 
on  surfaces,  such  as  corrosion.  A  textured  hydrophobic  surface  immersed  in  oil  allows  a 
water  droplet  to  rest  on  the  surface  with  the  oil  impregnating  the  texture  underneath  the 
water,  Fig.  3.2  and  a  large  contact  angle  is  achieved  (6*)  as  the  contact  line  is  pinned 
to  the  texture.  The  static  and  dynamic  behavior  of  such  droplets  have  not  been  studied  as 
extensively  as  the  superhydrophobic  configuration  in  air.  But  most  importantly,  a  significant 
reduction  in  the  contact  angle  hysteresis  relative  to  a  smooth  surface  should  be  realized  (for 
posts  in  particular),  again  because  the  interface  drags  over  the  posts  in  a  similar  manner 
whether  advancing  or  receding  in  oil  or  air.  We  also  note  that,  within  the  droplet,  while  the 
friction  of  the  aqueous  phase  moving  along  the  saturated  microtextured  surface  is  reduced 
because  of  the  liquid  held  in  the  structure,  the  effective  slip  generated  (as  measured  by  the 
slip  coefficient)  is  not  as  great  as  occurs  when  the  texture  is  filled  with  air.  The  principal 
effect  that  the  Cassie-Baxter  wetting  affords  as  drops  move  in  oil  over  the  oil-saturated 
microtexture  is  an  elimination  of  the  contact  angle  hysteresis.  We  note  that  Aizenberg  et  al77 
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and  Varanasi  et  al78  in  designing  alternative  superhydrophobic  surfaces  for  allowing  water 
droplets  to  move  unimpeded  over  a  surface  in  air  have  also  used  hydrophobic  microstructures 
saturated  with  an  oil,  which  in  air  supports  a  thin  film  of  oil  over  the  surface.  Droplets  in 
air  then  ride  over  this  thin  lubricating  layer,  rather  than  on  the  solid  surface;  the  difference 
with  the  study  here  is  that  the  water  droplets  are  completely  submerged  in  oil. 


oil 


(a)  (b) 


Figure  3.3:  (a)  A  mixture  of  nanoparticles  and  binder  cures  to  form  a  microtextured  surface  (b)  Atomic 
force  microscopy  image  of  the  particle  surface  in  tapping  mode. 


3.2  Microtexture  Fabrication  Using  Nanoparticles 

To  fabricate  the  textured  surfaces,  we  disperse  particles  in  a  polymer  resin,  coat  thin  layers 
of  the  mixture  on  a  substrate,  and  then  cure  the  polymer  resin.  If  the  volume  fraction  of 
the  particle  relative  to  the  polymer  is  large  enough,  then  as  the  resin  cures  a  roughened 
surface  forms  from  the  thin  coats  of  the  mixture  which  functions  as  they  have  a  hydrophobic 
texture  (Fig.  3.3(a)).  This  method  for  texturizing  the  surface  is  scalable,  allowing  for  large 
surface  area  modification.  Hydrophobically  modified  spherical  silica  nanoparticles,  20-30 
nm  in  diameter,  were  first  dispersed  in  a  mixture  of  tetrahydrofuran  (THF)  and  isopropyl 
alcohol  (IPA)  (40.60  v/v)  at  a  concentration  of  10  mg/ml  by  sonication.  After  dispersal, 
150  mg  of  methylphenyl  silicone  resin  binder  (Momentive  Performance  Materials,  SR355S) 
was  dispersed  in  the  mixture  and  sonicated.  After  sonication,  the  mixture  was  diluted  with 
fresh  THF/IPA.  To  coat  a  substrate  with  this  mixture,  cleaned  glass  slides  were  vertically 
dipped  and  withdrawn  from  the  mixture  at  a  velocity  of  10  cm/min  using  a  dip-coater  over 
repeated  cycles  to  coat  the  glass  slide  with  the  resin.  After  coating  the  sample  was  put  in  an 
oven  at  40°C  for  10  min  to  cure  the  resin  binder.  The  important  consideration  in  this  study 
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is  that  the  nanoparticle/binder  adhere  strongly  to  the  surface  so  that  the  impregnating  oil 
does  not  remove  them.  For  this  purpose,  typically  six  cycles  of  dip-coating  are  required. 

In  Fig.  3.3(b)  are  shown  atomic  force  microscopy  topology  images  in  tapping  mode  of  the 
surface  made  with  six  dip-coats  after  curing.  The  topology  shows  peaks  and  valleys  of  the 
order  of  a  few  hundred  nanometers,  and  lateral  variations  along  the  surface  of  a  few  microns. 


3.3  Equilibrium  Contact  Angles  in  Air  and  Under  Oil 

Equilibrium  contact  angles  ( Oeqjtexture )  of  water  and  oil  droplets  on  the  fabricated  microtex- 
tured  surfaces  in  air,  and  water  droplets  under  oil  are  undertaken  using  standard  contact 
angle  goniometry  to  measure  the  angle  directly  from  video  images  of  sessile  drops.  Fig. 
3.4(a)  shows  a  side  view  of  a  water  droplet  in  air  on  the  textured  surface;  contact  angles 
equal  to  165°  are  recorded.  When  a  droplet  of  oil  (hexadecane  )  is  placed  on  the  surface, 
the  hydrophobic  nature  of  the  texture  allows  the  oil  to  infuse  into  the  texture.  The  inter¬ 
face  becomes  oleophilic  and  a  contact  angle  of  approximately  8°  is  measured  (Fig.  3.4(b)). 
Finally  a  water  droplet  placed  on  the  surface  under  hexadecane  shows  a  large  contact  angle, 
also  approximately  160°  (Fig.  3.4(c)) 


(a)  (b)  (c) 


Figure  3.4:  (a)  A  water  droplet  in  air  on  the  nanoparticle  microtexture,  (b)  a  hexadecane  oil  droplet  in 
air  on  the  microtexture,  and  (c)  an  equilibrium  water  droplet  under  hexadecane  on  the  texture. 


The  low  contact  angle  of  the  oil  droplet  on  the  surface  in  air  suggests  the  oil  is  wetting 
into  the  microstructure,  and  the  large  contact  angles  of  the  water  droplets  in  oil  indicate  a 
Cassie-Baxter  wetting  state  since  the  contact  angles  of  water  on  methyl  silanated  surfaces  are 
typically  140°  in  oil.  We  have  also  measured  the  contact  angle  hysteresis  of  water  droplets 
on  these  surfaces  and  find  them  only  to  be  a  few  degrees,  also  confirming  a  Cassie-Baxter 
wetting  state,  in  contrast  to  water  drops  in  oil  on  a  planar  silanated  surface  of  approximately 
15  degrees.  A  clear  demonstration  of  the  mobility  of  the  nanoparticle  textured  surface  to 
water  droplets  in  oil  is  shown  in  Fig.  3.5,  in  which  a  water  droplet  is  released  onto  an  inclined 
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surface  and  easily  rolls  on  the  surface. 


(a)  (b)  (c) 


Figure  3.5:  Movement  of  a  water  droplet  under  hexadecane  along  a  surface  inclined  a  few  degrees.  The 
droplet  is  a  few  hundred  microns  in  diameter 
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Chapter  4 


Hydrodynamics  of  Large  Interfacial 
Slip  at  Porous  (Air-Permeable) 
Polymer  Surfaces  For  Drag  Reduction 

4.1  Introduction 

Recently,  new  attention79  has  been  paid  to  the  possibility  of  hydrodynamic  slip  at  an  in¬ 
terface  between  a  stationary  solid  surface  and  a  “simple”  (non  polymeric)  liquid  moving 
over  the  surface.  The  consideration  of  boundary  slip  began  long  ago  with  the  continuum 
level  formulation  of  the  Navier  slip  condition  for  a  Newtonian  fluid  which  equated  the  fluid 
velocity  tangent  to  the  surface,  vs,  to  the  boundary  tangential  stress,  rs,  by  the  slip  coef¬ 
ficient  A,  i.e.  vs  =  jTs  where  ji  is  the  fluid  viscosity  and  A  has  units  of  length.  Since  this 
formulation,  experimental  studies  have  made  clear  that  the  “no-slip”  condition  of  A  =  0  is 
sufficient  to  accurately  model  most  macroscopic  flows  with  length  scales  in  the  range  of  mil¬ 
limeters  to  meters.  However,  recent  molecular  dynamics  simulations  on  atomically  smooth 
surfaces  have  demonstrated  slip  on  the  molecular  scale,  and  calculated  A  as  a  function  of 
the  strength  of  the  liquid-solid  interaction80,81.  For  strong  liquid-solid  interactions  which 
characterize  complete  or  strong  wetting  of  the  liquid  on  the  solid  surface,  slip  lengths  are  of 
the  order  of  only  a  few  molecular  diameters  (0(1  nm)),  while  relatively  weaker  interactions 
of  partially  wetting  fluids  have  slip  lengths  extending  tens  of  diameters  (0(10  nm)).  Current 
tools  for  measuring  A  include  particle  image  velocimetry  (PIV),  image  velocimetry  enhanced 
with  evanescent  near-field  illumination  at  the  surface,  and  atomic  force  microscope  (AFM) 
and  surface  force  apparatus  (SF)  measurements. 82  Several  experiments  are  consistent  with 
the  MD  calculations.  For  smooth  surfaces  of  a  wetting  liquid,  e.g.  water  flowing  along  a 
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hydrophilic  surface,  either  zero  slip  or  coefficients  less  than  a  few  nanometers  are  recorded, 
while  water  over  a  partially  wetting  hydrophobic  surface  (e.g.  self  assembled  octadecyl  silane 
(OTS)  monolayers)  obtain  slip  lengths  of  tens  of  nanometers,  e.g.83-86.  For  nonpolar  liquids 
wetting  smooth  hydrophobic  surfaces  (the  weakest  liquid/surface  interactions),  slip  lengths 
are  of  the  order  of  a  few  tens  of  nanometers87-89. 

“Giant”  slip  (A  order  1  /mi  or  larger)  is  the  subject  of  great  interest  for  its  application  to 
reducing  surface  friction  in  micro  and  nanofluidic  channel  flows.  Large  slip  can  be  achieved 
when  low  friction  air  layers  are  situated  between  the  liquid  and  the  surface  a  circumstance 
which  arises  when  a  population  of  nanoscopic  gaseous  domains  adhere  to  a  surface,  or  air 
becomes  trapped  in  a  micro  or  nano-textured  surface  that  is  not  wet  by  the  liquid  (superhy- 
drophobicity) . 90  AFM  studies  have  provided  direct  evidence  of  nanoscopic  gaseous  domains, 
primarily  at  the  interface  between  water  and  a  hydrophobic  surface.91  AFM  measurements 
have  also  verified  the  reduction  in  surface  friction92,  and  nanoscopic  gas  domains  have  been 
suggested93  as  one  reason  why  some  measurements  of  slip  at  the  water /hydrophobic  solid 
surface94,95  obtain  one  micron  or  larger  slip  lengths.  Hydrophobic,  textured  surfaces  filled 
with  air  provide  a  more  reproducible  method  for  generating  large  slip,96-99,  and  theoretical 
MD  and  continuum  studies  of  model  textures  (see  for  example100,101)  demonstrate  that  large 
slip  requires  the  solid  fraction  of  the  surface  to  be  a  few  percent. 

Little  attention  has  been  paid  to  obtaining  giant  slip  when  nonpolar  liquids  slip  over  a 
surface  ( “oil  slippery  surfaces” ) ,  although  significant  interest  is  developing  due  to  the  emer¬ 
gence  of  dropwise  microfluidic  platforms  which  are  based  on  water  droplets  moving  in  a 
continuous  oil  stream102.  Generating  gaseous  domains  at  the  surface  is  the  key  to  large  slip, 
and  one  method  for  generating  a  significant  coverage  of  nanobubbles  on  surfaces  submerged 
in  oil  is  to  nucleate  them  spontaneously  on  oil  contact.  Spontaneous  nonequilibrium  for¬ 
mation  of  nanobubbles  at  a  surface  occurs,  for  example,  when  an  air  saturated  liquid  (e.g. 
ethanol)  is  displaced  by  a  second,  miscible,  liquid  (e.g.  water),  also  saturated  with  air,  but 
with  a  lower  solubility91  and  nanobubbles  nucleate  to  accommodate  the  reduced  solubility 
of  the  displacing  phase.  In  this  letter  we  demonstrate  that  a  nanoporous,  hydrophobic  poly¬ 
mer,  polydimethylsiloxane  (PDMS),  which  has  a  significant  permeability  to  air  because  of 
the  hydrophobicity  of  the  polymer  and  a  relatively  large  free  volume103,  permits  oil  to  slip 
over  its  surface  with  an  order  one  micron  slip  length.  When  oil  contacts  the  surface  of  an 
air-equilibrated  PDMS  substrate,  air  is  released  as  nanobubbles  to  the  oil/PDMS  surface. 
The  driving  force  for  this  release  derives  from  the  fact  that  the  PDMS  is  also  permeable  to 
oil  because  of  its  hydrophobic  nanoporosity,  and  diffusion  of  oil  into  the  nanopores  on  the 


63 


Figure  4.1:  Measurement  of  microchannel  slip  at  an  oil/PDMS  surface  by  observing  the  dielectrophoretic 
merging  of  water  droplets  in  oil  moving  in  close  proximity  to  the  PDMS  channel  wall. 


contact  of  oil  to  the  surface  displaces  the  air  to  the  interface. 

We  demonstrate  in  particular  oil  slip  at  the  inside  surface  of  a  microchannel  formed  in  a 
PDMS  monolith  fabricated  by  soft  lithography.  We  use  a  viscous  mineral  oil  (/iG^  approxi¬ 
mately  102  times  that  of  water),  a  mix  of  high  and  low  molecular  weight  olefins,  which  does 
not  macroscopically  swell  the  PDMS104  and  distort  the  surface  from  the  atomically  smooth 
topology  normally  evident  in  AFM  measurements  but  can  still  displace  air  in  PDMS  by  the 
solubilization  of  small  linear  alkanes  in  the  oil. 


4.2  Experimental  Measurement  of  Droplet  Approach 

The  microfluidic  arrangement  is  shown  in  Fig.  4.1.  (Additional  details  are  given  in  the  Ap¬ 
pendix  A.)  Using  optical  microscopy,  we  measure,  in  a  PDMS  channel  of  rectangular  cross 
section  with  height  h  (100  /.an)  and  wide  width  w  (300  /.an) ,  the  edge-to-edge  separation  dis¬ 
tance  (s(t))  of  pairs  of  nearly  occluding  water  droplets  (“A”  and  “B”)  which  are  entrained 
in  a  mineral  oil  stream  and  are  driven  together  by  a  dielectrophoretic  (DEP)  force  of  attrac¬ 
tion,  FDEp.  This  force  is  due  to  an  electric  field  £  applied  parallel  to  the  bottom  wall  of 
the  channel  and  along  the  flow  direction.  The  merging  droplet  pair  is  part  of  a  single  file 
droplet  train,  formed  upstream  by  flow  focusing105  of  oil  and  water  streams  (flow  rates  equal 
to  0.4-1  jd / min  and  0.04-0.1  j)l / min,  respectively),  through  an  orifice  50  //in  in  width.  The 


64 


focusing  forms  droplets  with  radius  a  equal  to  approximately  40  //m,  which  are  separated  by 

a  few  radii,  and  flow  at  the  average  stream  velocity  V  of  approximately  250  -  500  //in /sec. 

The  electric  field  is  applied  to  switchback  flow  lanes  at  the  downstream  end  of  the  chip  by  a 

voltage  V  set  across  parallel  copper  strip  electrodes  inserted  through  the  PDMS  (a  dielectric) 

to  insure  a  uniform  field  across  the  flow  channels.  Relative  to  the  aqueous  phase  which  is 

de-ionized  water,  the  oil  is  nonconducting  and  the  field  polarizes  each  of  the  droplets  of  the 

pair  into  dipoles.  In  the  flow  lanes  parallel  to  £,  the  polarized  drops  are  aligned  with  £ 

and  attract  each  other,  creating  the  dielectrophoretic  force.  Prior  to  the  application  of  the 

field,  the  train  flow  is  observed,  and  when  a  pair  are  observed  to  pair-off  to  a  relatively  close 

separation  (less  than  one  radius)  due  to  flow  disturbances,  the  field  is  applied  to  merge  the 

pair,  and  a  video  recording  is  made  with  a  high  speed  camera.  At  the  time  of  application 

of  the  field,  the  droplets,  heavier  than  the  oil  (pW)ater=103Kg/m3  and  p0u= 8.75xl02Kg/m3) 

have  settled  to  a  separation  distance  d  from  the  bottom  wall  of  the  microchannel  which  is 

determined  by  their  settling  velocity  and  the  transit  time,  r,  from  their  formation  at  the 

orifice  until  application  of  the  field  and  is  of  the  order  of  a  few  hundred  nanometers.  At 

ds 

these  distances  from  the  wall,  the  approach  velocity  Va  —  Vb  =  ~r  due  to  the  DEP  force 

dt 

is  affected  (as  we  show  by  numerical  solution  of  the  hydrodynamic  equations)  by  the  drag, 
against  the  bottom  wall,  of  the  intervening  oil  between  the  droplet  and  the  wall,  and  the 
slip  on  this  wall.  From  comparison  of  s(t)  with  numerical  solutions,  A  is  obtained  with  a 
precision  which  can  distinguish  a  micron  size  slip  length.  We  also  measure  the  slip  when  the 
bottom  surface  is  a  glass  slide,  functionalized  with  octadecyltrichlorosilane  (OTS),  which  we 
do  not  expect  to  nucleate  nanobubbles  at  the  surface  because  of  low  air  permeability,  and 
we  find  zero  slip  to  the  precision  of  the  measurement. 

One  illustrative  data  set  is  given  in  Fig.  4.2  An  edge  detection  routine  is  used  to  de¬ 
termine  the  droplet  perimeters  from  which  the  size  of  the  merging  droplets  a  and  the  pair 
separation  distance  s(t)  are  computed.  The  separation  distance  is  shown  in  Fig.  4.2,  and  the 
relative  approach  velocity  is  approximately  40  //m/sec.  To  obtain  different  data  sets  of  s(t) 
corresponding  to  different  droplet  radii  or  droplet- wall  separations  d  (transit  times  r),  the 
oil  and  water  flow  rates  at  the  flow  focusing  orifice  are  changed,  or  the  merging  at  different 
switchback  lanes  is  observed. 

A  nonionic  surfactant  (Span  80,  sorbiton  monoleate)  is  dissolved  in  the  mineral  oil  at 
a  concentration  C  =  2.3 xlO-2  M,  well  above  its  critical  micelle  concentration  (CMC) 
(Ccmc— 2.3  xlO-4  M).  At  and  above  the  CMC  the  equilibrium  tension  7  is  3  mN/m  and  the 
equilibrium  surface  concentration  r<7MC=3.6x  I0_6mole/m2;  dynamic  tension  measurements 
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Figure  4.2:  Dielectrophoretic  merging  of  40  pm  radius  droplets  at  the  oil/PDMS  surface:  Frame  captures 
of  the  pairwise  merging  at  time  intervals  of  .12  sec,  flow  direction  from  bottom  to  top  and  time  (t)  as  function 
of  the  measured  edge-to-edge  scaled  separation  s/a  from  the  images.  The  continuous  line  is  a  fit  for  a  value 
for  the  droplet-wall  drag  coefficient,  a. 


indicate  the  desorption  rate  constant,  kd ,  is  of  order  10-3  sec-1.106  Since  the  time  scale  for 
convection  of  surfactant  along  the  droplet  a/V,  order  10-1  sec,  is  much  shorter  than  the  des¬ 
orption  time,  1/kd,  surfactant  collects  at  the  trailing  edge  of  the  droplet  causing  the  tension 
to  be  larger  at  the  front  than  the  back.  This  Marangoni  gradient  opposes  the  surface  flow  and 
immobilizes  the  interface  since  the  ratio  of  the  characteristic  scale  for  the  retarding  tension 
gradient  RTTcmc / a  to  the  oil  viscous  stress  on  the  droplet  surface,  ji.0,i  V /a,  the  Marangoni 


number  (Ma  = 


RTF, 


CMC < 


VoilV 

of  interfacially  rigid  droplets. 


is  order  102.  Hence  the  pairwise  hydrodynamic  interaction  is  one 


To  compare  the  data  sets  of  s(t)  to  a  hydrodynamic  model  of  the  merging  process,  the 
applied  field  and  droplet-wall  separation  distance  d  have  to  be  determined.  Since  the 
PDMS  and  oil  are  dielectric  phases,  they  act  as  capacitors  in  series,  and  therefore  £  = 
V{2 LPDms{£ou/ £pDMs)  +  LcY1  with  Spdms  and  s0u  the  dielectric  constants  for  the  PDMS 
and  oil  (2.65  and  2.18,  respectively),  and  Lc  and  Lpdms  are,  respectively,  the  length  of  the 
channel  and  the  distance  between  the  electrode  and  the  channel  (6  mm  and  4.5  mm,  re¬ 
spectively).  An  AC  electric  potential  V  (500  Hz,  5  kV  sine)  is  applied  across  the  electrodes, 
resulting  in  an  average  field  strength  £=365  V/mm.  While  an  alternating  potential  is  used  to 
prevent  any  residual  charge  accumulation  in  the  PDMS  and  oil,  the  oscillation  does  not  affect 
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the  merging  process  since  FDEP  is  proportional  to  the  square  of  the  field,  and  the  oscillation 
period  is  much  faster  than  than  time  scale  for  the  merging.  To  obtain  d,  we  assume  that  at 
the  flow-focusing  orifice  the  droplets  detach  symmetrically  from  the  top  and  bottom  walls  of 
the  channel  and  are  therefore  initially  centered  at  the  midplane,  hence  the  initial  separation 
di  is  h/2  —  a.  Clearly,  each  drop  is  not  precisely  released  at  the  center  of  the  channel,  but 
statistically  above  and  below  the  midplane.  However,  our  data  analysis  will  average  over 
several  data  sets  at  a  nominal  value  of  d  which  should  account  for  this  statistical  variation. 
The  distance  d,  accounting  for  only  the  resistance  of  the  lower  wall,  is,  for  d/a  <C  1,  given 
by107  hi  |  y:  |  =  —  ^sahv^er^Poii) T  where  g  js  the  acceleration  of  gravity.  (Values  for  d/a  are 
less  than  0.25,  for  which  this  expression  is  accurate.  The  equation  for  d  can  be  corrected  to 
include  wall-slip  and  the  effect  of  the  resistance  due  to  the  opposite  channel  wall,  but  these 
corrections,  if  included,  can  be  shown  to  be  negligible  in  the  determination  of  d. 108 ) 


4.3  Finite  Element  Simulations  of  the  Droplet  Hydro¬ 
dynamics 

The  droplet  hydrodynamics  is  in  the  Stokes  regime  of  negligible  inertia.  The  Reynolds  num¬ 
ber,  as  given  by  Re=poiiVh/ /j,oii=O(10~4),  and  the  droplets  remain  spherical  until  the  onset 
of  coalescence  (cf.  Fig.  4.2)  as  the  viscous  forces  are  smaller  than  the  tension  force  (capillary 
number,  Ca  =  [x0i\V / 7=  O(10-2))  and  the  Maxwell  electrical  stresses  are  smaller  than  the 
tension  force  (electric  Bond  number,  Be  =  £0£0uE2a/ 7=O(10-2)),  where  za  is  the  permit¬ 
tivity  of  free  space.  The  total  fluid  drag  exerted  on  each  of  the  droplets  as  they  merge  can 
be  calculated  from  the  sum  of  the  fluid  drags  (formulated  as  a  drag  coefficient  /  multiplied 
by  67 Tfj,0ua  and  a  velocity)  in  three  flow  configurations  (Fig.  4.3):  the  Poiseuille  flow  over 
fixed  droplets  (fp),  and  the  motions,  stationary  in  the  farfield,  of  A  with  B  fixed  or  B  with 
A  fixed,  with  the  later  two  configurations  each  divided  into  a  mutual  approach  (fm)  and 


2  Fdep 


.  The 


an  in-tandem  motion  (fu).  The  total  fluid  drag  balances  Fdep,  Va  —  VE  = 

67 qu0a/m 

slip  is  obtained  by  comparison  of  a  theoretical  calculation  of  fm  (a  function  of  h/a ,  s/a , 
d/a  and  A/a)  to  an  experimental  value  calculated  through  the  measurement  of  V4  —  VE  (or 
equivalently  s(f))  and  Fdep- 


To  obtain  Fdep,  the  electric  field  in  the  mineral  oil  around  the  merging  (uncharged)  water 
droplets  in  the  microchannel  is  approximated  by  the  bispherical  harmonic  solution109  for 
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Figure  4.3:  The  drag  coefficient  fm  as  a  function  of  s/a  for  no-slip  (A=0)  and  A=1  /jm  for  d= 1  pm  (top) 
and  d=100  nm  (bottom)  for  h=100  pm  and  a  =  40  pm. 


the  electrostatic  field  in  an  unbounded  insulating  dielectric  surrounding  a  pair  of  perfectly 
conducting  spheres  (zero  net  charge)  due  to  a  farfield  electric  field  applied  in  the  direction  of 
the  line  of  centers  between  the  spheres,  from  which  FDEP  =  £0£0na28,2^s(s/a).  Q(s/a)  is  an 
infinite  series  function  of  s/a  evaluated  here  with  polynomial  interpolation  (see  Appendix  B). 


The  coefficient  fm  is  calculated  numerically  from  a  COMSOL  (4.2a)  finite  element  simu¬ 
lation  using  the  experimental  parameters.  The  simulation  is  first  verified  by  computing  the 
drag  coefficient  for  a  single  sphere  in  a  plane-parallel  channel  moving  parallel  to  the  wall, 
denoted  as  a(d/a,X/a,h/a),  and  comparing  for  A/a=0,  to  an  interpolating  formulae  for  a 
for  no  slip  obtained  from  a  multipole  solution  by  Feuillebois  et  al110  (see  Appendix  B).  Fig. 
4.3  shows  fm  (symbols)  as  a  function  of  s/a  for  no  slip  and  A=1  //in,  for  a=40  //m  and  for 
a  separation  distance  d  =  1  /im  (top)  and  100  nm  (bottom).  As  expected,  the  closer  the 
droplet  pair  to  the  wall,  the  greater  is  the  influence  of  the  slip,  and  it  is  clear  that  when  d 
of  the  order  of  a  few  hundred  nanometers  a  order  ten  percent  reduction  in  fm  is  achieved 
from  the  no-slip  case  for  A=1  /im,  and  this  change  is  the  basis  of  our  measurement  of  a 
order  one  micron  slip  length.  A  slip  coefficient  can  be  obtained  from  the  experimental  pro¬ 
files  s(t)  by  comparison  to  the  integration  of  —  =  Va  —  VE  = - PE^  using  the  numerical 

dt  67T/i0a/m 

calculations  of  fm  as  a  function  of  s.  We  avoid  this  extended  calculation  by  an  approxi¬ 
mation  for  fm  as  the  sum  of  the  drag  on  a  single  (rigid)  droplet  moving  at  a  distance  d  in 


68 


a  channel  of  height  h  with  slip  A,  a(d/a,  X/a,h/a)  (obtained  by  COMSOL  calculation  and 
independent  of  s),  the  drag  on  a  (rigid)  droplet  pair  mutually  approaching  at  a  distance  s 
from  each  other  in  an  infinite  medium,  (k  (Jeffreys  solution  for  which  there  is  a  correlation, 
IR(^)  =  {l  +  jl  +  .38e_^n“+'68^  /6'3j,  see  supplementary  material)  and  a  correction 

+  j  4*  D  j  — j.  The  correction  factor  T> (s/a)  accounts 

for  the  double  counting  evident  in  f£PProx  when  s/a  tends  to  infinity  and  should  tend  to  -1; 
we  use  D(s/a)  =  —  |  ////,  which  allows  congruence  of  the  approximate  formulation  over  the 
entire  range  of  s/a,  cf.  Fig.  4.3.  A  theoretical  prediction  for  s(t)  can  be  constructed  by 

4a,-|-KR|  -j+£>j  a 

e0eoii£23{f  }/37rpoii 

the  initial  separation  of  the  droplet  pair.  This  prediction  is  easily  fit  to  a  data  set  by  adjust¬ 
ing  a  (which  is  independent  of  s),  and  the  fit  is  shown  for  the  illustrative  data  in  Fig.  4.2. 
(In  practice,  to  include  the  forces  on  the  droplets  due  to  immediate  neighbors  of  the  train 
at  an  assumed  distance  st  from  the  pair,  we  have  addended  their  dielectrophoretic  dipole 
contribution  ( ,  ttta  )  1°  ^  arid  the  leading  order  droplet- droplet  interaction  1  +  a/ (2s/)  to 
(k.)  All  data  sets  corresponding  to  different  d  and  a  are  fit  in  this  way.  To  correlate  the  fitted 
values  of  a  to  A,  we  first  bin  all  the  data  sets  into  groups  in  which  in  each  group  the  radii 
differ  by  at  most  3  percent.  Each  binned  group  in  diameter  is  then  further  binned  into  wall 
separation  distances  d  which  differ  by  no  more  than  3  percent.  The  results  are  plotted  as 
the  symbols  with  error  bars  (from  the  standard  deviation  from  the  average)  in  Fig.  4.4  for  a 
droplet  radii  bin  a=37.5±  1.25  /im  and  a  bin  a=40±1.25  /urn  as  a  function  of  the  (binned) 
values  of  d/a.  As  each  droplet  radius  bin  corresponds  to  a  fixed  value  of  h/a,  the  theoretical 
value  of  a  in  these  bins  is  only  a  function  of  d/a  and  A/a.  Plotted  in  Fig.  5.4  as  symbols 
(for  the  two  radii  bins)  are  the  theoretical  values  as  a  function  of  d/a  for  values  of  A  equal 
to  zero  (the  Feuillebois  correlation110),  1  and  4  [im.  (from  COMSOL  calculation).  For  the 
PDMS  bottom  channel  wall,  the  comparison  of  the  theoretical  and  experimental  values  of  a 
show  clearly  a  micron-sized  slip,  while  for  a  bottom  microchannel  wall  made  of  glass,  no-slip 
is  obtained  as  expected. 
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4.4  Conclusions 

Our  demonstration  of  order  one  micron  slip  at  the  interface  of  an  oil  and  a  polymeric  surface 
(PDMS)  which  releases  to  the  surface  air  retained  in  the  material  to  form  a  lubricating 
layer,  may  serve  as  method  for  enabling  giant  slip  without  having  to  modify  the  surface  with 
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5 


Normalized  Sphere-Wall  Separation  Distance  (d/a) 


(a) 


(b) 

Figure  4.4:  Hydrodynamic  drag  coefficient  a  as  a  function  of  separation  d/a  for  droplets  of  radius  37.5 
(i m  (a)  and  40  /am  (b)  for  fixed  height  h=100  /rm.  Symbols  with  error  bars  are  the  experimentally  fitted 
coefficients,  the  remaining  symbols  are  from  numerical  simulation  (the  accompanying  dotted  lines  are  a 
guide),  and  the  no-slip  line  is  from  the  Feuillebois  et  al110  correlation. 
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an  air-sequestering  texture.  These  results  are  particularly  relevant  to  microfluidics  where 
PDMS  is  the  standard  material,  and  the  use  of  oil  streams  with  reagent  water  droplets  have 
become  a  dominant  lab  on  a  chip  platform. 

4.5  Appendix  A:  Materials  and  Cell  Fabrication 

Aqueous  Phase:  DI  water,  from  Millipore  ultrafiltration  unit,  conductivity,  18  MO  cm-1, 
and  from  handbook  values  pw=  103  kg  m-3,  and  pw=  10-3  kg  m-1s-1  at  20  °C.  Mineral 
oil:  Sigma- Aldrich,  eQu  =  2.18  (parallel  plate  capacitor  measurement),  p0u=  8.75xl02  kg 
m-3  (manufacturer’s  data  sheet)  at  20  °C.  The  oil  viscosity  poli  is  temperature  sensitive  in 
the  range  in  which  the  experiments  were  undertaken  (20-23°C),  and  was  measured  (with  the 
surfactant,  Span  80,  C=2.3x  10-2  M)  by  a  temperature-controlled  rheometer  (see  data  below, 
Figure  1)  and  fit  to  a  correlation  p(T)  =  2.20  x  10_5T3  —  1.63  x  10_3T2  +  3.58  x  10_2T  —  0.139. 
A  thermocouple  in  the  PDMS  measured  the  temperature  for  each  experiment,  and  in  the 
fitting  for  the  drag  coefficient  a  the  viscosity  at  the  temperature  measured  for  the  experiment 
was  used  by  calculation  from  the  correlation. 


Figure  4.5:  Viscosity  of  mineral  oil  as  a  function  of  temperature. 


The  microfluidic  cell  is  fabricated  via  the  methods  of  soft  lithography111  from  two  slabs 
of  polymerized  and  cured  PDMS,  one  slab  containing  the  inscribed  channels  and  a  second, 
flat  slab  bonded  to  the  first  to  enclose  the  channels.  The  PDMS  prepolymer  and  curing 
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agent  is  obtained  as  a  kit,  Sylgard  184  (Dow  Corning).  The  top  half  of  the  cell  containing 
the  inscribed  fluidic  circuitry  is  molded  from  a  master  fabricated  by  photolithography  using 
a  transparency  mask  and  a  negative  tone  SU-8  photoresist  and  developer  (Microchem),  and 
the  lower  half  is  a  cast  layer  of  PDMS  which  is  bonded  to  the  top  half  by  plasma  treating  the 
surfaces  and  then  conformally  contacting  the  halves.  Dielectric  constant  of  cured  elastomer 
is  2.65  (manufacturer’s  material  data  sheet).  Fluids  are  delivered  to  the  cell  from  syringe 
pumps  (Harvard  PHD)  via  polyethylene  tubing  (1.5  mm  ID)  through  punch-hole  entry  ports. 

4.6  Appendix  B:  Interpolating  Functions  and  Valida¬ 
tion  Simulations  of  the  Finite  Element  Calculations 

The  calculation  of  Davis109  for  the  dielectric  force  on  a  pair  of  zero  charge  perfectly  conduct¬ 
ing  spheres  in  an  unbounded  insulating  dielectric  subject  to  a  uniform  electric  field  £  at  infin¬ 
ity  applied  along  the  line  of  centers  of  the  spheres  is  given  as  an  infinite  series  solution  of  the 
edge-to-edge  separation  distance  s  divided  by  a,  from  which  the  dielectrophoretic  force  be¬ 
tween  the  spheres  is  obtained  as  FDEP  =  e0£0ji£2a29:(s/a)  where  ^s(s/a )  is  an  infinite  series. 

We  approximate  the  series  by  the  polynomial  S(s/a)  =  —2.02  x  10_7(s/a)3  +  7.33  x  10_7(s/a)2  —  9.99  x  1( 
and  a  comparison  of  this  interpolation  formulae  and  the  Davis  expression  is  given  in  4.6. 


Figure  4.6:  Normalized  electrostatic  force  as  a  function  of  normalized  separation  distance  as  given  by  the 
exact  bispherical  calculation  and  the  interpolating  equation. 
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The  solution  for  the  hydrodynamic  drag  exerted  on  a  pair  of  two  mutually  approaching 
(interfacially  rigid  or  solid)  spheres  in  an  infinite  medium  is  the  Stimson  and  Jeffrey  solu¬ 
tion  in  bispherical  coordinates  and  is  expressed  as  an  infinite  series  and  the  drag  coefficient 
fR(s/a)  is  expressed  as  an  infinite  series  (see  eq.  4.1);  for  numerical  calculation  we  use  the 
interpolation  as  given  by  Ivanov  et  al112  and  written  in  the  manuscript 


fR(s/a)  = 


,  n(n  +  1)  (Anexp(2e)  Anexp(— 2e)  exp(— 2eAn)  1 

smhe£ -V^  A'  +  -XT — 1  +  a  '  - 1 

L  ZAn- 1  Z^n+ 1  ^n-l^n+1  J 


S  1 

cosh  e  —  1  +  — ,  A n  =  n  +  -,  An  =  sinh(2Ane)  —  Xn  sin(2e) 

A  CL  A 


(4.1) 


Detailed  multipole  solutions  for  the  drag  coefficient  for  a  singie  solid  sphere  translating 
parallel  and  between  two  parallel  walls  was  given  by  Feuillebois  et  al,  as  cited  in  the  chapter. 
They  construct  an  interpolating  formulae  (eq.  29)  for  this  drag  coefficient,  and  in  Fig.  3  we 
compare  the  resuits  of  this  correlation  with  our  COMSOL  results  as  a  means  of  validating 
our  numerical  simulations. 
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Figure  4.7:  Hydrodynamic  drag  force  on  a  single  sphere  translating  between  two  parallel  walls  for  no  slip, 
A  =  0  fim  as  a  function  of  the  sphere/wall  separation  distance  d  (edge-to-edge) ,  comparing  the  Feuillebois 
correlation  and  our  COMSOL  calculation  for  two  values  of  the  channel  height  h  relative  to  the  sphere  radius 
a. 
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Chapter  5 


Hydrodynamic  Capture  of 
Microbeads  in  Microtexture  Obstacle 
Courses  for  Biomolecular  Screening  in 
Pathogen  Sensing  and  Point  of  Care 
Diagnostics 

5.1  Introduction 

Analytical  methods  which  allow  a  target  protein  to  be  screened  against  a  probe  library  of 
potential  binding  partners  (e.g.  other  proteins  or  small  biomolecular  ligands)  to  determine 
binding  affinities  are  central  to  fundamental  studies  in  cell  and  molecular  biology,  where 
they  can  be  used  to  map  the  web  of  interactions  by  which  proteins  orchestrate  biological 
activity.  Screening  tools  are  also  important  in  applied  research,  where  they  are  used  in  drug, 
antibody  and  enzyme  discovery113-115,  in  the  identification  of  disease  markers  and  in  the 
development  of  biosensors  for  environmental  surveillance  and  food  monitoring116-119.  The 
standard  method  for  screening  protein  binding  interactions  is  the  microtitre  well  plate,  in 
which  the  library  of  probe  molecules  is  spotted  into  an  array  of  wells,  and  each  well  is 
then  incubated  with  a  target  and  interrogated  individually  (typically  with  an  enzyme  linked 
sandwich  immunoassay)  to  determine  target-probe  conjugation.  Many  applications,  partic¬ 
ularly  those  using  combinatorial  synthesis  for  the  development  of  new  biomolecules,  require 
a  large  number  of  probe  molecules  to  be  displayed,  and  for  these  applications  the  screen¬ 
ing  platform  should  be  miniaturized  to  enable  higher  library  probe  densities  than  can  be 
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accommodated  with  microtitre  plates.  Flat  microarrays  (see  for  example  the  reviews120-129) 
are  one  innovation,  in  which  probes  are  robotically  spotted  into  small  areas  or  patches  (  a 
few  hundreds  of  microns  in  diameter)  on  a  substrate  to  form  a  library.  Probe  biomolecules 
either  adsorb  directly  to  the  substrate  surface,  or  the  surface  is  functionalized  to  covalently 
bind  the  printed  probe  molecules.  The  library  is  then  incubated  with  a  target,  and  bind¬ 
ing  is  detected,  in  a  single  (parallel)  step,  usually  by  fluorescently  labelling  the  target  and 
examining  the  array  for  spots  which  fluoresce  upon  excitation.  In  the  flat  microarray  (as 
with  the  microtitre  plate  assay)  target  first  diffuses  through  a  quiescent  analyte  solution 
to  the  surface  and  subsequently  kinetically  binds  to  the  surface  probe.  Binding  is  detected 
when  enough  target  binds  to  the  probes  to  register  a  fluorescent  signal.  When  transport 
through  the  analyte  phase  is  only  by  diffusion,  this  assay  time  can  be  relatively  long  even  if 
the  kinetic  step  is  fast130,131.  Streaming  analyte  across  the  flat  microarray  can  enhance  the 
transport  by  bringing  analyte  solution  directly  to  the  probe  surface  in  a  convective  boundary 
layer.  The  detection  time  would  then  be  determined  by  only  the  kinetic  rate  (a  property 
of  the  probe  and  target),  and  would  not  be  limited  by  a  diffusion  barrier.  With  this  as  a 
goal,  biosensor  designs  have  incorporated  a  printed  microarray  of  probes  (or  a  single  patch  of 
one  probe)  as  the  bottom  surface  of  a  microfluidic  channel,  and  have  convectively  screened 
the  array  with  target  (for  reviews  see132-140).  In  an  alternative  format,  separate  channels 
in  a  microfluidic  cell,  each  displaying  one  probe  at  the  channel  bottom,  have  been  used  to 
construct  a  library141-144.  The  use  of  microfluidics  to  convectively  screen  patch  arrays  of 
probe  libraries  as  a  lab  on  a  chip  not  only  reduces  the  assay  time  for  screening  the  library 
by  enhancing  transport  of  the  target  to  the  probe  (and  thereby  increases  the  throughput), 
it  has  the  added  advantage  of  reducing  the  volumes  of  targets  and  reagents  that  are  needed 
for  the  assay. 

The  microfluidic  screening  platform  in  which  the  probe  library  is  an  array  of  patches 
situated  on  the  floor  of  the  flow  channel  is  one  miniaturized  design.  However,  the  assembly 
of  this  platform  is  complicated  because  the  microarray  has  to  first  be  constructed  ex-situ  of 
the  microfluidic  cell,  and  then  incorporated  into  the  cell,  a  procedure  which  becomes  very 
difficult  when  the  spotted  probes  need  to  be  continually  immersed  in  water,  as  for  example 
cell  membrane  proteins  which  are  spotted  sequestered  in  lipid  bilayers.  In-situ  methods  for 
assembling  the  probe  library  directly  inside  the  flow  channel  are  under  current  study,  and 
are  based  on  using  microbeads  (10-100  /rm  in  diameter)  to  host  the  probe  biomolecules,  and 
facilitate  their  arrangement  into  an  array  library  within  the  microfluidic  cell145-148.  Probe 
biomolecules  are  first  covalently  linked  to  the  microbead  surface  in  sets,  with  each  set  pre¬ 
pared  separately  and  displaying  a  single  probe  on  its  surface.  (Similar  to  flow  cytometric 
assays  platforms128,149-156.)  A  microbead  suspension  is  then  streamed  through  the  microflu- 
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idic  flow  channel,  and  patterning  techniques  are  used  to  position  and  fix  the  microbeads  in 
separable  positions  along  the  channel  to  form  a  probe  array.  A  library  of  multiple  probes 
is  constructed  by  sequentially  streaming  the  microbeads  in  sets  with  each  set  displays  only 
one  probe.  After  a  set  is  introduced,  the  locations  of  the  members  of  the  set  in  the  channel 
are  noted  before  the  next  set  is  introduced  so  a  registry  of  the  probe  on  each  arrayed  mi- 
crobead  can  be  compiled.  Fluorescently  labelled  target  can  then  be  streamed  through  the 
cell  and  over  the  microbeads  to  undertake  a  screening  assay,  with  binding  events  recognized 
by  identification  of  fluorescing  microbeads,  and  the  binding  probe  identified  by  referencing 
the  library  registry. 

The  patterning  used  to  fix  the  microbeads  in  array  positions  in  the  microfluidic  cell  is  the 
key  to  the  assembly  process.  The  most  straightforward  method  for  fixing  beads  in  position 
on  a  surface  would  be  to  pattern  the  surface  in  an  array  of  adhesive  domains,  i.e.  islands 
with  functional  groups  which  stick  to  the  microbeads  when  the  beads  are  deposited  on  the 
surface.  Surfaces  patterned  with  functionalities  which  allow  beads  to  stick  by  electrostatic 
forces  or  covalent  bonding  have  been  used  to  assemble  a  microbead  array  on  a  substrate 
outside  of  a  flow  cell 157-161 ,  and  these  arrays  have  been  utilized  to  screen  biomolecules  by 
incubating  static  analyte  solutions  above  the  array  in  analogy  with  the  a  flat  microarray. 
Similarly,  a  microbead  array  can  be  assembled  within  a  microfluidic  flow  channel  by  stream¬ 
ing  a  suspension  of  beads  through  though  the  channel  and  allowing  them  to  stick  to  adhesive 
domains  on  the  floor  of  the  channel.  However,  these  immobilized  microbeads  are  typically 
displaced  by  the  target  flow  because  the  binding  of  the  microbead  to  the  surface  is  not  over 
enough  bead  area  to  insure  a  strong  retaining  force.  For  this  reason,  studies  have  focused 
on  using  retaining  structures  to  position  the  microbeads  in  fixed  locations.  The  two  designs 
that  are  the  most  common  are  wells  patterned  into  the  bottom  surface  of  the  microfluidic 
cell162,163  and  traps  arranged  as  a  microfluidic  obstacle  course  (DiCar lo  et  al 164,165 ,  Cooper 
et  al166-171,  Huebner  et  al. 172  and  Nehorai  et  al173,174.  Well  designs  size  the  cavity  to  be 
only  slightly  larger  than  the  microbead  to  be  retained,  so  that  only  one  microbead  is  cap¬ 
tured  at  a  location.  Well  deposition  has  also  been  enhanced  by  using  electric  and  magnetic 
fields  to  assist  in  the  capture175-177,  or  by  using  holes  placed  in  the  well  and  connected  to 
a  drain  to  provide  fluid  suction  (see  McDevitt  et  al 178-187  and  Ketterson188).  Trap  obstacle 
courses  in  microfluidic  cells  have  been  used  predominantly  to  hydrodynamically  array  cells, 
or  droplets  containing  cells,  but  a  few  reports  have  used  the  trapping  geometry  to  array  mi- 
crobeads,  including  ones  which  has  studied  optimal  geometries  for  the  trap  shape  and  course 
arrangement  to  insure  high  capture  efficiency173,174,  and  our  effort  which  demonstrated  the 
arraying  of  microbeads  with  a  lipid  bilayer  (lipobeads)  for  screening  membrane  receptors189. 
Traps  designed  to  capture  a  single  microbead  are  usually  half-open,  “V”  shaped  retaining 
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structures  which  span  the  height  of  the  channel  and  are  oriented  with  the  open  part  of  the 
cavity  facing  the  flow.  These  designs  typically  have  an  opening  at  the  back  end  of  the  trap, 
to  facilitate  flow  through  the  trap  enabling  easy  capture  and  allowing  target  to  stream  over 
a  trapped  microbead  (see  Fig.  5.1  used  in  our  study189  and  similarly  in173,174). 

The  major  question  which  emerges  when  retaining  structures  are  used  to  assemble  probe- 
functionalized  microbeads  in  a  microfluidic  cell  to  construct  a  probe  library  is  the  effect  of 
the  retaining  structure  on  the  mass  transfer  of  the  target  to  the  probes  on  the  microbead 
surface.  As  a  reference,  consider  first  the  target  transport  in  the  idealized  case  in  which  the 
microbeads  are  (fictitiously)  attached  to  the  bottom  of  the  flow  channel  without  a  retaining 
structure.  The  target  transport  occurs  in  two  steps.  Target  diffuses  from  the  stream  passing 
over  the  microbeads  to  the  bead  surface  by  crossing  the  flow  streamlines,  and  then  kinetically 
conjugates  to  the  probes  on  the  bead  surface.  As  a  result,  a  diffusion  zone  is  created  around 
the  microbead  in  which  the  concentration  falls  below  the  inlet  concentration.  Within  this 
zone,  a  target  concentration  gradient  is  formed  normal  to  the  surface  with  the  concentration 
decreasing  from  the  inlet  concentration  at  the  top  of  the  zone  to  a  reduced  concentration  in 
the  sublayer  of  liquid  immediately  adjacent  to  the  surface  (the  sublayer  concentration).  As 
binding  proceeds  and  the  surface  becomes  saturated  with  target,  the  kinetic  conjugation  rate 
decreases,  causing  the  sublayer  concentration  to  increase  as  the  gradients  in  concentration 
in  the  diffusion  zone  relax  and  eventually  disappear  at  equilibrium  in  which  case  the  concen¬ 
tration  around  the  bead  becomes  equal  to  the  inlet  concentration.  Hence  the  overall  binding 
rate  of  the  target  to  the  probe  is  determined  by  an  in-series  transport  process  of  diffusion 
and  surface  kinetics.  When  the  convective  flow  along  the  surface  is  fast  relative  to  the  diffu¬ 
sion  normal  to  the  surface,  this  transport  picture  is  retained,  but  the  diffusion  zone  becomes 
a  thin  boundary  layer  and  the  diffusive  flux  to  the  surface  becomes  large,  and  the  surface 
quickly  saturates.  For  any  flow  regime,  it  is  the  kinetic  conjugation  step  of  the  target  to  the 
surface  probe  that  actually  populates  the  surface  with  the  target.  This  step  is  described  by 
a  surface  kinetic  equation  that  relates  the  rate  of  conjugation  to  the  instantaneous  surface 
concentration  of  the  target,  and  the  target’s  sublayer  concentration  adjoining  the  surface, 
as  for  example,  the  Langmuir  equation.  If  this  surface  kinetic  rate  is  much  slower  than  the 
diffusion  rate  (due  to  the  concentration  gradient  normal  to  the  surface  in  the  diffusion  or 
boundary  layer  zone),  the  concentration  gradient  decreases  and  the  concentration  of  target 
in  the  zone  approaches  the  inlet  concentration.  In  the  limit  of  a  very  slow  kinetic  conjugation 
relative  to  diffusion,  the  overall  binding  rate  is  simple  to  describe  since  it  is  only  determined 
by  integrating  the  surface  kinetic  equation  with  the  sublayer  concentration  equal  to  the  inlet 
concentration  (the  ideal  kinetic  limit).  However,  if  the  kinetic  rate  is  of  order  of  the  rate  of 
bulk  diffusion,  then  the  diffusion  gradients  reduce  the  overall  binding  rate  from  the  simple 
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Figure  5.1:  Microbeads  hosting  surface  probes  are  arrayed  in  a  microfluidic  cell  using  a  trapping  course. 

kinetic  limit  and  the  transport  process  is  mixed  with  the  sublayer  concentration  obtained 
from  the  solution  of  the  convective  diffusion  mass  transfer  equations. 

When  retaining  structures  such  as  traps  and  wells  are  used  to  array  microbeads  in  fixed 
positions,  the  confinement  reduces  the  convective  flow  around  the  microbead,  and  decreases 
the  mass  transfer  rate  of  target  to  the  surface  probes  by  increasing  the  size  of  the  diffusion 
zone  or  the  boundary  layer  thickness.  A  reduction  in  the  target  transport  rate  affects  the 
performance  of  the  screening  assay,  as  it  increases  the  time  necessary  to  achieve  the  par¬ 
ticular  level  of  target  concentration  necessary  to  detect  a  signal.  The  focus  of  this  study  is 
to  examine  this  effect  for  the  case  in  which  microbeads  are  arrayed  by  an  obstacle  course 
of  traps  (see  Fig.  5.1).  While  the  analysis  of  the  binding  of  targets  to  surface  probes  as  a 
screening  platform  has  been  studied  in  detail  for  the  case  of  probe  “patches”  situated  on  a 
microchannel  wall133,190-199,  the  binding  of  target  to  microbeads  captured  in  traps  has  re¬ 
ceived  very  limited  attention.  Bau  et  al200  studied  the  geometry  of  microbeads  sandwiched 
between  the  top  surface  of  a  flow  channel  and  a  shallow  well  at  the  bottom  of  the  channel 
(the  microbeads  were  preassembled  in  the  shallow  wells  before  closing  the  cell) ,  and  obtained 
solutions  for  the  target  concentration  on  the  microbead  surface  as  a  function  of  the  stream 
velocity  and  target-probe  kinetic  rate  constants.  These  simulations  represent  approximately 
the  idealized  case  of  unobstructed  microbeads.  Our  aim  is  to  construct  numerical  solutions 
for  the  hydrodynamic  flow  and  target  mass  balance  equations  to  determine  the  surface  con¬ 
centration  on  the  bead  surface  as  a  function  of  time  when  the  microbeads  are  entrapped.  The 
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Figure  5.2:  (a)  Schematic  of  the  trapping  array  used  in  the  experiments  and  modeled  in  the  simulation, 
and  an  inset  detail  of  an  individual  trap  with  a  microbead.  The  analyte  solution  is  streamed  from  bottom 
to  top.  (b)  Schematic  of  the  unit  cell  computational  domain  with  dimensions. 

solutions  are  compared  to  the  binding  rates  of  target  to  probes  on  unconfined  microbeads 
in  a  channel  (as  the  Bau  simulations)  to  determine  the  effect  of  the  trap.  The  simulations 
will  also  be  used  to  identify  regimes  of  slow  surface  kinetics  for  which  the  overall  binding 
rate  approaches  the  ideal  kinetic  limit.  This  criteria  is  of  interest  in  experimental  studies  of 
target-probe  binding  kinetics  since  the  ideal  kinetic  binding  rate  is  a  simple  expression  easily 
compared  to  experimental  measurements  of  the  binding  rate  to  determine  the  appropriate 
kinetic  model  or  to  evaluate  surface  kinetic  rate  constants.  We  will  also  undertake  screening 
experiments  on  the  microbead  array  to  illustrate  how  the  simulations  can  be  used  to  select 
experimental  conditions  on  the  imposed  flow  rate  so  that  the  experimentally  measured  over¬ 
all  binding  rate  of  the  target  to  the  surface  probes  is  controlled  by  the  surface  conjugation 
step  so  the  overall  binding  rate  can  be  easily  simulated.  For  this  demonstration,  we  will  use 
the  protein  NeutrAvidin  as  the  target,  and  biotin  (its  binding  partner)  as  the  probe.  Biotin 
will  be  bound  to  the  surface  of  the  microbeads,  and  the  NeutrAvidin  will  be  fluorescently 
labelled  so  that  its  accumulation  on  the  microbead  surface  can  be  quantitatively  measured 
to  obtain  the  overall  binding  rate. 
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5.2  Hydrodynamic  and  Mass  Transfer  Numerical  Sim¬ 
ulations 

5.2.1  Formulation 

The  simulation  of  the  target  analyte  flow,  and  the  convective  diffusion  and  binding  of  the 
target  to  the  probes  on  the  microbeads  arrayed  in  traps  in  the  flow  cell  is  constructed 
as  a  three  dimensional  simulation  which  models  directly  the  trapping  geometry  used  in 
our  earlier  publication  on  lipobead  arraying189  and  the  Neutravidin-biotin  binding  assay 
experiments  which  are  described  later  on.  A  schematic  of  the  experimental  layout  of  the 
trap  arrangement  in  the  obstacle  course,  and  the  detail  of  an  individual  trap,  is  shown  in 
Fig.  5.2a.  The  traps,  in  a  wide  channel,  are  arranged  in  parallel  rows  (perpendicular  to  the 
flow)  which  are  separated  by  200  //m,  and  columns  (along  the  flow)  which  are  separated  by 
260  /im  and  staggered  symmetrically  with  an  offset  of  130  /tin  (Fig.  5.2b).  The  channel 
is  of  rectangular  cross  section,  with  an  overall  channel  width  of  3  mm,  a  length  of  6.1  mm 
and  a  height,  h  is  60  /im.  A  total  of  144  traps  are  arranged  in  a  cluster  at  approximately 
the  center  of  the  flow  channel.  (The  open  area  downstream  of  the  cluster  accomodates  an 
entry  port  for  the  introduction  of  the  microbeads,  and  both  ends  are  connected  to  channels 
of  smaller  width(300  //in  from  which  the  analyte  flow  enters  and  exits.)  The  “V”  shaped 
configuration  of  the  trap  contains  an  aperture  at  the  back  end  which  is  20  /im  in  width, 
permitting  the  microbeads  which  have  a  diameter  2 a  equal  to  42.3  /an  to  be  captured.  The 
trap  height  is  equal  to  the  channel  height.  The  microbeads,  which  are  made  of  glass  in 
the  experiments,  settle  to  the  bottom  of  the  flow  cell  during  the  entrapment  process,  and 
therefore  an  approximately  20  /im  gap  exists  between  the  microbead  and  the  top  wall  of  the 
microchannel,  permitting  flow  over  a  microbead  when  it  is  localized  in  the  trap.  To  simplify 
the  calculations,  we  ignore  the  edge  effects  of  the  side  walls  of  the  channel,  and  undertake 
calculations  on  a  unit  cell  (Fig.  5.1b)  with  symmetry  boundary  conditions  for  the  flow  and 
mass  transfer  along  the  side  walls  of  the  unit  cell.  A  cartesian  coordinate  system  is  located 
with  an  origin  at  at  the  center  of  the  entrance  and  the  bottom  wall  of  the  channel,  with  y 
in  the  flow  direction,  x  perpendicular  to  the  side  walls  and  z  perpendicular  to  the  bottom 
wall.  The  unit  cell  consists  of  4  of  the  twelve  rows  of  the  trapping  course  bounded  by  inlet 
and  exit  cross-section  planes  at  y  =  0  and  L,  respectively.  (Calculations  demonstrate  that 
the  addition  of  a  further  row  did  not  change  either  the  flow  or  the  mass  transfer.)  The 
microbeads  in  the  unit  cell  are  symmetrically  placed  at  the  centerline  of  the  traps,  and  are 
fixed  in  a  position  touching  at  one  point  the  bottom  of  the  flow  cell.  While  in  experiments 
the  analyte  flow  pushes  the  microbead  directly  against  the  trap  walls,  in  the  simulations,  the 
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microbeads  are  separated  by  a  distance  of  1  p  to  avoid  complexities  in  the  mesh  construction 
for  the  numerical  solution.  Reducing  this  separation  distance  did  not  change  the  results.  A 
uniform  velocity  U  (in  the  y  direction)  and  a  uniform  target  concentration  cQ  is  imposed  at 
the  upstream  inlet  (y  =  0);  at  the  exit  (y  =  L)  a  zero  pressure  condition  is  imposed,  and 
the  derivative  of  the  concentration  in  the  y  direction  is  set  equal  to  zero.  The  symmetry 
boundary  conditions  along  the  side  walls  (x  =  ±.W)  of  the  unit  cell  are  that  the  velocity 
in  the  x  direction  is  equal  to  zero,  and  the  derivative,  in  the  y  direction  of  the  velocity  in 
the  y  and  z  direction  and  the  concentration  is  equal  to  zero.  At  the  initial  time,  the  target 
concentration  in  the  unit  cell  is  assumed  to  be  uniform  and  equal  to  cQ. 

The  target  solution  is  assumed  to  be  an  incompressible  Newtonian  fluid  with  density  p  and 
viscosity  p  of  water  (p=103  kg  m-3  and  p  =  10-3  kg  m_1sec_1)  independent  of  the  analyte 
concentration  which  is  assumed  to  be  dilute.  The  steady  hydrodynamic  flow  is  described  by 
the  solution  of  the  continuity  and  Navier-Stokes  equations,  written  in  nondimensional  form 
in  terms  of  the  velocity  v  scaled  by  the  velocity  U. 

V  •  v  =  0  (5.1) 

9te  [v  •  Vv]  =  — Vp  +  V2v  (5.2) 

where  the  nondimensional  variables  V  and  V2  are  the  gradient  and  Laplacian  operators 
(scaled  by  h),  p  is  pressure  (nondimensionalized  by  pU/h ),  and  91c  =  is  the  Reynolds 

p 

number.  The  continuity  and  Navier-Stokes  equations  are  solved  in  cartesian  coordinates 
with  the  inlet  and  outlet  conditions  detailed  above,  and  boundary  conditions  of  no  slip  on 
the  interior  walls  of  the  trap  and  the  bead  surface  and  symmetry  conditions  along  the  unit 
cell  border  x  =  ±W.  The  solution  is  obtained  numerically  using  finite  elements  implemented 
with  the  COMSOL  Multiphysics  simulation  package  (4.2),  using  both  triangular  and  quadri¬ 
lateral  meshes. 

The  concentration  field  of  the  target  is  obtained  by  solving  the  convective  diffusion  equa¬ 
tion 

dc 

— — |-  Pew-  Vc  =  V2c  (5.3) 

or 

In  the  above,  c  is  the  concentration  of  target  (nondimensionalized  by  the  inlet  concentration 
cD),  v  is  the  steady  velocity  obtained  above,  and  Pe  is  the  Peclet  number  defined  as  Pe  = 
Ua/T>,  with  T>  the  diffusion  coefficient  of  the  target  and  time  is  scaled  by  the  diffusion  time 
r©  =  a2/T>  the  time  required  for  target  to  diffuse  over  the  length  scale  a  of  the  bead.  Pe  is  the 
ratio  of  the  time  required  for  fluid  to  convect  over  the  microbead  ( a/U )  to  the  diffusive  time 
scale  r©.  Eq.  5.3  is  solved  with  the  inlet  and  outlet  conditions  detailed  above,  and  assuming 
zero  flux  of  target  on  the  surfaces  of  the  trap  and  channel  walls  and  symmetry  conditions  on 
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the  borders  (x  =  ±VU)  of  the  unit  cell.  The  boundary  condition  on  the  microbead  surface 
models  the  binding  kinetics,  and  is  described  in  detail  in  the  following. 

The  kinetic  binding  of  target  from  the  sublayer  of  liquid  immediately  adjacent  to  the 
microbead  surface  to  probes  on  the  microbead  surface  is  a  bi-molecular  process,  which  is 
described  generally  by, 

afc}=w_i4fc}  <5'4) 

where  T  is  the  surface  concentration  of  target,  Too  is  the  maximum  concentration  of  target 
which  can  bind,  t  is  time,  cs  is  the  sublayer  concentration  of  target  at  the  surface  and  ka  and 
kd  are  the  association  and  disassociation  rate  constants,  respectively  and  <p  is  the  fraction  of 
sites  available  when  the  target  concentration  is  T.  To  calculate  cf>,  we  specify  first  the  density 
of  probes  fixed  on  the  surface  as  Tp,  and  the  projected  area  per  molecule  of  the  target  as  At. 
When  the  surface  density  of  probes  (Tp)  is  small  enough  such  that  bound  target  does  not 
overlap  unbound  sites  (TpAt  <  1),  the  maximum  concentration  of  bound  target  is  equal  to 
the  probe  density  Too  =  Tp  and  the  fraction  of  sites  available  is  b  =  1  —  j— ,  leading  to  a  simple 
Langmuir  kinetic  scheme.  Under  many  circumstances,  the  target’s  projected  area  is  much 
larger  than  the  probe  area,  and  surface  densities  can  be  large  enough  so  that  TpAt  >  1.  In 
this  case,  as  binding  proceeds,  the  binding  of  targets  to  probes  becomes  sterically  hindered, 
and  the  fraction  of  sites  available  is  not  given  by  0  =  1  —  Models  for  0  to  account  for 
steric  hindrance  based  on  Monte  Carlo  simulations  of  random  sequential  adsorption  have 
been  developed  (e.g.201),  and  these  also  calculate  AtT^  as  a  function  of  TpAt.  Here,  for 
our  simulations,  we  generalize  the  Langmuir  equation  and  assume  that  b  =  1  —  tI—  where 
the  maximum  packing  can  be  obtained  from  the  random  sequential  model  calculations.  In 
the  experiments,  the  surface  density  is  low  enough  such  that  TpAt  <C  1  so  the  Langmuir 

r  k 

model  is  applicable  with  T ^  =  Tp.  The  equilibrium  surface  density  (reg)  is  — —  =  - - 

I  oo  1  T  k 
k  c 

where  k  =  -  °.  Many  biomolecular  binding  processes  are  nearly  irreversible  (i.e.  k  W  1), 

kd 

and  in  our  screening  experiments  on  the  binding  of  Neutravidin  to  surface  biotin  previous 
experiments  show  that  this  is  certainly  the  case.  Therefore  in  our  simulations,  we  will 
neglect  the  disassociation  step.  The  boundary  conditions  on  the  microbead  surface  equate 
the  diffusive  of  target  to  the  kinetic  rate  of  binding,  and  sets  the  rate  of  accumulation  of 
target  at  the  microbead  surface  to  the  flux: 


{n  •  Vc}bead  =  Z3a  [c,  { 1  —  f } 

(5.5) 

^  =  {n'Vc}bK>J 

(5.6) 
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Figure  5.3:  Simulation  of  the  flow  around  the  upstream  center  microbead  in  the  unit  cell  (microbead  “A” 
in  Fig.  5.2)  in  a  trapped  (a)  and  untrapped  (b)  configuration  for  a  flow  Reynolds  number  9tc  =  0.2.  The 
magnitude  of  the  velocity  (divided  by  the  average  velocity  U)  is  plotted  in  the  x  —  z  plane  corresponding  to 
the  equatorial  section  of  the  bead  perpendicular  to  the  flow,  y  =  Y 


where  n  is  the  outward  normal  to  the  microbead  surface,  cs  is  the  nondimensional  sublayer 

r  k  F  a 

concentration,  Fl  =  — Da  is  the  Damkohler  number,  Da  =  -  °°  and  F  is  the  surface 
cQa  D 

concentration  scaled  by  the  maximum  binding  density  Too.  The  parameter  Q  is  the  ratio  of 
the  adsorption  depth  the  distance  above  the  surface  which  contains  (per  unit  area) 

enough  target  to  saturate  the  surface  to  the  microbead  radius  a.  The  Damkohler  number 
is  the  ratio  of  the  characteristic  kinetic  flux  to  the  surface  {kac0F ^  to  the  characteristic 
diffusive  flux  on  the  particle  length  scale  (Dc0/a).  The  surface  concentration  (T)  is  a  function 
of  the  position  on  the  bead  surface,  and  we  denote  by  f  the  average  value  on  the  bead  surface. 
F  —  j  JA  FdA  where  A  is  the  area  of  the  microbead  surface. 

The  mass  transfer  equations  are  integrated  using  the  COMSOL  simulation  package,  using 
forward  marching  in  time.  An  adaptive  time  step  is  used  in  which  for  the  earliest  times  for 
large  Pe  the  time  step  is  set  to  one  order  of  magnitude  smaller  then  the  time  required  for 
formation  of  the  boundary  layer  around  the  microbead,  S2/T>  and  for  order  one  Pe  the 
diffusion  time  scale  a2 /T> 
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5.2.2  Hydrodynamic  Flow  Simulations 

In  a  typical  microfluidic  trapping  geometry,  channel  heights  and  widths  are  of  order  102 
jim  and  103  //in,  respectively  and  for  flow  rates  of  order  10-1  -  102  /A/ min,  the  average 
velocity  U  ~  10  —  102  pm/s.  Hence  IHe,  based  on  the  microbead  radius  (typically  10  //in)  is 
usually  less  than  10_1,  and  the  motion  is  predominately  a  Stokes  inertialess  flow.  In  these 
simulations,  we  choose  Ate  =  0.2,  which  is  the  experimental  value.  The  hydrodynamic  flow 
pattern  around  a  microbead  situated  at  the  bottom  of  the  downstream  trap  at  the  center  of 
the  unit  cell  (“A”  in  Fig.  5.2  with  center  at  position  y  =  Ya)  is  shown  in  Fig.  5.3a.  The  flow 
pattern  is  the  magnitude  of  the  velocity  (divided  by  the  average  velocity  U )  in  the  x  —  z  plane 
(  y  =  Ya)  perpendicular  to  the  flow  and  through  the  bead  equator.  (The  flow  is  symmetric 
in  x  since  x  =  0  is  the  the  symmetry  plane  of  the  unit  cell.  The  projection  of  the  figure  is 
an  oblique  view,  accounting  for  the  elliptical  appearance  of  the  spherical  bead.)  Due  to  the 
staggered  configuration  of  the  traps,  most  of  the  flow  circulates  around,  rather  than  through 
the  trap,  as  the  hydrodynamic  resistance  through  an  occupied  trap  is  much  larger  than  the 
resistance  in  the  open  space  between  the  traps.  Within  the  trap,  the  figure  makes  clear  that 
the  maximum  velocity,  as  would  be  expected,  is  in  the  large  gap  (~  20  /irn)  between  the  top 
of  the  channel  and  the  bead  [h  <  z  <  2a,  y  =  Ya,  x  =  0).  Analyte  flows  through  this  space, 
and  then  out  the  back  of  the  trap.  The  lateral  gap  thickness  between  the  trap  sidewalls  and 
the  microbead  ( Lmin  —  a  >  x  —  a  >  0,  y  =  Ya,  z  =  a)  in  the  equatorial  plane  of  Fig.  5.3a,  is 
much  smaller  (Lmin  =  10  //in)  and  the  velocity  in  this  gap  is  less  than  between  the  bead  and 
the  top  wall  of  the  channel.  The  flow  velocity  in  the  same  plane  around  microbead  “A”  in 
the  untrapped  configuration  (the  same  array  of  microbeads  as  in  Fig.  5.2  but  without  traps) 
is  given  in  Fig.  5.3b,  and  as  expected,  the  flow  around  the  lateral  sides  of  the  microbead  is 
larger  due  to  the  absence  of  the  trap.  In  this  case  the  analyte  stream  freely  flows  over  most 
of  the  microbead,  except  the  lower  half  of  the  bottom  hemisphere  which  is  obstructed  by 
contact  with  the  bottom  wall.  As  fluid  is  no  longer  directed  around  the  blunt  edges  of  the 
traps,  the  fluid  velocity  in  the  space  occupied  between  the  symmetry  plane  and  the  position 
formelry  occupied  by  the  outer  trap  wall  is  reduced  relative  to  when  the  trap  is  present.  In 
fact  in  this  space,  spanning  a  distance  of  ~  2 a  from  the  microbead  to  the  symmetry  plane, 
approaches  two  dimensional  Poiseuille  flow  with  the  maximum  velocity  at  the  midplane  of 
the  channel  in  this  space  equal  to  3/2  of  the  upstream  average  velocity  U  in  front  of  the  trap 
cluster. 

As  stated  earlier,  in  the  regime  in  which  the  convective  rate  of  transport  of  target  along 
the  microbead  is  large  relative  to  the  diffusive  flux  to  the  surface,  the  mass  transfer  is  through 
a  thin  boundary  layer  around  the  particle.  The  thickness  of  this  boundary  layer  decreases 
with  the  local  velocity  gradient  (or  shear  rate)  at  the  microbead  surface,  and  the  target 
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Figure  5.4:  Simulation  of  the  flow  velocity  in  the  y  (streamwise)  direction  (divided  by  U)  between  (a)  the 
top  of  the  channel  and  the  microbead  (0  <  z  —  2a  <  h  —  2a,  y  =  Y a  ,  x  =  0)  and  (b)  the  side  of  the  microbead 
and  the  channel  wall  (0  <  x  —  a  <  Lmm,  y  =  Ya,  z  =  a)  and  the  trap  wall  and  the  symmetry  plane  of  the 
unit  cell  for  the  upstream  center  microbead  in  the  unit  cell  (microbead  “A”  in  Fig.  5.2)  and  comparison  to 
the  untrapped  microbead  for  91  c  =  0.2.  The  inset  shows  the  distribution  in  the  gap  between  the  microbead 
and  the  trap  wall. 


flux  becomes  correspondingly  larger.  Hence  an  examination  of  the  surface  shear  rate  is 
important  in  understanding  the  target  transport  in  the  convective  regime.  Detailed  profiles 
of  the  velocity  in  the  y  (stream)  direction,  as  a  function  of  the  distance  across  the  gaps  at 
the  top  and  sides  of  the  trapped  microbead,  are  shown  in  Fig.  5.4(a)  and  compared  to  the 
profiles  for  a  microbead  without  a  trap  (Fig.  5.4(b).  For  the  trapped  microbead,  the  flow  in 
the  lateral  gap  is  very  small  (velocities  are  of  order  10_2U),  and  the  shear  rate  at  the  wall 
is  also  small.  For  the  top  gap,  the  velocity  (of  order  U )  and  shear  rate  are  much  larger,  and 
the  target  mass  transfer  to  the  upper  hemisphere  of  the  microbead  should  be  expected  to  be 
larger  than  the  lower  hemisphere  (see  the  following  section).  In  contrast,  in  the  case  of  the 
untrapped  microbead,  while  the  shear  rate  at  the  top  is  comparable  to  that  of  the  trapped 
microbead,  the  shear  rate  along  the  lateral  side  is  much  larger  as  the  trap  no  longer  blocks 
the  flow,  and  the  mass  transfer  rate  should  be  correspondingly  more  symmetrical  between 
the  top  and  bottom  hemispheres.  Interestingly,  the  flow  rate  through  the  top  gap  is  less  for 
the  untrapped  microbead  than  the  trapped  one,  as  fluid  streams  more  easily  laterally  around 
the  microbead. 
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Figure  5.5:  The  binding  of  a  target  to  the  surface  of  a  microbead  in  a  trap  (microbead  “A”,  y  =  Ya)  for 
Pe  =  3.5  (Q  =  32,  94e  =  0.2):  (a)  The  average  nondimensional  surface  concentration  on  the  surface  of  the 
microbead,  T/Too,  as  a  function  of  r  for  Da  =1,  10  and  102  alongside  the  kinetic  limits  for  these  values  of 
Da.  (b)  Target  concentration  diffusion  zones  for  Da=  10  in  the  equatorial  plane  z  =  a  and  the  projection  of 
the  surface  concentration  averaged  for  the  upper  (z  >  a)  and  lower  ( z  <  a)  hemispheres  of  the  microbead 
for  r  =  .136,  20.4  and  30.26  and  Da  =  10,  Q  =  32,  SRe  =  2xl0“4. 


5.2.3  Mass  Transfer  Simulations 


The  objective  of  the  mass  transfer  calculations  is  to  determine  the  nondimensional  aver¬ 
age  concentration  on  the  microbead  surface  (f  /Too)  as  a  function  of  nondimensional  time 
t.  Aside  from  the  geometric  parameters  of  the  unit  cell,  the  average  concentration  is  a 

f 

function  of  Cl,  Da  and  Pe,  and  therefore  — — (r;  Da.  Pe,  Cl).  In  most  screening  applications 

t  OO 

the  concentration  of  the  target  is  low  enough  or  the  binding  capacity  large  enough  so  that 
Too/co  o,  and  hence  Cl  ;§>  1.  As  a  result,  the  accumulation  of  target  on  the  surface,  as  is 
clear  from  eq.  5.6,  occurs  in  accordance  with  a  slower  time  scale  r*  =  r /Cl,  as  the  larger  the 
value  of  O,  the  slower  is  the  binding  rate.  Rescaling  eq.  5.6  with  r*  reformulates  the  balance 

dr 

to  become  independent  of  Cl,  i.e.— —  =  (n  •  Vc}b  d.  The  second  boundary  condition,  eq. 
5.5,  is  independent  of  Cl  (as  is  the  zero  flux  boundary  condition  on  the  trap  and  channel 


walls),  and  the  convective  diffusion  equation  (eq.  5.3)  becomes 


1  dc 


+  Pe  v  •  Vc  =  V2c. 


Cl  dr* 

Solutions  of  the  rescaled  equations  are  therefore  independent  of  Cl  to  leading  order  for  Cl  1, 
as  the  transport  dynamics  in  the  bulk  becomes  quasi-stationary  and  dictated  in  time  only  by 
the  change  in  flux  on  the  surface  in  r*.  We  find  that  the  time  derivative  term  in  the  bulk  is 
only  important  in  the  earliest  times  in  which  the  diffusion  zone  around  the  particle  develops 
from  a  uniform  concentration.  Otherwise  the  bulk  time  derivative  term  is  negligible,  and  as 
such,  simulations  integrated  in  r  done  for  a  particular  value  of  Da  and  Pe  and  f2  (assumed 
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much  larger  than  one)  can  be  used  for  another  large  value  of  Q  by  just  rescaling  time  in  the 
completed  solution  by  the  ratios  of  the  two  values  of  D.  Therefore  for  Q  1,  simulations 
varying  fl  are  not  necessary.  In  the  general  set  of  simulations  presented  here  Q  is  fixed  at  32 
which  is  smaller  than  the  experimental  value,  but  is  used  to  accelerate  the  integration  time. 

The  target  surface  concentration  on  the  microbead  as  a  function  of  time  (T(r))  is  de¬ 
termined  by  a  balance  between  the  Peclet  (Pe  =  Ua/T> )  and  Damkohler  ( Da  =  kaT^a/D) 
numbers.  Pe  is  the  ratio  of  the  time  scale  for  diffusion  over  a  distance  a  normal  to  the  sur¬ 
face  (  a2 /T>)  to  the  nominal  scale  for  fluid  to  convect  over  the  microbead  surface  (estimated 
as  a/U ).  Target  proteins  or  smaller  biomolecular  ligands  have  molecular  weights  of  order 
103  — 104  and  corresponding  diffusion  coefficients  of  ~  102/rm2/s;  for  U  ~  102  —  104/im/s  and 
a  of  order  10  /im,  Pe  is  of  order  1  —  103,  which  encompasses  different  transport  behavior  in 
the  analyte  stream.  For  Pe  of  order  1-10,  diffusion  zones  (regions  where  the  concentration 
is  less  than  the  inlet  concentration  c0)  develop  normal  to  the  microbead  surface  and  are  of 
thickness  of  order  a.  However,  for  increasingly  larger  Pe  (101  -  103)  the  nominal  convection 
time  is  much  faster  and  the  diffusion  zone  takes  the  form  of  asymtotically  thin  boundary 
layers.  In  the  simulations  we  consider  an  order  one  value  (Pe  =  3.5)  with  order  a  diffu¬ 
sion  zones,  and  Pe  =  3500,  characteristic  of  thin  diffusion  zones  and  representative  of  our 
experiments.  Values  for  the  Damkohler  number  Da,  the  ratio  of  the  characteristic  kinetic 
flux  to  the  surface  ( kac0T ,*,)  to  the  characteristic  diffusive  flux  on  the  particle  length  scale 
(D cQ/a),  are  dependent  on  the  surface  kinetics.  For  Da  of  order  1  or  smaller,  the  kinetic 
rate  step  is  of  the  same  order  or  slower  than  diffusion,  and  the  gradients  in  concentration  in 
the  diffusion  zone  are  small  as  the  zone  is  nearly  at  the  bulk  concentration  ca.  In  the  limit 
Da  — >■  0,  the  bulk  concentration  around  the  microbeads  becomes  equal  to  cQ,  and  the  trans¬ 
port  becomes  kinetically  limited  with  the  surface  concentration  independent  of  the  position 
on  the  microbead  surface  and  given  by  the  solution  to  eqs.  5.5  -  5.6  after  elimination  of  the 
diffusive  flux  term. 


=  1  —  e" 


(5.7) 


For  larger  Da,  the  kinetic  step  becomes  faster  and  the  concentration  gradient  is  large  in 
the  diffusion  zone.  In  the  simulations,  we  study  values  of  Da  =  1,  characteristic  of  the  ex¬ 
periments,  and  two  larger  values  of  10  and  100  representing  large  gradients  in  the  diffusion 
zones. 

Consider  first  a  low  value  of  Pe  =  3.5.  The  average  nondimensional  target  surface  con¬ 
centration  f  (scaled  by  Too)  on  microbead  “A”  in  a  trap  as  a  function  of  nondimensional  time 
for  Pe  =  3.5  and  Da  =  1,  10  and  102  (Q  =32)  is  given  in  Fig.  5.5(a)and  compared  to  the 


87 


Figure  5.6:  The  binding  of  a  target  to  the  surface  of  a  microbead  in  a  trap  (microbead  “A”,  y  =  f'i) 
for  Pe  =  3500  (Q  =  32,  Die  =  0.2):  (a)  The  average  nondimensional  surface  concentration  on  the  surface 
of  the  microbead,  T/Too,  as  a  function  of  r  for  Da  =1,  10  and  102  alongside  the  kinetic  limits  for  these 
values  of  Da.  (b)  Target  concentration  boundary  layer  for  Da= 10  in  the  equatorial  plane  z  =  a  and  the 
projection  of  the  surface  concentration  averaged  between  the  upper  (z  >  a)  and  lower  (z  <  a)  hemispheres 
of  the  microbead  for  r  =  .0.012,  2.9  and  5.4  and  Da  =  10. 


kinetic  limits,  eq.  5.7,  for  each  of  the  values  of  Da.  Due  to  the  fact  that  the  convective  and 
diffusive  fluxes  are  of  the  same  order,  large  diffusion  zones  develop  around  the  microbead, 
and  this  slows  the  transport  rate  of  the  target  to  the  surface.  As  a  result,  for  every  value  of 
Da  in  the  figure,  the  binding  rate  is  much  slower  than  the  corrsponding  kinetic  limit  which 
assumes  the  concentration  is  uniform  and  equal  to  c0  around  the  microbead.  The  diffusion 
zones  around  the  microbead  are  depicted  for  Da  =10  in  Fig.  5.5(b)  for  three  nondimensional 
times,  r  =  .136,  20.4  and  30.26  in  the  equatorial  plane  z  —  a.  Diffusion  brings  target  to 
the  surface  where  it  kinetically  binds  to  the  receptors,  and  the  bulk  concentration  decreases 
in  the  zone  in  the  direction  normal  to  the  microbead  surface.  As  time  progresses,  target 
depletion  becomes  more  pronounced  and  the  gradients  become  sharper.  This  is  especially 
evident  at  the  downstream  end  of  the  microbead,  where  the  concentration  drops  to  zero 
through  the  opening  at  the  back  of  the  trap  by  r  =  30.26.  At  this  time  the  region  in  back 
of  the  trap  is  also  depleted.  This  significant  depletion  is  due  to  the  reduced  convective  flow 
at  the  back  of  the  microbead,  relative  to  the  front,  and  as  a  result  the  back  side  of  the  mi¬ 
crobead  lags  behind  the  front  in  adsorption.  This  is  shown  in  Fig.  5.5(b)  which  depicts  the 
nondimensional  surface  concentration  on  the  microbead  surface  (T /Too)  as  a  projection  on 
the  equatorial  plane  z  =  a  of  the  surface  concentrations  of  the  top  and  bottom  hemispheres. 
For  longer  times  (not  shown),  the  saturation  of  the  surface  decreases  the  binding  rate  and 
the  diffusive  flux  repopulates  target  in  the  diffusion  zone.  As  a  result  the  concentration 


gradients  in  the  zone  decrease,  and  eventually  fade  away.  The  surface  concentration  at  the 
back  end  of  the  microbead  catches  up  with  the  front  end  as  the  surface  concentration  over 
the  entire  microbead  surface  approaches  Too. 

The  binding  curves  in  Fig.  5.5(a)  show  clearly  the  effect  of  Da  at  this  low  (and  fixed) 
Pe.  For  the  highest  value  (Da  =  102),  the  binding  rate  is  fast  relative  to  diffusion  and  the 
sublayer  concentration  at  the  microbead  surface  (cs)  falls  far  below  ca,  and  sharp  concentra¬ 
tion  gradients  develop  in  the  diffusion  zone.  As  a  result,  the  binding  rate  is  much  reduced 
very  far  from  the  kinetic  limit  where  cs  =  cQ.  As  Da  decreases  to  10,  the  kinetic  limit  slows 
down  (in  accordance  with  eq.  5.7),  as  does  the  simulated  binding  rate.  At  Da  =  10,  the 
binding  rate  curve  is  however  closer  to  its  kinetic  limit  because  the  reduced  kinetic  rate  does 
not  deplete  the  sublayer  of  target  as  dramatically  as  when  Da  =  102.  Finally  when  Da  =1, 
the  simulated  curve  is  even  closer  to  the  kinetic  limit,  but  is  still  apart  due  to  the  diffusion 
barriers  due  to  the  large  diffusion  zonesfor  Pe  =  3.5.  In  fact,  a  Da  less  than  0.1  is  required 
for  the  simulated  curves  to  overlap  the  kinetic  curve  for  Pe= 3.5. 

When  Pe  is  increased  to  3500,  a  thin  concentration  boundary  layer  develops  over  the 
microbead  surface  due  to  the  higher  convective  rate  relate  to  the  diffusion  flux,  which  de¬ 
creases  the  diffusional  resistance  to  the  target  transport  by  reducing  the  thickness  of  the 
diffusion  zones,  and  increases  the  rate  of  binding  (Fig.  5.6).  The  diffusion  zone  is  shown  in 
Fig.  5.6(b),  again  for  Da  =  10,  but  for  earlier  nondimensional  times  r  =  .0.012,  2.9  and  5.4 
than  those  shown  for  the  smaller  Pe.  The  boundary  layer  is  very  thin,  except  at  the  very 
backend  of  the  microbead  where  the  slow  convection  through  the  trap’s  backend  opening 
leads  to  some  depletion.  Because  of  the  thin  boundary  layer  circumscribing  most  of  the 
microbead,  the  kinetic  binding  on  the  microbead  surface  is  more  uniform  than  in  the  case 
of  Pe  =  3.5.  The  backend  adjoining  the  trap  opening  lags  only  slightly  behind.  The  surface 
concentration  on  the  lateral  sides  of  the  microbead  which  adjoin  the  trap  side  walls  also  lag 
behind,  due  to  the  reduced  convection  through  the  side  gaps  between  the  trap  walls  and  the 
microbead  which  enlarges  the  boundary  layer.  The  boundary  layer  on  the  underside  is  also 
enlarged  by  the  reduced  convection  between  the  microbead  and  the  bottom  wall  of  the  flow 
cell.  Because  of  this,  and  the  fact  that  the  gap  between  the  top  of  the  channel  and  the  bead 
is  large  (~  20  /im),  and  the  boundary  layer  fits  within  the  gap,  the  top  (z  >  a)  and  bottom 
hemispheres  have  different  average  concentrations  as  shown  in  Fig.  5.7  for  Pe  =  3500  and 
Da  =102,  where  the  more  unobstructed  flow  in  the  top  gap  results  in  a  larger  binding  rate. 
This  difference  becomes  reduced  as  Da  decreases,  and  the  kinetic  limit  is  approached,  since 
the  cublayer  concentration  becomes  uniform.  In  the  case  of  small  Pe  (3.5),  the  difference  in 
binding  rates  between  the  hemispheres  is  not  very  pronounced  because  of  the  large  diffusion 
zone  which  wraps  around  the  trap  creating  an  overall  barrier  to  transport. 
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Figure  5.7:  Asymmetric  binding  to  the  top  and  bottom  hemispheres  of  a  trapped  microbead  (“A”)  due  to 
reduced  convection  at  the  bottom  half  of  the  bead  which  rests  on  the  channel  floor.  Pe  =3500,  Da  =  102, 
Q  =  32,  <Re  =  2xl0-1. 


Nondimensional  Time 


Figure  5.8:  The  average  nondimensional  sublayer  concentration  cs/c0  at  the  surface  of  a  trapped  microbead 
“A”,  as  a  function  of  nondimensional  time  for  Pe  =3500,  Da  =  1,  10,  and  102,  Q  =  32  and  =  2xl0-1. 
The  inset  details  the  sublayer  concentration  immediately  after  the  target  analyte  is  introduced  for  Da  = 
100. 
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(a)  (b) 


Figure  5.9:  The  binding  of  a  target  to  the  surface  of  an  untrapped  microbead  (“A”)  for  Pe  =  3.5  and  Da 
=10  (11  =  32,  93e  =  2.0xl0-4):  (a)  The  average  nondimensional  surface  concentration  on  the  surface  of  the 
microbead,  V /T^,  as  a  function  of  r  for  Da  =10  and  comparison  to  a  trapped  microbead  for  the  same  Pe, 
Da,  fl  and  91c.  (b)  Target  concentration  boundary  layer  for  Da=  10  in  the  equatorial  plane  z  =  a  and  the 
projection  of  the  surface  concentration  averaged  between  the  upper  (z  >  a)  and  lower  (z  <  a)  hemispheres 
of  the  microbead  for  r  =  0.136,  20.4  and  32.6. 


In  Fig.  5.6(a),  r(r)/roo  is  given  for  Da  =  1,  10  and  102  and  Pe  =3500,  and  it  is  clear, 
on  comparison  with  Fig.  5.5(a)  for  Pe  =  3.5  that  binding  rates  are  much  faster  for  the  same 
values  of  Da.  In  this  case,  because  of  the  reduced  diffusional  resistance,  the  binding  curve  at 
Da  =  1  now  overlaps  the  kinetic  curve.  This  overlap  represents  the  case  where  the  binding 
is  controlled  solely  by  the  target-probe  kinetics,  and  is  a  regime  which  has  the  particular 
advantages  when  these  platforms  are  used  to  measure  binding  constants,  since  the  kinetic 
binding  data  only  has  to  be  regressed  against  the  analytic  kinetic  curve  (eq.  5.7),  without  the 
need  to  model  the  convective  diffusion.  In  this  high  Peclet  number  regime,  a  general  criteria 
can  be  developed  through  scaling  arguments  to  identify  the  value  of  Da  below  which  the 
transport  becomes  kinetically  limited.  When  a  thin  concentration  boundary  layer  develops 
around  the  microbead,  the  boundary  layer  thickness,  8  at  the  back  end  can  be  estimated  as 
the  distance  the  target  should  be  away  from  the  surface  for  which  the  time  to  move  along 
the  microbead,  a/ Us  (where  Ug  is  the  velocity  at  a  distance  5  from  the  microbead)  becomes 
equal  to  the  time  required  for  the  target  to  diffuse  to  the  surface,  82 /CD.  The  velocity  can  be 
estimated  as  the  average  velocity  gradient  (or  shear  rate)  evaluated  at  the  microbead  sur- 

D 


face,  7  multiplied  by  8,  and  therefore  =-  =  —  or  -  ~ 

7  o  T>  a 


n  i/3 


a2  7 


For  the  well  studied  case 


of  a  fixed  microbead  in  an  infinite  medium  subject  to  a  uniform  flow  with  concentration 
cQ  far  from  the  bead,  the  average  shear  rate  7  ~  U^/a202.  The  average  shear  rate  over 
the  surface  of  microbeads  in  a  microchannel  either  situated  in  a  trap  or  (ideally)  fixed  to  a 
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surface  is  a  more  complicated  flow  than  the  unbounded  case  due  to  the  effect  of  the  channel 
and  trap  walls.  However,  to  a  first  approximation  in  order  to  construct  a  simple  criterion, 


we  assume  that  7  ~  U/a  so  that 


for  Pe  1.  The  diffusive  flux  to  the  sur¬ 


face,  for  Pe  1,  through  the  boundary  layer  scales  as 


P  {c0  -  cs} 


(where  cs  is  the  average 


sublayer  concentration  on  the  microbead  surface);  equating  this  flux  to  the  the  maximum 
kinetic  flux  (/7T00)  defines  a  scale  for  the  (minimum)  sublayer  concentration  at  the  earliest 
times,  Cs,mm  ~ .  When  Da/Pe1^3  <C  1,  the  sublayer  concentration  is  equal  to  the 

Co  1  +  pel/3 

farfield  concentration,  the  process  is  kinetically  controlled  and  the  binding  is  given  by  eq. 
5.7.  To  validate  this  criterion,  note  first  Fig.  5.8  which  plots  for  microbead  “A”  the  average 
nondimensional  sublayer  concentration  cs  as  a  function  of  nondimensional  time  for  Pe  =3500 
for  Da  =  1,  10,  and  102.  For  the  largest  value  of  Da  the  large  kinetic  binding  rate  relative 
to  diffusion  leads  to  a  rapid  drop  in  the  sublayer  concentration  to  a  value  of  approximately 
0.1,  and  the  sublayer  slowly  recovers  to  one  as  the  surface  saturates  with  target  and  the 
concentration  in  the  diffusion  zone  relaxes.  As  Da  decreases,  the  slowing  of  the  kinetic  rate 
decreases  the  binding,  and  the  initial  reduction  in  the  sublayer  concentration  in  the  earliest 
times  is  much  smaller  than  the  case  for  Da  =  102.  The  scaling  evaluation  for  the  minimum 


sublayer  concentration, 


I  A _ hA. 

1  +  Pe1/ 3 


gives  0.9,  0.6  and  0.1  for  Da  =  1,  10,  and  102, 


respectively,  which  is  in  very  good  agreement  with  the  simulation  minima  as  given  in  Fig. 
5.7.  The  criteria  for  kinetic  control  Da/Pe1^3  <C  1  is  also  validated  in  the  simulations  in 
Fig.  5.6,  where  it  is  clear  that  for  the  three  values  of  Da,  1,  10,  and  102,  Da/Pe1/3  C  1  = 
0.06,  0.6  and  6.6  respectively  and  only  the  Da  =  1  value  satisfies  the  criteria  and  is  in  fact 
the  only  curve  coincident  with  the  kinetic  limit  eq.  (5.7). 

Alternatively  when  Da/Pe1^3  W  1,  the  sublayer  concentration  tends  to  zero,  at  least 
initially,  and  the  target  flux  to  the  surface  is  controlled  solely  by  the  diffusive  mass  trans¬ 
fer.  The  characteristic  time  for  the  target  to  bind  to  an  equilibrium  surface  density  (teq><s) 
T>c0  teq  d  Da 

is  given  by  tep  p  — =-  ~  Too  or  re„  p  =  -  0  ’  ~  where  reQ  p  denotes  a  nondimensional 

0  h^/T)  Pe1'6 

completion  time  scaled  by  the  diffusion  time  across  the  channel  (fp).  For  the  case  of  Da  = 

102  in  Fig.  5.6,  which  is  the  only  value  of  the  Damkohler  number  large  enough  to  satisfy 

Da/Pe1'3  1,  Teq‘T>  ~  6.6  which  is  in  excellent  agreement  with  the  simulation  equilibration 

time.  Bruus  et  al199,  for  a  patch  of  probes  on  one  wall  of  a  two  dimensional  microfluidic 

channel,  have  examined,  by  comparison  to  numerical  simulations,  the  validity  of  reproducing 

the  entire  binding  curve  of  the  target  to  the  patch  for  Pe  1  by  a  boundary  layer  approxi- 
. .  f  cn  -  A 1  .  „  (T  f  )  .  „  ...  .  .  .  ... 


mation  in  which  T 


k  r 

"a*-  00 


—  kdT  with  the  boundary  layer  thickness 


given  by  Pe  1  /3  (see  also  the  two  compartment  models  for  transport  to  surface  probes  in  a 


92 


channel203-211  ).  Bruus  et  al  find  for  large  0.  small  Da  and  large  Pe  the  approximation  is 
most  accurate.  We  find  that  simulating  the  binding  curves  of  Fig.  5.6  for  binding  of  target 
to  microbeads  in  traps  in  our  parameter  space  gives  only  qualitative  agreement,  although  the 
equilibration  is  captured  accurately.  This  is  probably  due  to  the  more  complicated  nature 
of  the  flow  in  the  case  of  the  entrapped  microbeads. 

To  assess  the  effect  of  the  trap  on  the  mass  transfer  of  the  target  to  the  microbead  probe 
surface,  simulations  of  the  binding  are  presented  for  microbeads  situated  in  the  same  array 
as  in  Fig.  5.2,  but  without  the  traps.  Consider  first  the  case  of  Pe  =  3.5,  where  the  diffusion 
zone  is  very  large,  and  in  the  presence  of  a  trap,  extends  beyond  the  trap  itself  as  shown  for 
Da  =  10  (Fig.  5.5(b)).  For  the  very  same  conditions  ( Pe  =  3.5,  Da  =  10,  12  =  32,  and  Tie 
=  2.0  x  10-4)  and  the  same  nondimensional  times  (0.136,  20.4  and  32.6),  Fig.  5.9(a)  shows 
that  the  diffusion  zone  around  the  unobstructed  microbead  is  more  compact.  The  concen¬ 
tration  gradients  at  the  later  times  are  not  as  sharp,  and  less  depletion  is  evident  at  the  back 
end  of  the  microbead.  What  becomes  clear  in  comparing  Figs.  5.5(a)  and  5.9(a)  is  that  the 
trap  obstructs  the  transport  to  the  back  end  completely  at  low  Pe,  leading  to  a  complete 
depletion  in  the  aperture  at  the  back  end  of  the  trap.  When  unobstructed,  the  transport 
is  not  unlike  a  bead  in  an  infinite  medium,  where  some  depletion  at  the  back  end  would  be 
expected  due  to  the  larger  boundary  layer  at  the  downstream  face  of  the  microbead.  The 
nondimensional  surface  concentration  on  the  microbead  surface  (again  a  projection  of  the 
surface  concentration  averaged  between  the  top  (z  >  a)  and  bottom  (z  <  a)  hemispheres) 
shows  a  simple  asymmetry  with  the  back  end  lagginf  behind  the  front  because  of  thinner 
boundary  layer  at  the  upstream  face  of  the  microbead. 

The  more  compact  diffusion  zone  in  the  case  of  open  rather  than  trapped  microbeads 
also  leads  to  a  significantly  more  rapid  binding  rate,  as  shown  in  Fig.  5.9(a)  which  compares 
the  two  transport  geometries  for  Pe  =  3.5.  In  fact,  this  difference  increases  with  the  Da 
number,  as  the  diffusion  limitation  becomes  more  rate  determining  due  to  a  decrease  in  the 
sublayer  concentration.  In  the  opposite  limit,  as  Da  — >  0  and  the  sublayer  concentration 
approaches  the  inlet  concentration,  both  the  trapped  and  untrapped  microbeads  approach 
the  kinetic  limit.  For  this  value  of  Pe,  Da  is  required  to  be  as  small  as  0.1  for  overlap,  so 
differences  exist  for  Da  larger  than  this  value. 

When  the  Pe  number  is  increased  to  3500,  the  diffusion  zone  around  a  trapped  microbead 
takes  the  form  of  a  thin  concentration  boundary  layer  which  hugs  the  microbead  and  lies 
within  the  trap  except  at  the  back  face  and  bottom  part  of  the  microbead  where  the  bead 
contacts  the  trap  walls.  The  presence  of  the  enclosing  trap  is  therefore  not  as  important  in 
the  high  Peclet  number  regime,  and  as  shown  in  5.10,  for  identical  conditions  ( Pe  =  3500, 
12  =  32,  and  Tic  =  2.0  x  10-1),  the  kinetic  binding  rate  to  the  unobstructed  microbeads 
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Figure  5.10:  Comparison  of  the  binding  rate  for  trapped  and  untrapped  microbeads  “A”  in  the  array  in 
the  high  Peclet  number  regime  (Pe  =3500)  for  (a)  Da  =  102  and  (b)  Da  =1.  ($1  =  32  and  =  2.0xl0_1. 


is  only  marginally  faster  at  Da  =  102  (5.10(a))  and  becomes  equal  to  the  binding  rate  for 
the  trapped  microbead  by  Da  =1  (Fig.  5.10(b)),  at  which  point  both  agree  with  the  kinetic 
equation  (eq.  5.7). 

5.3  Experimental  Binding  Assay 

The  aim  of  the  experiments  is  to  first  demonstrate  a  binding  assay  with  the  hydrodynami- 
cally  assembled  microbead  library  platform  of  Fig.  5.1,  using  as  a  prototype  the  conjugation 
of  fluorescently  labelled  NeutrAvidin  as  the  target  analyte  to  biotin  (its  binding  partner)  dis¬ 
played  on  the  microbead  surface,  and  secondly  to  measure  the  surface  kinetic  rate  constant. 
We  illustrate  how  the  flow  rate  of  the  assay  (or  equivalently  the  Pe  number)  can  be  selected, 
using  the  simulations  as  a  guide,  so  that  the  transport  of  target  is  kinetically  controlled, 
and  a  kinetic  rate  constant  can  be  obtained  by  comparing  the  measurements  to  a  simple 
kinetic  expression  (e.g.  eq.  5.7)  without  having  to  solve  for  the  target’s  convective-diffusion 
equations. 

5.3.1  Materials  and  Methods 

The  fluorescently  labeled  protein,  fluorescein  conjugated  NeutrAvidin  (Neutravidin-FITC), 
was  obtained  from  Thermo  Fisher  Scientific,  Inc.  NHS-PEG4-biotin  and  NHS-PEG4  (N- 
hydroxysuccinimide  esters  linked  to  four  polyethylene  glycol  (PEG,  -(CH2CH2)20-)  groups 
and  capped  with  biotin  or  a  methyl  group)  were  purchased  from  Pierce  Corp.  Toluene,  ace¬ 
tone,  isopropanol,  methanol  (99.8%)  and  dimethylformamide  (DMF)  were  obtained  from  Alfa 
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Aesar  (HPLC  Grade).  Phosphate  buffered  saline  buffer  (PBS,  pH=7.4),  hydrogen  peroxide, 
ammonium  hydroxide  and  hydrochloric  acid  were  obtained  from  Fisher  Scientific.  Amino- 
propyltrimethoxysilane  (NH2(CH2)3Si(OCH3)3,  APS)  and  PEG  silane,  2-(methoxy(poly- 
ethyleneoxy) propyl  trimethoxysilane  (CH3(CH2CH20)n(CH2)3Si(0CH3)3,  n=6-9)  were  ob¬ 
tained  from  Gelest.  DI  water  with  a  resistivity  of  18  M-cm  was  obtained  from  a  Millipore 
water  filtration  system.  Negative  tone  SU8-2050  photoresist  and  developer  were  purchased 
from  Microchem.  Sylgard  184  elastomer  base  and  curing  agent  was  purchased  from  Dow 
Corning  as  a  kit.  All  chemicals  were  used  as  received  except  for  the  PBS  which  was  filtered 
with  0.2  /im  syringe  filters.  Soda  lime  glass  calibration  microbeads,  42.3  /mi  in  diame¬ 
ter  with  a  1.5  /rm  standard  deviation  in  the  size  distribution,  were  obtained  from  Thermo 
Fisher  Scientific,  Inc.  Test  grade  single  crystal  polished  silicon  wafers  (3  inch  in  diameter) 
were  obtained  from  Silicon  Quest  Inc.  Polyethylene  tubing  (1.5  mm  OD,  1  mm  ID)  was 
obtained  from  VWR.  Dye  removal  resin  column  kits  were  obtained  from  Thermo  Fisher 
Scientific  (Pierce  Protein  Biology  Products). 

5.3.2  Microbead  Functionalization 

N-hydroxsuccinimide  ester  (NHS)-conjugated  reagents  are  commonly  used  to  label  proteins, 
reacting  with  primary  amine  groups  to  form  an  amide  bond  between  the  protein  and  the 
reagent212,213,  and  this  reaction  is  used  here  to  functionalize  the  surface  of  glass  microbeads 
with  mixed  monolayers  of  PEG  and  biotin  by  first  using  a  silane  reaction  to  attach  amine 
groups  to  the  hydroxylated  surface  of  the  glass  microbeads.  Soda  lime  glass  microbeads 
(42.3  /im  in  diameter)  are  first  washed  in  an  aqueous  solution  consisting  of  4  weight  percent 
hydrogen  peroxide  and  4  weight  percent  ammonium  hydroxide  to  remove  hydrophobic  surface 
contaminants.  The  solution  is  heated  for  10  min.  on  a  hotplate  with  vortexing,  decanted,  and 
the  beads  are  resuspended  in  DI  water.  This  is  repeated  three  times,  and  the  beads  are  then 
washed  similarly  with  an  aqueous  solution  consisting  of  4  weight  percent  hydrogen  peroxide 
and  0.4M  HC1.  After  resuspension  in  DI  water,  the  microbeads  are  dried  in  a  convection  oven 
for  12  hours  at  85°(7.  To  attach  amine  groups  to  the  microbead  surface,  the  washed  and  dried 
microbeads  (30  mg)  are  added  to  10  ml  of  a  5mM  solution  of  aminopropyltrimethoxysilane 
(APS)  in  toluene  for  lhr  30min  under  stirring.  (This  concentration  is  100  times  larger  than 
the  concentration  needed  to  graft  to  a  coverage  of  225  A2  per  amine,  corresponding  to  the 
measured  maximal  amine  grafting  densities,  see  the  discussion  below.)  In  the  presence  of 
water,  either  on  the  surface  or  in  toluene,  the  trijmethoxy  groups  attached  become  hydrolyzed 
to  hydroxyl  groups  which  react  with  the  surface  silanols  to  form  siloxane  bonds  producing 
a  self-assembled  monolayer  (SAM)  terminated  with  a  primary  amine.  To  remove  bounded 
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but  unreacted  APS  from  the  microbead  surface,  the  microbeads  are  washed  and  sonicated 
with  toluene.  They  are  then  centrifuged  and  washed  twice  with  dimethylformamide  (DMF). 

The  probe  molecule,  biotin,  is  attached  to  the  grafted  amine  surfaces  of  the  microbeads 
by  reacting  the  amines  on  the  microbead  surface  with  NHS  derivatized  biotin  (NHS-PEG4- 
biotin)  in  solution  to  form  an  amide  linkage  of  the  biotin  to  the  amine.  The  surface  density 
of  the  biotin  probe  (rp)  can  be  varied  by  reacting  with  mixtures  of  NHS-PEGzrbiotin  and 
NHS-PEG4.  The  resulting  surface  of  the  microbead  displays  PEG-biotin  and  PEG,  with 
equal  subunits  of  PEG,  so  that  the  biotin  extends  just  beyond  the  PEG  layer  making  it 
accessible  to  the  NeutrAvidin  analyte  in  solution.  NHS-PEG  linkers  are  utilized  because 
the  resulting  background  PEG  layer  resists  the  non-specific  adsorption  of  protein214,215  (in 
particular  the  target  NeutrAvidin),  so  that  the  target  recorded  on  the  microbead  surface  is 
only  due  to  the  direct  binding  of  the  biotin  to  NeutrAvidin.  For  mixed  functionalization,  30 
mg  of  microbeads  are  suspended  in  10  ml  of  a  solution  of  NHS-PEGzrbiotin  and  NHS-PEG4 
in  DMF,  with  a  total  concentration  of  the  NHS  ester  of  1  mg/ml  solution,  and  with  the 
molar  ratio  of  NHS-PEG4-biotin  and  NHS-PEG4  varied  with  stirring  overnight.  (This  total 
concentration  is  enough  to  cover  all  the  amine  sites  by  a  factor  of  100.)  The  microbeads  are 
then  washed  three  times  with  DI  water  and  stored  in  a  refrigerator  at  4°C  until  used. 

5.3.3  Device  Design  and  Fabrication 

The  trapping  course  incorporated  in  the  microfluidic  cell,  and  used  to  array  the  microbeads 
into  a  probe  library,  consists  of  a  staggered  array  of  “V”  shaped  enclosures  as  shown  in 
Fig.  5.2  and  directly  modeled  in  the  simulations.  The  trapping  course  is  set  in  a  wide 
channel  of  height  60  /mi  which  contains  an  open  space  upstream  for  an  inlet  port  to  inject 
the  microbeads.  The  wide  channel  is  connected  upstream  by  a  narrow  channel  to  an  inlet 
for  introducing  the  analyte,  and  connected  downstream  to  an  opposing  narrow  channel  to  an 
outlet  for  collecting  microbeads  and  analyte.  The  cell  is  fabricated  in  two  layers,  a  bottom 
part  containing  the  features  of  the  trapping  array,  connecting  network  and  inlet  and  outlet 
ports  which  is  imprinted  on  a  monolith  by  molding  against  a  master.  The  second  part  is  a 
glass  slide  which  is  bonded  to  the  face  of  the  molded  bottom  part  containing  the  features  to 
enclose  the  channels  and  the  trapping  array.  The  molded  part  is  made  from  the  elastomer 
polydimethylsiloxane  (PDMS)  using  soft  lithography.  The  master,  with  the  negative  features 
of  the  fluidic  network,  is  made  from  a  photoresist  using  lithography.  The  bottom  part  is 
molded  from  the  master  using  the  two  component  heat  curing  system  from  Dow  Corning 
Sylgard  184;  the  ratio  of  base  to  elastomer  is  10:1  by  weight.  The  photoresist  master  is  placed 
at  the  center  of  a  Petri  dish,  and  the  base-curing  agent  mixture  is  poured  over  the  master, 
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and  allowed  to  sit  under  vacuum  for  30  min  to  release  bubbles  sticking  to  the  bottom  of  the 
master  holes  that  would  eventually  become  the  PDMS  trap.  The  mixture  is  then  cured  by 
heating  in  an  oven  at  65° C  for  one  hour.  The  resulting  elastomer  is  subsequently  peeled  off 
of  the  master,  and  cut  into  individual  cells.  Inlet  and  exit  ports  are  punched  through  the 
molded  face  to  the  backside  of  the  monolith,  and  an  access  hole  is  also  punched  in  the  space 
upstream  of  the  trapping  array  for  the  introduction  of  the  microbeads.  The  elastomer  and 
glass  slide  are  sealed  by  bringing  the  faces  in  conformal  contact  to  form  siloxane  linkages 
after  oxygen  plasma  treating  their  open  faces  (Harrick  Plasma  Cleaner).  The  inside  surfaces 
of  the  flow  channels  and  the  interior  surfaces  of  the  traps  are  functionalized  with  a  PEG 
silane  to  prevent  nonspecific  adsorption. 

5.3.4  Microbead  Arraying 

Microbeads  are  arrayed  in  the  microfluidic  device  by  introducing  them  through  the  injection 
port  in  front  of  the  traps.  The  prepared  device  is  submerged  in  a  beaker  of  water,  which  is 
placed  in  a  dessicator  chamber  under  mild  vacuum  to  extract  air  bubbles  in  the  device.  After 
the  air  bubbles  have  been  removed  by  visual  inspection,  the  device  is  kept  under  water.  A 
2.0  ml  aliquot  of  a  microbead  suspension  is  loaded  in  a  glass  gas  tight  syringe,  and  tubing  is 
slipped  over  a  needle  at  the  end  of  the  syringe.  The  microbeads  are  allowed  to  settle  back 
against  the  plunger  of  the  syringe  with  the  syringe  orientated  horizontally,  and  the  tubing 
extending  from  the  syringe  is  then  filled  with  water  by  depressing  the  plunger.  The  tubing 
is  then  inserted  into  the  access  hole  above  the  staging  area  in  front  of  the  trapping  array. 
The  syringe  is  maintained  in  a  horizontal  position  so  no  beads  are  yet  introduced  into  the 
device.  The  inlet  port  of  the  cell  is  then  connected  to  a  target  analyte  filled  syringe  by 
inserting  the  tubing  into  the  access  port  on  the  device  and  slipped  over  a  needle  attached 
to  the  syringe.  Tubing  is  inserted  into  the  exit  hole.  When  all  the  ports  are  connected  to 
the  tubing,  the  device  is  gently  removed  from  the  water,  dried  with  a  wipe  and  placed  on 
the  stage  of  a  optical  microscope.  The  syringe  containing  the  analyte  target  is  placed  in  a 
syringe  pump  (Harvard  Apparatus  PHD).  The  microbeads  are  introduced  into  the  device  by 
raising  the  syringe  to  a  vertical  orientation  to  allow  microbeads  to  settle  at  the  bottom  of  the 
syringe.  The  syringe  plunger  is  then  depressed  to  push  the  microbeads  through  the  obstacle 
course,  delivering  0.5  ml  in  approximately  30  sec  for  a  flow  rate  of  103  /d/ min.  The  filling 
of  the  course  with  the  microbeads  is  observed  through  optical  microscopy  by  visual.  After 
injection,  the  syringe  is  placed  back  in  a  horizontal  position,  the  beads  return  to  the  side 
and  back  end  of  the  syringe  and  the  plunger  is  depressed  slowly  to  flow  DI  water  through 
the  course  to  remove  the  few  beads  that  are  untrapped. 
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5.3.5  Screening  Assay  and  Fluorescence  Imaging  of  NeutrAvidin- 
FITC  Binding  to  the  Biotinylated  Microbeads 

NeutrAvidin  is  a  homotetrameric  protein  (MW  =  60000),  consisting  of  four  identical  sub¬ 
units,  each  of  which  contains  a  receptor  pocket  site  for  the  insertion  of  the  small  molecule 
biotin  ligand  (MW  =  244).  NeutrAvidin  is  the  deglycosylated  form  of  the  naturally  occur¬ 
ring  protein  Avidin,  with  a  near  neutral  pi  equal  to  6.3.  X-ray  studies  indicate  that  the 
sites  are  barrel  shaped  enclosures  constructed  from  a  helices,  and  arranged  in  the  tetramer 
so  that  two  sites  are  on  one  face  of  the  protein,  and  two  sites  are  on  an  opposite  face,  with 
each  face  having  an  area  of  At  =  3025  A.  Because  of  the  large  area  of  the  binding  face  of 
the  protein,  the  insertion  of  the  surface  biotin  into  the  barrel  site  can  cover  other  surface 
ligands  and  make  them  unavailable  for  binding  another  protein  if  the  biotin  probe  density 
Tp  is  larger  than  the  inverse  of  the  protein  binding  face. 

The  biotin  probe  density  on  the  surface  of  the  microbeads  is  determined  firstly  by  the 
density  of  the  amine  groups  grafted  to  the  surface  by  the  silane  reaction,  and  secondly  by 
the  ratio  of  the  reagents  NHS-PEG4-biotin  and  the  NHS-PEG4  which  are  mixed  in  solution 
and  react  with  the  surface  amines  to  link  either  PEG4  or  PEG4-biotin  to  the  surface  via  an 
amide  bond.  The  surface  density  of  APS-silane  grafted  amines  on  glass  which  are  available 
for  reaction  with  NHS  has  been  measured  in  2.19(b)  by  first  reacting  NHS,  linked  to  a 
disulfide,  to  surface  grafted  amines  on  glass  slides,  and  then  attaching  a  fluorescent  group 
to  the  disulfide.  This  group  is  then  removed  into  solution  by  reduction,  and  the  fluorescence 
in  solution  is  measured  to  obtain  the  number  of  amine  groups  which  reacted  with  the  NHS- 
disulfide.  They  find  that  aminated  surfaces  prepared  by  APS  silanization  in  a  toluene  solvent 
(as  here)  results  in  surface  densities  of  approximately  20  amine  sites  per  100  nm2  (500 
A2 /site).  In  a  second  study,  the  surface  density  of  APS-silane  grafted  surface  amines  on  silica 
nanoparticles  was  measured  by  directly  reacting  the  amine  with  a  isothiocyanate  linked  to  a 
fluorescent  group,  and  then  measuring  the  fluorescence  of  the  nanoparticles  to  quantify  how 
many  fluorescent  molecules  had  bound  to  the  surface  of  the  particles.  In  this  case  the  silanes 
where  grafted  onto  the  particles  from  an  ethanol  solvent,  and  they  find  a  larger  surface 
density  of  44  amine  sites  per  100  nm2  (227  A2 /site).  In  either  case,  the  surface  density  is 
much  smaller  than  the  theoretical  limit  of  the  maximum  packing  for  the  grafted  amines, 
which,  based  on  the  cross  section  of  the  propyl  chain  of  APS  in  an  all  trans  configuration, 
is  20  A2,  or  500  amine  sites  per  100  nm2.  (This  is  also  approximately  the  surface  density 
of  hydroxyl  groups  on  silica  2.19(b).)  The  smaller  surface  densities  of  the  grafted  amines 
compared  to  the  maximum  packing  is  attributed  to  the  reduced  intermolecular  cohesion  of 
the  propyl  group  which  does  not  allow  for  a  high  surface  density  self  assembled  monolayer 


98 


(SAM)  to  form,  and  hydrogen  binding  of  grafted  amines  to  surface  hydroxyl  groups  which 
blocks  the  grafting  of  additional  APS  to  the  surface  2.19(b). 

In  our  silane  functionalization  of  the  glass  microbeads,  the  concentration  of  the  APS 
silane  solution  uses  is  100  times  larger  than  the  concentration  necessary  to  cover  the  mi¬ 
crobeads  assuming  the  maximum  attainable  coverage  of  44  sites  per  100  nm2  or  227  A2  per 
amine,  so  we  assume  this  maximum  amine  coverage  is  achieved.  Although  this  coverage  is 
low,  the  area  occupied  per  amine  molecule  is  still  much  larger  than  the  binding  face  of  Neu- 
trAvidin,  so  that  if  every  amine  site  was  linked  to  a  biotin  by  the  hydroxysuccinimide  reaction 
with  NHS-PEC>4-biotin,  protein  bound  to  a  surface  biotin  would  block  or  stericaly  hinder 
other  unoccupied  biotin  probe  sites  from  binding  protein.  (Note  that  the  NHS-PEC^-biotin 
links  to  the  surface  amine,  the  oligomeric  ethylene  oxide  chain  forms  a  helical  configuration 
from  the  surface  amide  bond  with  a  cross  section  of  approximately  20  A2  per  molecule, 
which  is  much  smaller  than  the  area  per  grafted  amine,  so  the  PEO4  does  not  overlap  other 
surface  amines.)  Thus  to  model  the  binding  kinetics  of  NeutrAvidin  to  biotinylated  glass 
microbeads  with  the  highest  probe  densities  (Tp  =20  -  44  biotin  molecules/100  nm2),  the 
effect  of  steric  hindrance  should  be  incorporated  in  eq.  5.4,  as  the  fraction  of  biotin  sites 
available  for  binding  ,  0,  is  not  equal  to  the  fraction  of  unbound  sites.  If  the  biotin  probe 
density  is  reduced  so  that  the  area  per  biotin  probe  is  much  larger  than  the  binding  face  of 
NeutrAvidin  (i.e.  l/rp  >  At),  then  when  the  protein  binds  to  the  surface  biotin  it  does  not 
obstruct  unbound  probes,  and  the  available  fraction  for  binding  is  equal  to  the  fraction  of 
unbound  probes  (Langmuir  kinetics,  0  =  1  —  and  the  maximum  surface  concentration  of 
the  target,  Too,  equals  the  probe  density  Pp.  Here  we  study  this  low  probe  density  limit  to 
insure  Langmuir  kinetics,  and  lower  the  biotin  probe  density  by  linking  the  amines  on  the 
surface,  through  the  hydroxysuccinimide  reaction,  to  PEG4  as  well  as  PEG4-biotin  by  re¬ 
acting  the  microbeads  with  a  mixture  of  NHS-PEG4  and  NHS-PEG4-biotin  in  a  molar  ratio 
of  1:1000(NHS-PEG4  to  NHS-PEG4-biotin).  These  molecules  are  identical,  except  for  the 
biotin  moiety,  which  is  distal  to  the  NHS  reacting  center,  and  therefore  the  reaction  rates  to 
the  surface  amines  should  be  proportional  to  the  bulk  concentrations,  As  such,  the  ratio  of 
the  surface  concentration  of  PEG4  linked  to  the  surface  amines,  to  the  surface  concentration 
of  linked  PEG4-biotin  should  be  in  the  ratio  of  their  bulk  concentrations.  Therefore  the 
surface  probe  density  of  the  biotin  is  1/1000  of  the  density  of  the  amine  sites,  or  a  (biotin) 
area  per  molecule,  which  is  1000  times  larger  than  the  amine  area  per  molecule  of  227  -  500 
A2,  which  in  turn  is  much  larger  than  the  binding  face  of  NeutrAvidin. 

The  target  analyte,  NeutrAvidin,  fluorescently  labeled  with  fluorescein  (NeutrAvidin- 
FITC),  is  flowed  through  the  cell,  and  the  accumulating  fluorescence  on  the  surface  of  the 
microbeads  functionalized  with  the  biotin  surface  probes  are  obtained  using  epifluorescence 
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measurements  on  an  inverted  Nikon  Eclipse  TIE  microscope  with  a  motorized  stage  for 
scanning  different  areas  of  the  trapping  array,  a  10x/NA=0.3  air  objective,  a  filter  cube  for 
measurements  of  FITC  luminescence  (excitation  filter  470  ±  40  nm,  emission  (barrier)  filter 
525  ±  50  nm)  and  a  high  sensitivity  digital  camera  (DigiSight,  Nikon).  The  microfluidic  cell 
is  placed  with  the  glass  slide  down,  so  that  the  array  is  observed  through  the  optically 
more  transparent  thin  glass  slide  than  the  PDMS  half  of  the  monolith.  Images  of  the 
microbead  array  fluorescence  at  different  regions  in  the  assembled  array  were  obtained  in 
regular  10  min  time  intervals  from  a  starting  point  (f=0)  when  the  analyte  flow  is  first 
introduced  into  the  array  for  a  period  of  approximately  1.5  hr.  The  NeutrAvidin-FITC 
flow  rate  is  100  /il/min  and  the  concentration  of  protein  is  500  ng/ml,  and  images  were 
recorded  for  approximately  1.5  hr  which  is  the  observed  time  necessary  for  the  protein 
binding  to  saturate  and  the  fluorescence  on  the  microbead  surface  to  become  steady.  NIH 
imaging  software  is  used  on  the  time  sequence  of  images  of  the  array  luminescence  to  encircle 
the  fluorescence  on  individual  microbeads  in  the  array,  and  sum  the  recorded  intensities  of 
the  pixels  in  the  encircled  regions  to  obtain  the  individual  fluorescence  of  each  encircled 
microbead.  For  each  microbead  whose  fluorescence  is  recorded  in  this  way,  the  integrated 
intensity  at  each  time  is  divided  by  the  steady  luminescence  of  the  microbead,  obtained 
by  averaging  the  integrated  intensity  of  the  microbeads  over  the  last  few  minutes  of  the 
1.5  hr  of  the  experiment,  where  the  intensity  appeared  approximately  constant.  Assuming 
the  surface  concentration  of  the  labelled  protein  is  linear  in  the  measured  fluorescence,  the 
quotient  of  the  measured  microbead  fluorescence  at  time  t  to  the  asymptotic  saturation  value 
for  long  time  is  equal  to  the  average  surface  concentration  (f  (t))  divided  by  the  maximum 
packing  concentration  (Too).  Each  experiment  involves  the  recording  of  tho  the  fluorescence 
for  several  microbeads  at  each  time  interval;  these  individual  recordings  (divided  by  the 
asymptotic  (saturation)  values  for  each  microbead)  are  averaged  and  plotted  as  Too  as  a 
function  of  time,  with  error  bars  established  from  the  standard  deviations  of  the  individual 
microbead  measurements.  The  result  for  the  0.1%  NHS-PEG4-biotinylated  microbeads  (from 
three  parallel  experiments)  is  plotted  in  Fig.  5.11  (blue  solid  line)  along  with  the  error  bars 
obtained  from  the  standard  deviations  of  the  measurements. 

5.3.6  Comparison  with  Simulations 

Simulations  are  constructed  to  be  in  congruence  with  the  experiments  and  follow  the  rep¬ 
resentation  of  the  array  as  presented  earlier  in  which  one  unit  cell  of  four  rows  is  modelled 
(Fig.  2).  The  physical  parameters  of  the  analyte  solutions  are  set  to  the  same  as  the  exper¬ 
iment:  NeutrAvidin  diffusion  coefficient,  D  =  6xl0-11m2sec-1,  Cq  =  8.33  xlO-6  mol  m-3 
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Figure  5.11:  A  comparison  of  simulation  data  with  experiment  data  for  the  dynamic  response  of  the 
microbead  array.  The  blue  solid  line  with  error  bars  shows  the  normalized  intensity  of  the  fluorescence  of 
NeutrAvidin-FITC  bound  to  0.1  mol  %  NHS-  PEG4-Biotin  microbeads  as  a  function  of  time  for  a  solution 
concentration  of  500  ng/ml  and  a  flow  rate  of  100  /A/ min.  The  red  line  shows  the  simulation  fit,  for  a  value 
of  ka  =  180  m3  mole-1  sec-1,  and  the  kinetic  limit  (red)  overlies  the  finite  element  simulation. 

and  screening  assay  perfusion  flow  rate  of  102/il/min  (average  velocity  U  =  10_2m/sec).  As 
we  noted,  for  a  1/1000  molar  ratio  of  NHS-PEG4  —  biotin  to  NHS-PEG4  the  molecular  sat¬ 
uration  coverage  of  the  protein  is  equal  to  the  surface  density  of  the  biotin,  and  the  kinetics 
can  be  described  by  a  Langmuir  model  in  which  the  bound  proteins  do  not  interact.  For  0.1 
mol  %  NHS-PEG4-Biotin,  assuming  an  area  per  molecule  of  the  amine  binding  sites  of  200 
AA2  ,  r qq  8.3  X  10  TfXoli  *  TTl  .  TThe  Peclet  number  P C  based  OH  the  mdlUS  of  the  bead 
(Pg  =  Ua/D)  is  equal  to  3500,  and  0  based  on  the  initial  concentration  of  NeutrAvidin-FITC 
(Q  =  is  equal  to  4.8.  The  Damkohler  number  Da  =  fcrg’ct  is  then  varied  until  they  are  in 
agreement  with  the  experiments  and  we  find  Da  =  0.05  which  corresponds  to  ka  =  180  m3 
mol-1  s_1,  (see  Fig.  5.11  (red  dot  line)).  Also  in  the  figure  is  shown  the  purely  kinetically 
controlled  curve,  which  as  expected  from  the  general  simulations  for  a  value  of  Da  much  less 
than  one  and  Pe  of  order  103,  overlap  the  mixed  simulation  curve  as  the  process  under  these 
conditions  is  kinetically  controlled. 

5.3.7  Summary  and  Conclusions 

In  this  study,  we  have  examined  the  mass  transfer  of  a  target  molecule  to  a  microbead 
probe  library  which  was  incorporated  in  a  wide  channel  of  a  microfluidic  cell.  The  library  is 
constructed  by  conjugating  the  probe  molecules  to  the  microbead  surface,  and  assembling  the 
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microbeads  in  an  array  layout  by  hydrodynamically  capturing  the  microbeads  in  a  trapping 
course  situated  in  the  channel.  The  traps  are  V-shaped,  half-open  enclosures  with  the  open 
end  facing  the  flow,  and  are  arranged  in  a  staggered  configuration.  The  enclosures  span  the 
entire  height  of  the  channel,  which  is  just  slightly  larger  than  the  microbead  diameter,  and 
are  sized  to  sequester  only  one  microbead,  so  that  the  trap  course  arrangement  serves  as  the 
template  for  arraying  the  microbeads,  and  keeping  them  in  place.  The  traps  have  an  aperture 
opening  at  the  back  end  to  allow  flow  over  the  surface  of  the  microbead.  Assaying  the  binding 
interaction  of  a  target  with  the  probe  library  is  undertaken  by  flowing  a  fluorescently  labelled 
target  analyte  through  the  course,  and  identifying  microbeads  which  are  luminescent. 

The  central  issue  which  we  address  is  the  influence  of  the  retaining  structure  of  the  trap 
on  the  accessibility  of  the  surface  probes  on  the  microbead  surface  to  the  target  analyte 
streaming  over  the  microbeads.  The  dynamic  response  of  the  microbead  array  -  i.e.  the  time 
required  to  develop  a  recognizable  signal  (or  equivalently  the  number  of  fluid  cell  volumes 
for  this  signal)  is  dependent  on  the  ability  of  the  target  to  diffuse  rapidly  enough  across 
streamlines  to  the  probe  surface.  Once  at  the  surface,  the  target  is  required  to  kinetically 
conjugate  to  the  probes  and  this  rate,  as  an  in-series  process,  also  contributes  to  the  overall 
binding  rate  of  the  target  to  the  probe.  We  undertook  numerical  simulations  of  the  mass 
transfer,  using  a  Langmuir  model  to  describe  the  surface  conjugation  step,  and  computed 
the  surface  concentration  of  the  target  on  the  microbead  surface  as  a  function  of  time.  This 
overall  binding  rate,  in  nondimensional  form,  was  characterized  by  two  groups,  a  Peclet 
number  (Pe)  describing  the  rate  of  target  solute  convection  across  the  microbead  (which  is 
proportional  to  the  analyte  flow  rate)  to  the  rate  of  diffusion  to  the  particle  surface,  and  a 
Damkholer  number  (Da)  describing  the  rate  of  target-probe  surface  conjugation  to  the  rate 
of  bulk  diffusion  to  the  surface.  We  find  that  the  effect  of  the  trap  on  the  transport  rate 
can  be  described  through,  and  is  critically  dependent  on,  the  Peclet  number.  For  fixed  Da 
and  order  one  values  of  Pe  (corresponding  to  weak  convection),  as  shown  in  Fig.  5.5,  large 
diffusion  zones,  where  the  concentration  is  below  the  inlet  concentration,  develop  around 
the  microbead  and  envelop  the  trap.  In  this  case,  because  of  the  presence  of  the  trap,  the 
weak  convection  cannot  bring  sufficient  amounts  of  fresh  analyte  to  the  probe  surface  at 
a  rate  comparable  to  the  rate  of  diffusion  to  the  surface.  As  a  result  the  overall  rate  of 
target-probe  binding  is  slow,  and,  as  shown  in  Fig.  5.9  is  much  slower  that  the  rate  at  which 
target  binds  to  the  surface  in  the  idealized  case  of  an  unobstructed  microbead  at  the  floor  of 
the  microchannel.  Solutions  for  the  fluid  flow  also  demonstrate  the  attenuation  of  the  flow 
around  the  microbead  -  especially  in  the  space  between  the  sidewalls  of  the  trap  and  the 
microbead  -  relative  to  the  idealized  unobstructed  microbead.  For  the  same  value  of  Da, 
and  Pe  of  order  103,  we  find  the  very  interesting  result  that  the  larger  convection  at  higher 
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Pe  leads  to  the  development  of  diffusion  zones  which  take  the  form  of  thin  boundary  layers 
which  encircle  the  beads  within  the  trap  (Fig.  5.6).  Because  the  boundary  layers  are  thin, 
the  diffusion  fluxes  are  large  and  the  overall  binding  rate  increases  relative  to  Pe  of  order 
one.  In  fact,  the  binding  rate  approaches  the  rate  for  the  idealized  case  of  an  unobstructed 
microbead  positioned  at  the  bottom  of  the  channel.  Thus  if  flow  rates  are  large  enough, 
the  presence  of  the  trap  retaining  structure  necessary  to  hold  the  microbeads  in  place  does 
not  necessarily  have  to  hinder  the  target  transport  to  the  surface.  In  the  trap  geometry 
used  here,  the  aperture  at  the  back  end  of  the  trap  (and  the  gap  between  the  microbead 
and  the  channel  top)  allows  flow  over  the  microbeads  (see  the  flow  simulation,  Fig.  5.3 
and  Fig.  5.4)  and  it  is  this  flow,  when  Pe  is  large  enough,  that  allows  the  development  of 
the  thin  boundary  layers  which  sit  inside  the  trap  and  around  the  microbeads  and  enhance 
the  diffusion  rate  of  the  target  to  the  surface.  In  other  studies,  wells  at  the  floor  of  the 
microchannel  are  used  to  sequester  and  array  microbeads  to  construct  the  probe  library. 
In  this  geometry,  target  only  streams  over  the  top  part  of  the  microbead,  and  at  high  Pe 
thin  boundary  layers  cannot  encircle  a  large  portion  of  th  emicrobead.  As  a  result  diffusive 
transport  is  hindered,  and  one  solution  which  is  proposed  is  to  create  drains  in  the  wells  to 
enhance  convection  around  the  microbead.  The  trapping  geoometry  studied  here  uses  an 
aperture  at  the  back  end  to  accomplish  the  same  goal  (as  well  as  providing  a  gap  between  the 
microbead  and  the  top  surface  of  the  channel),  and  is  much  easier  to  implement.  Another 
method  for  constructing  a  microbead  array/  probe  library  in  a  microfluidic  cell  is  to  place 
the  functionalized  microbeads  in  shallow  wells  arrayed  at  the  bottom  of  the  channel,  and 
then  place  the  channel  top  on  the  microbead  array.  The  height  of  the  channel  is  the  same 
size  as  the  microbead  diameter,  so  that  the  channel  top  holds  the  microbeads  in  place.  The 
binding  rate  in  this  geometry  approaches  the  idealized  case  of  an  unobstructed  microbead 
in  the  channel  for  all  Pe,  but  the  microbead  array  assembly  is  undertaken  ex-situ,  which  is 
more  demanding  than  the  in-situ  hydrodynamic  entrapment  of  the  trap  array  studied  here. 

The  surface  kinetic  parameter  Da  is  determined  by  the  nature  of  the  target  and  the 
probe,  and  the  probe  surface  density,  and  is  independent  of  the  trapping  geometry.  We  find 
that  as  Da  decreases  for  a  fixed  Pe,  the  reduced  kinetic  conjugation  rate  relative  to  diffusion 
results  in  more  uniform  surface  concentrations  around  the  microbead  as  the  diffusion  zones 
disappear.  We  identified  for  each  Pe  values  of  Da  for  which  the  surface  kinetics  controls  the 
overall  binding  rate,  and  the  surface  concentration  can  be  described  by  a  simple  exponential 
solution  which  represents  the  integration  of  the  Langmuir  kinetic  equation  for  a  constant 
sublayer  concentration  (eq.  5.7).  The  larger  Pe,  the  larger  is  the  critical  value  of  Da  to 
achieve  this  kinetic  limit,  as  the  larger  Pe  have  reduced  diffusion  limitations.  We  find  in 
particular  for  the  trapping  geometry,  that  for  Pe  of  order  103,  values  of  Da  less  than  or  equal 
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to  one  insures  the  kinetic  limit. 
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Chapter  6 

Electrocoalescence  of  Water-in-Crude 
Oil  Emulsions  in  Two  Dimensions 


Direct  observation  of  electrocoalescence  of  water  droplets  in  crude  oil  emulsions  is  demon¬ 
strated.  Two  dimensional  configurations  of  monodisperse  water  droplets  are  subjected  to  a 
uniform,  external  electric  field  and  coalescence  of  individual  droplet  pairs  is  observed.  By 
calculating  the  electrostatic  force  between  droplet  pairs  using  finite  element  simulations,  the 
identity  of  coalescing  pairs  can  be  accurately  predicted. 

6.1  Introduction  and  Background 

Emulsion  stability  is  relevant  to  a  wide  range  of  applications,  including  foods,  cosmetics, 
petroleum  and  other  industrial  processes.  The  fundamental  issue  in  the  understanding  of 
emulsion  stability  is  the  interaction  between  droplets  (or  bubbles)  of  the  dispersed  phase, 
specifically  the  time  required  for  them  to  approach  and  coalesce.  This  can  be  understood 
in  terms  of  the  time  required  to  drain  the  continuous  fluid  phase  from  between  adjacent 
droplets.  To  predict  this  time,  knowledge  of  the  forces  between  the  droplets  is  required 
to  make  the  coalescence  time  a  well-defined  hydrodynamics  problem.  Validation  of  such  a 
prediction  is  most  accurately  realized  by  direct  observation  of  individual  droplet  pair  inter¬ 
actions,  but  this  can  be  difficult  to  achieve  from  experiments  performed  on  bulk  emulsions. 

Microfluidics  offers  the  ability  to  generate  emulsions  of  monodisperse  droplets  and  observe 
interactions  between  droplets  on  a  pairwise  basis 102,216-218 .  Droplet  separation  distances  can 
be  precisely  measured,  enabling  calculation  of  droplet  interaction  forces  on  a  pairwise  basis. 
The  origin  of  these  interaction  forces  depends  on  the  emulsion  system  being  studied  and  its 
application.  In  some  instances,  emulsion  stability  is  desired,  while  other  applications  are 
predicated  on  efficient  separation  of  the  emulsified  phases.  In  these  applications,  an  external 
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field  is  often  applied  to  facilitate  droplet  coalescence  and  separation. 

A  primary  example  of  this  field-assisted  separation  is  the  use  of  dielectrophoretic  forces  to 
electrocoalescence  conducting  water  droplets  in  a  insulating  continuous  oil  phase219,220.  The 
importance  of  electrocoalescence  is  primarily  due  to  its  utilization  in  the  oil  industry,  where  it 
is  used  to  separate  water  from  crude  oil  in  a  process  colloquially  known  as  ’desalting’.  Water 
droplets  intentionally  introduced  into  the  crude  oil  to  extract  salt  are  subsequently  removed 
via  application  of  an  electric  field,  which  polarizes  the  conducting  droplets,  producing  attrac¬ 
tive  (and  repulsive)  forces  between  droplets  that  result  in  coalescence  of  droplet  pairs221-230. 
Successive  coalescence  events  between  droplets  produce  progressively  larger  droplets  that 
settle  under  gravity  into  a  bulk  water  phase  that  can  be  easily  removed.  The  simplest  model 
of  electrocoalescence  assumes  that  the  water  droplets  behave  as  infinitely  conducting  spheres 
in  a  perfectly  insulating  oil,  and  that  the  interfacial  rheology  of  the  water-crude  oil  interface 
plays  no  role  in  droplet  coalescence.  In  practice,  many  crude  oils  have  significant  conductiv¬ 
ity,  reducing  droplet  polarization  when  the  electric  field  is  applied  at  0  Hz.  More  importantly, 
crude  oils  contain  numerous  polar  compounds  known  as  asphaltenes,  which  adsorb  at  the 
oil-water  interface  and  act  to  stabilize  the  droplets,  inhibiting  coalescence.  Due  the  unique 
nature  of  the  asphaltene  constituents  in  a  specific  crude  oil,  the  effect  of  these  surface-active 
compounds  cannot  be  ascertained  independently.  A  valid  experimental  investigation  should 
therefore  be  conducted  using  the  crude  oil  under  study  to  accurately  predict  electrocoales¬ 
cence. 

Numerous  studies  of  various  aspects  of  electrocoalescence  have  been  published.  Sev¬ 
eral  examine  the  physics  of  the  coalescence  event  itself.  Others  examine  the  behavior  of 
droplet  pairs  initially  separated  by  some  distance.  Bibette  et  al.231  measured  the  effect  of 
field  strength  and  droplet  separation  distance  on  the  coalescence  of  microfluidically  gener¬ 
ated  surfactant-stabilized  water  droplet  pairs  in  hexadecane,  mapping  out  a  phase  diagram 
which  includes  three  regimes:  coalescing,  non-coalescing  and  partial  coalescing.  Chiesa  et 
ap  226-228  moc[eiec[  the  trajectory  of  a  small  droplet  coalescing  into  a  larger  one  as  a  force 
balance  between  the  attractive  electrostatic  and  resistive  hydrodynamic  forces,  showing  good 
agreement  between  experiment  and  theory.  However,  no  current  study  examines  electrocoa¬ 
lescence  in  two  dimensions,  with  interactions  between  multiple  droplets  taken  into  account. 
Furthermore,  all  of  the  existing  studies  are  performed  in  model  systems,  which  are  of  limited 
applicability  to  water  in  crude  oil  emulsions  where  the  heterogeneous  nature  of  the  crude  oil 
may  have  a  dramatic  impact  on  electrocoalescence. 

This  study  reports  results  of  electrocoalescence  experiments  performed  on  two  dimen¬ 
sional  configurations  of  monodisperse  water  droplets  in  a  crude  oil.  The  objective  was  to 
first  use  microfluidics  to  visualize  the  electrocoalescence  process,  and  secondly  to  accurately 
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(a)  (b) 


Figure  6.1:  Schematic  of  the  microfluidic  cell  used  in  the  experiments  (a)  Arrangement  of  the  flow  focusing 
orifice  and  the  chamber  to  collect  the  droplets  and  apply  the  electric  field,  (b)  Detail  of  the  flow  focusing. 


predict  electrocoalescence  between  droplet  pairs  based  on  a  calculation  of  the  pairwise  elec¬ 
trostatic  force  as  determined  by  finite  element  simulations. 

6.2  Experimental 

To  enable  direct  observation  of  water  droplets  in  crude  water-in-oil  emulsions,  experiments 
were  performed  in  a  microfluidic  cell  composed  of  transparent  PDMS.  The  emulsion  was 
formed  by  using  flow  focusing 105,232,233  to  create  a  monodisperse  train  of  water  droplets 
in  crude  oil,  which  were  then  directed  through  a  feeding  channel  to  a  holding  chamber  in 
which  they  arrange  themselves  in  an  arbitrary  configuration,  and  in  which  the  electric  field 
is  applied  (Fig.  6.1).  Droplets  flow  out  of  the  chamber  through  an  exit  port.  The  height 
of  the  channels  in  the  cell  is  100  //m,  and  the  droplets  produced  from  the  50  //in  orifice  are 
approximately  80  /an,  so  the  droplets  move  through  the  feeding  channel  and  the  chamber  as 
a  single  layer.  The  microchannels  were  fabricated  using  standard  soft  lithography  techniques 
to  produce  two  layers  of  PDMS,  one  containing  the  fluidic  channel  and  a  second,  flat  layer 
to  seal  the  channel.  The  two  layers  were  bonded  together  following  exposure  to  an  oxygen 
plasma  and  mounted  on  a  standard  glass  microscope  slide.  Access  ports  cored  into  the  top 
layer  using  a  biopsy  punch  allowed  introduction  fluids  of  interest  via  polyethylene  tubing 
connected  to  syringe  pumps.  To  generate  a  uniform  electric  field  £  across  the  droplets, 
planar  electrodes  of  aluminum  were  inserted  into  the  PDMS  perpendicular  to  the  lateral 
plane  of  the  cell,  and  sited  externally  from  the  fluid.  The  electrodes  were  connected  to  an 
amplifier  controlled  by  a  frequency  generator. The  electric  field  is  applied  across  the  channel 
and  the  droplet  coalescence  events  are  recorded  using  a  high-speed  camera.  To  prevent  the 
electrical  conductivity  of  the  crude  oil  from  reducing  the  electric  field  over  time,  the  electric 
field  is  applied  at  a  frequency  of  500  Hz.  The  resulting  coalescence  process  in  the  chamber 
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(a) 


(b) 


Figure  6.2:  (a)  Flow  focusing  production  of  water  droplets  in  oil  and  (b)  assemblage  of  droplets  in  the 
chamber. 


is  visualized  by  optical  microscopy  using  an  inverted  microscope  in  brightfield  mode  and  a 
high  speed  camera. 

Optical  micrographs  of  the  flow  focusing  production  of  the  water  droplets  in  oil  at  the 
flow  focusing  orifice,  and  a  snapshot  from  a  fast  video  recording  of  the  droplet  collection 
in  the  chamber  are  given  in  Fig.  6.2.  The  experiments  were  undertaken  in  a  continuous 
flow  through  mode  in  which  the  droplet  train  flows  continually  into  the  chamber,  and  the 
droplets  exit  through  the  port  on  the  opposite  side  of  the  chamber.  Because  the  width  of 
the  chamber  is  much  larger  than  the  width  of  the  microchannel  which  feeds  the  droplets  into 
the  chamber,  the  flow  in  the  chamber  is  relatively  slow,  and  the  droplets  only  slowly  change 
their  configuration  with  respect  to  one  another. 


6.3  Theoretical 

To  predict  electrocoalescence  between  the  water  droplets  in  the  two-dimensional  configura¬ 
tion  upon  application  of  the  electric  field,  the  electrostatic  forces  between  the  droplet  pairs 
must  be  calculated.  We  denote  by  £  the  magnitude  of  the  uniform  electric  field  applied 
to  the  droplets  in  the  chamber,  by  a  the  droplet  radius,  by  eon  and  e water  the  dielectric 
constants  of  the  oil  and  water  phases,  respectively,  by  a  on  and  crwater  the  conductivities  of 
the  oil  and  water  phases,  respectively,  and  by  7  the  interfacial  tension  of  the  droplets.  The 
magnitude  of  the  applied  electric  field  is  less  than  that  required  to  distort  the  shapes  of  the 
droplets  (the  electric  Bond  number,  Bo  =  is  less  than  one),  therefore  the  droplets 

are  modeled  as  monodisperse  spheres  with  known  diameter  and  position.  The  period  of  the 
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Figure  6.3:  (a)  Two  dimensional  configuration  of  water  droplets  in  crude  oil,  (b)  the  associated  simulation 
cell  and  overlay  of  color  map  of  the  magnitude  of  the  electric  field  and  (c)  Detail  of  the  electric  found  in  the 
boxed  area  of  (a)  and  (b). 


applied  field  is  much  greater  than  the  polarization  time  of  the  droplets,  r  =  ewater/^Water- 
This  produces  an  electric  field  distribution  around  the  droplets  that  is  quasi-static  and  can 
be  calculated  from  the  solution  to  the  Laplace  equation  for  the  electric  potential  Von  within 
the  oil  phase,  and  within  each  of  the  droplets  i  in  the  assembly  V^ater  : 

V2V0j,  =  0  (6.1) 

V2Vj,„(er  =  0,  i  =  1, 2...N  (6.2) 

where  N  is  the  number  of  droplets.  The  boundary  conditions  on  the  surface  of  the  droplets 
are  continuity  of  the  electric  potential  and  the  balance  of  the  electrical  conductive  flux  (no 
charge  accumulation  at  the  surface  at  steady  state), 

V(Dil  =  ^Wetter-*  *  =  lj  2...IV  (6-3) 

n*  '  {aOiNVOil  ~  VWater'Wwater}  =0  *  =  1,  2. ..TV  (6.4) 

where  n'  is  the  unit  normal  to  droplet  i. 

Although  the  external  field  is  uniform,  the  field  around  the  droplets  is  non-uniform  due 
to  their  polarization,  with  larger  field  strengths  at  the  two  poles  aligned  with  the  external 
field  and  lower  field  strengths  at  the  two  poles  aligned  perpendicular  to  the  external  field. 
In  the  case  of  a  single  droplet,  the  magnitudes  of  the  field  strengths  at  these  poles  will  be 
equal,  resulting  in  no  net  force  on  the  droplet  and  therefore  no  motion.  However,  droplets 
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in  proximity  to  one  another  will  experience  non-symmetric  local  fields  and  will  consequently 
experience  net  attractive  or  repulsive  forces  depending  on  their  orientation  relative  to  the 
external  field.  In  the  simple  case  of  a  single  droplet  pair,  this  net  force  has  been  theoretically 
calculated  using  bispherical  coordinates.  In  the  case  of  a  configuration  of  multiple  droplets 
as  studied  here,  finite  element  simulations  are  used  to  solve  the  above  set  of  equations  using 
the  COMSOL  Multiphysics  4.2  simulation  package  for  an  instantaneous  configuration  of  the 
droplets.  To  illustrate  in  Fig.  6.3(a)  is  shown  a  configuration  of  droplets  in  the  assembly  as 
captured  from  the  high  speed  recording.  This  image  is  used  to  construct  a  simulation  domain 
of  the  droplets  in  the  chamber  in  which  spheres  are  constructed  at  the  locations  of  the  droplets 
as  dictated  by  the  image  to  form  a  two  dimensional  assembly.  (Fig.  6.3(b)).  The  simulation 
domain  is  closed  by  top  and  bottom  walls,  separated  a  distance  h  apart,  corresponding  to 
the  height  of  100  pm  of  the  microfluidic  chamber.  These  walls  are  modeled  as  conductors 
with  surface  charge  to  apply  the  electric  field.  The  three  dimensional  simulation  domain  is 
then  closed  by  side  walls  in  which  the  potential  is  assumed  to  be  linear,  corresponding  to 
the  applied  field  £  applied  across  the  electrodes. 

Superimposed  on  the  simulation  domain  of  Figure  6.3(b)  is  a  color  map  of  the  magnitude 
of  the  electric  field  within  each  droplet  and  in  the  continuous  oil  phase.  The  large  electric 
field  strengths  between  the  droplets  in  the  chamber  which  align  along  the  field  direction  are 
evident,  and  indicative  of  the  strong  dipolar  attractive  force  which  causes  the  particles  to 
coalesce.  Note  that  the  field  between  the  droplets  which  are  aligned  perpendicular  to  the 
field  do  not  develop  the  high  electric  field  strengths  and  therefore  do  not  coalesce. 

This  electric  field  computed  by  the  COMSOL  simulation  is  used  to  calculate  the  net 
force  on  each  droplet.  In  this  quasistatic  formulation,  the  force  on  the  droplet  i  is  composed 
of  the  electrical  Maxwell  stress  tensor  force  on  the  droplets,  and  the  pressure. 

F Droplet  =  j n'-T'oudA  -  j)  '[/ou-ndA  (6.5) 

where  n1  is  the  surface  normal  and  TV>,/  is  the  Maxwell  stress  tensor  and  p'ou  is  the  pressure 
exerted  on  droplet  i  by  the  oil  phase.  The  general  expression  for  the  Maxwell  stress  tensor 
in  the  oil  phase,  where  Eon  is  the  electric  field  vector  in  the  oil  phase  is  given  by: 

Ton  =  cou^onEou  —  -^\Eou\fi  (6.6) 
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(a)  (b) 


Figure  6.4:  (a)  Forces  exerted  on  two  droplets  (in  red),  and  the  position  vectors  of  the  droplets  and  (b) 
details  of  computing  the  algorithm  for  identifying  nearest  neighbors  for  coalescence. 


where  5  is  the  identity  matrix.  Applying  the  Maxwell  equations  \7-eEou  =  0  and  VxEoi/  =  0 
and  a  vector  identity  simplifies  the  divergence  of  T on,  which  is  equal  to  the  pressure  gradient: 

V-T  on  =  £°^~£oV|E  ou\2  =  Vpou  (6.7) 

At  the  surface  of  the  droplet,  E on  has  no  tangential  component,  therefore  Eo,; | 2  =  Enm2 
and  the  net  force  on  the  droplet  can  be  calculated: 

?'DrV,e,  =  f'-fK^dA  (6.8) 

From  the  finite  element  simulations  of  the  electric  field,  the  force  on  each  of  the  droplets  is 
computed  from  eq.  6.8. 

This  net  force  is  expected  to  predict  the  droplet  motion  upon  application  of  the  electric 
field.  However,  the  close  proximity  of  the  droplets  prior  to  the  application  of  the  field  means 
that  no  significant  motion  between  droplets  is  possible  prior  to  coalescence.  Therefore,  the 
net  force  on  each  droplet  must  be  deconstructed  to  identify  the  component  of  the  force 
due  to  each  adjacent  droplet.  In  other  words,  the  forces  on  the  droplet  must  be  calculated 
on  a  pairwise  basis  to  identify  which  droplet  pairs  experience  net  attractive  forces  and  are 
therefore  likely  to  coalesce.  This  is  done  by  simple  vector  projection  using  the  net  force  of 
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Figure  6.5:  (a)  and  (b)  Two  configurations  of  the  droplets  in  the  chamber  and  (c)  plot  of  pairwise  electo- 
static  force  between  droplet  pairs  versus  their  separation  distance  s  normalized  by  the  droplet  radius  a  for 
these  two  configurations. 
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each  droplet  and  the  center  to  center  vector  between  each  droplet  pair  (Fig.  6.4(a)): 


F  ■  = 

pair 


[F2  —  Fi]-[r2  —  1*1] 


r2 


(6.9) 


With  this  definition,  negative  values  of  the  pairwise  force  correspond  to  net  attraction 
between  the  droplets  (along  their  line  of  centers),  and  positive  values  to  net  repulsion.  Be¬ 
cause  each  droplet  only  has  a  certain  number  of  nearest  neighbors  which  block  its  potential 
coalescence  with  other  droplets,  an  algorithm  is  employed  to  identify  which  droplet  pairs 
are  not  obstructed  by  third  droplets  and  are  therefore  potentially  able  to  coalescence  (Fig. 
6.4(b)).  Thus  for  example  droplet  2  is  an  accessible  neighbor  of  droplet  1  if  the  distance 
d  from  droplet  3  center  to  the  line  of  centers  from  droplet  1  to  droplet  2  is  greater  than 
the  radius  of  droplet  3.  Using  these  calculations,  the  resulting  pairwise  droplet  forces  are 
plotted  as  a  function  of  the  separation  distance,  s  (edge  to  edge)  between  the  droplets  in 
the  pair  (normalized  by  the  droplet  radius  a)  for  two  (instantaneous)  configurations  as  given 
by  Fig.  6.5(c)  for  two  configurations  as  given  by  Fig.  6.5(a)  and  6.5(b).  The  electric  held 
is  computed  first  in  the  oil  phase  from  these  configurations,  and  then  eq.  6.8  is  used  to 
calculate  the  force. 

To  evaluate  Fig.  6.5(c)  for  its  ability  to  predict  which  droplets  will  coalesce  from  thee 
mapping  of  the  pairwise  droplet-droplet  forces  for  the  configuration,  we  note  that  the  coales¬ 
cence  of  individual  droplet  pairs  proceeds  upon  application  of  the  external  electric  held.  Due 
to  the  fact  that  the  droplet  conhguration  changes  subsequent  to  the  conhguration  depicted, 
only  droplets  pairs  coalescing  within  the  hrst  0.1  seconds,  which  are  therefore  assumed  to  be 
coalescing  due  to  the  forces  calculated  from  this  droplet  conhguration,  are  considered.  These 
coalescing  pairs  are  identihed  on  the  plot  of  pairwise  droplet  forces  versus  separation  dis¬ 
tance  in  red.  As  expected,  droplet  pairs  with  large  negative  (attractive)  forces  between  them 
coalesce,  while  droplet  pairs  with  positive  (repulsive)  forces  between  them  do  not  coalesce. 
At  larger  separation  distances,  droplet  pairs  do  not  experience  large  forces  and  therefore  do 
not  coalesce  due  to  the  initial  droplet  conhguration.  As  shown  in  6.5(c),  several  droplet  pairs 
with  negative  forces  do  not  coalesce.  This  can  be  explained  by  the  fact  that  one  droplet  in 
the  pair  is  also  half  of  a  pair  with  a  larger  negative  force,  and  therefore  coalesces  with  the 
droplet  with  which  it  has  the  larger  attractive  force,  as  expected. 

In  conclusion,  we  have  demonstrated  accurate  prediction  of  pairwise  electrocoalescence 
of  water  droplets  emulsihed  in  a  crude  oil.  The  unique  ability  of  microhuidic  geometries  to 
generate  monodisperse  emulsion  droplets  and  observe  them  in  a  macroscopically  opaque  con¬ 
tinuous  phase  such  as  crude  oil  should  encourage  further  work  using  microhuidics  to  study 
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emulsion  stability. 
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